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ABSTRACT 



% 

I. 


Title of Dissertation: Experlaental a. i Theoretical Studies of 

Tnttt' stellar 0; ains 


Joseph Andrew Noth III 

Dissertation directed by: Dr. Isidore Adler, Professor 

Geochenistry Program, Chemistry Department 


Steady state vibrational populations of SIO and CO in dilute 
black body radiation fields have been calculated as a function of 
total pressure, kinetic temperature and chemical composition of the 
gas. Approximate calculations for polyatomic molecules have also been 
presented. Vibrational disequilibrium becomes increasingly 
significant as total pressure and radiation density decrease. Many 
regions of postulated grain formation are found to be far from thermal 
equilibrium before the onset of nucleation. Calculations based upon 
classical nucleation theory or equilibrium thermodynamics are expected 
to be of dubious value in such regions. 

Laboratory measurements of the extinction of snail iron and 
magnetite grains were made from ic*5nm to P^Onm and found to be 
consistent with predictions based upon published optical constants. 
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This implies thst nsll iron particles are not responsible for the 
220nn interstellar extinction feature. A feature which begins near 
160nm in the extinction spectrum of HDtblTP la identified as due to 
water ice. 

Measurements have been made of the critical partial pressure of 

SiO (P ) necessary to initiate avalanche nucleatlon in the SiO-H. 
c t 

s/stem as a function of temperature (750K<T<10P0IO , The condensate 
produced by this process is rathw than Analysis of P^ 

versus T using classical nucleatlon theory yields a value of ^00 

p 

ergs/cm for the surface free energy of the initial clusters. Despite 
the fact that this value is in reasonable agreement with those from 
the literature, numerous inconsistencies in the analysis are noted. 

It is shown that classical nucleatlon theory is not applicable to this 
system. 

Measurements of P^ vs T for the Mg-SiO-H^ system have been 
obtained (750K<T<inisIO . These are compared to similar measurements 
for the SiO-H^ system. The presence of magnesium lowers the 
nucleatlon barrier for T<c?^K but does not effect condensation at 
higher temperatures. Infrared spectra of both the initial condensate 
and samples annealed in vacuo at lOnOK are presented and compared with 
Infrared observations of 0H26.S+0.6 . Much of the previously ignored 
fine structure in both the 10 and 20 micron features could be modeled 
by laboratory produced amorphous silicate smokes. 
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Introduction 


Tho Qh«iloal oanpoaltion of inter atollor dudb his been the 
«d>J«)t of debate for lore then forty years. Mlthin the pert twenty 
five years however, our oonoeptim of the interstellar aediun has 
undergone a drertio change. This environaent was previously se«fi as a 
benign, quiet oasis in which solids could leisurely nucleate and grow. 
Once foraed, a grain was believed to persist for eons. 

Ve now know that the interstellar aediun is an extrenely hostile 
environaent. It is permeated by cosmic, X and Y rays; bathed in both 
hard and soft ultraviolet radiation; heated and nixed by shock waves 
from supernova explosions and intercloud collisions. Such processes 
bring about the destruction of many of the previously proposed grain 
materials much faster than the rates at which they could be reformed. 
In fact, because of the extremely low density of condensible material 
and the hostility of the environment, it is doubtful that any type of 
grain could form in the general interstellar medium. 

Our view of the chemical composition of interstellar solids, once 
believed to consist of ices of water, methane and ammonia, has 
undergone a similar transformation. Qrains must now be composed of 
materials capable of withstanding high particle and radiation fluxes 
for long periods of time. Metals, metal oxides, graphite, metal 
carbides and silicates all fulfill such criteria and are therefore the 
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prftforrtd Mt«rial8 to nodal intarstallnr dust. Navarthalass. even 
though lees are no longer expected to be the najor constitumta of 
interstellar grains, they night still be inportant in the nore 
sheltered r^ions of the interstellar nediun. 

This nodel of interstellar naterials is supported to sane extent 
by neasurenents of the depletion patterns of the elenents in the 
inter^ellar nediun (Field, 1974; Hrivedi and Lariner, 19RD* In such 
studies, the absolute absorption strength in a resonance line of an 
elenent is used to calculate its abundance along the line of sight to 
a specific star. These abundances are then compared to the 'cosmic' 
abundance of the element and its depletion factor is determined. 
Although there are numerous observational and theoretical difficulties 
in such determinations, the average cosmic depletion pattern indicates 
that refractory elements such as Si, Mg, Ni, Fe, Ca, A1 and Cr are 
more depleted than less refractory species such as S, 0, N, Na, Zn or 
P. Elements found to be depleted from the gas phase are most probably 
incorporated into interstellar grains. 

Although models based on the measured depletion patterns can 
predict the average elemental composition of interstellar solids, such 
models are not able to predict the chemical structure of such 
material. A complete model of interstellar grains requires an 
understanding of the processes by which they are formed and 
subsequently modified in typical astrophysical environments. Such a 
model necessarily requires that the structure of the major grain 
components be well characterized. 
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A auQo««8ful ■ttaok upon this pr<Al«n mift thonriTM** Inoorporatt 
at laaat thraa llnaa of raaaaroh. Fir at « it la naoaaawry to 
undar stand the physical envlronaent In uhleh nuoleatlon prootaaaa 
occur. Seomd, It requires the laboratory ^eduction of Material a 
which emulate the <^serv«d properties of interstellar grains. 

Finally, It is essential to understand the nuoleatlon process Itaelf , 
eapeolally for the refractory materials iSiloh are the most likely 
components of the dust. 

Chapter 1 brietly reviews the astronomical literature tihloh 
pertains to the formation, composition or structure of interstellar 
grains. Chapter 2 reviews previous laboratory measurements of 
materials proposed as models for the dust. This chapter also contains 
a brief outline of classical nucleation theory and the manner in which 
it has been applied to astrophysical systems. 

The results of theoretical calculations to determine the steady 
state populations of SiO and CO in dilute black body radiation fields 
as a function of the total pressure, chemical composition and kinetic 
temperature of the ambient gas are reported in Chapter 3> An 
approximate calculation applicable to polyatomic molecules is also 
described in this chapter. Chapter A describes the production of 
small iron and magnetite grains. Both of these materials are expected 
to be present in the Interstellar medium. Tne optical extinction 
spectra of these grains was studied from 1 95nm to 830 nm and compared 
to astronomical observations. 
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Chapter 5 deaorlbea •eaauraaant of the taaparatura d^mdanoa of 
tha oritieal praaaura of SIO naoaaaary to initlata apontanaous 
nuolaatlon In a pura SiO-H^ ayaten. ThaoC leaauranenta are than uaad 
aa tha baala upon whioh to analysa tha tenperatura dapandanea of tha 
onaat of avalanoha nuolaatlon in a Ng-SlOoH^ systam, daaorlbad in 
Chapter 6. Chapter 6 raporta atudiea of tha infrared apectra of 
anorphoua magnasiun ailieatas and as a function of annealing 

time at 1000K in vacuo. These materials ware produced under carefully 
controlled conditions in the laboratory. The implications of these 
studies for the formation of interstellar silicates in the atmospneres 
of oxygen rich stars are discussed, as are the observational 
consequences of this assumption. 


References 

Field. G.B. , 197*», Ap.J. , W. II53 

Trivedi, B.M.P. and Larimer, J.W. , 1961, Ap.J., 246, *563 
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Chapur 1 Aatroooaioal 0i>— rvtlooa 


Iht astronoaloal littratur* ralatlng to Intarstallar grains has 
Paan axtansivaly ravlawad on nunarous occasions. Tha raviaw of Savaga 
and Hathis (1979) lists aany of tha books and raviaw artiolas which 
hava appaarad bafora 1979. This chaptar will ba dlvidad into four 
saotions, in tha first saotion 1 will raviaw tha various spactral 
faaturas obsarvad in tha intarstallar axtinction curva froa tha far 
infrarad to tha axtrwa ultraviolat. In tha sacond 1 will raviaw tha 
inforaation axtractad froa scattaring and polarisation aaasuraaants . 

In tha third saction I will discuss tha typas of grain aodals which 
hava baan propoaad to axplaln thasa obsarvations. In tha last saotion 
I will raviaw tha information avallabla In tha litaratura on raglons 
consldarad to ba iaportant sources of grains. I will also discuss tha 
nachanlaas by which thasa grains could ba aodlflad or destroyed. 


I. Spectral Faaturas of the Intarstallar Extinction Curva 

One of tha oldest and aost useful sources of inforaation about 
tha properties of Intarstallar grains Is tha Intarstallar extinction 
curva. TSbla 1 lists values for tha average extinction of 
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Tabl« 1 

An Awnrift intnrttnUnr Eitinotion Curvt 


Color 

x(aloron») 

l*V«loron»"' ) 

EU-V)/E(B-V) 

A^/E(B-V) 


m 

0.00 

-3.10 

0.00 

L 

3.M 

0.29 

-2.94 

0.16 

K 

2.2 

0.45 

-2.72 

0.38 

J 

1.25 

0.80 

-2.23 

0.87 

1 

0.90 

1.11 

-1.60 

1.50 

R 

0.70 

1.43 

-0.78 

2.32 

V 

0.55 

1.82 

0.00 

3.10 

B 

0.44 

2.27 

1.00 

4.10 


0.40 

2.50 

1.30 

4.40 


0.344 

2.91 

1.80 

4.90 


0.274 

3.65 

3.10 

6.20 


0.250 

4.00 

4.19 

7.29 


0.240 

4.17 

4.90 

8.00 


0.230 

4.35 

5.n 

8.87 


0.219 

4.57 

6.57 

9.67 


0.210 

4.76 

6.23 

9.33 


0.200 

5.00 

5.52 

8.62 


0.190 

5.26 

4.90 

8.00 


0.180 

5.56 

4.65 

7.75 


0.170 

5.88 

4.77 

7.87 


0.160 

6.25 

5.02 

8.12 


0.149 

6.71 

5.05 

8.15 


0.139 

7.18 

5.39 

8.49 


0.125 

8.00 

6.55 

9.65 


0.118 

8.50 

7.45 

10.55 


0.111 

9.00 

8.45 

11.55 


0.105 

9.50 

9.80 

12.90 


0.100 

10.00 

11.30 

14.40 
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Int«r8tttllar aattrlil as a function of wavolongth. Ihia table la 
taken frcB the review of Savtie and Mathia (1979)* It auat be noted 
however that thla average curve ia by no aeana univeraally invariant. 
Indeed there are many atara which are known to poaaeaa aaoeoloua 
extinction (aee for inatanoe Sitko. Savage and Hiade, 1961 or Snow and 
Saab, I960). Some of theae curvea are dlaouaaed at greater length in 
Chapter 4 of thla work. 

Table 2 liata noat of the known apeotral featurea of the 
interatellar extinction curve, the material currently held to be 
reaponalble for the feature and the region or regiona in which it ia 
obaerved. Huch of thia information ia available in the reviewa of 
Huffman (1977) or Merrill (1979). In what followa I will briefly 
diacuaa wnat ia 'known* about each featire in thia table. 1 ahall 
begin with thoae in the infrared and proceed through the 
electromagnetic apectrin into the ultraviolet. 

a. infrared Extinction 

The far infrared region (x > 30tio) la the "laat frontier of 
aatronomical apectroacopy "(Ikifflaan, 1977). Obaervationa muat be 
taken at very high altitudea in order to avoid much of the 
interference due to atmoapheric water vapor. Becauae the field ia 
relatively young, moat of the apectra available in the literature are 
taken at only low reaolution. Theae indicate that the thermal 
emiaaion from duat graina aurrounding hot aourcea falla off rot«hly aa 
. Ihia impliea either 
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Tabl* 2 

Major Faaturaa of tha intaratallar Extlnotlon Curva 


Faatura 

Carrier 

Where Obaerved 

Inoraaaa in axtinotion 
(X < 130na) 

Silicates 

1 

Incraaaa in Extinction 
(X « 159MB) 

Water ice 

HD<MM79 

Extinction Haxiaua (220m) 

Graphite 

1.*! 

Change in Slope (x 410na) 

Large Particles 

1 

Diffusa Bands (Table 2b) 

(impurities ?) 

1 

Vary Broad Structure 

Magnetite 

1 

Hid->Infrarad Signatures 

( See Table 2a) 


45 micron ilbsorption 

Water Ice 

KL Nebula 

Far-Infrared Falloff 
(Intensity « x”^) 

'Dirty* Silicates 
Magnetite 

3 


1. General Interstellar Medium 
2* Molecular Qouds/Oompact Hll 

3. Oxygen Rich Stars 

4. Carbon Rich Stars 

Legend 

5. Planetary Nebulae 
Regions 6. Eruptive Stars 

7* Miso. Galactic Sources 
8. Galactic Nuclei 
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that tha gralna hava aMlaalvitiaa whloh ara highar than axpaotad la 
tha infrarad or that thay Mra niieh hattar Aaorbwa than axpaotad on 
tha baalB of a alltoata nodal • 9uoh anhanoad anlaalon night ha dm to 
tha praaanoa of Inpuritiaa within tha grain ( Aawiaatad,l975) or to tha 
praaanoa of * giant* gralna ( Rowan-gohinaon , 1976). It should ha notad 
that Day 0976* 1981) has produoad morphous sllioatas, tha axtinotion 
of whloh falls off as 

Both tha sansitivity and raaolution of far infTarad dataotors ia 
rapidly inproving. This has rasultad in tha raoant dataotion of a 
faatura at S5|in attrihutad to watar ioa (Eriolcson at al* 1980; 
Papoular at al, 1978). I should note howavar that this faatura has 
only haan dataotad in tha KL Nabula in Orion, tha ragioit )ith tha 
strongast known 3. 07 in ioa faatura. 

Hid infrarad (2|tn< x < 30m) spaotra of nany olassas of objaots 
hava baan availabla at aadim to high raaolution for about tha laat 
fiva yaara in tha abort wavalangth and of thia ranga (X < 13m)* High 
rasolution atudias at longar wavalai^tha ara publiahad nuoh laaa 
fraquantly. Thia haa raaultad in tha idautifioation of nmaroua 
faaturaa at X < 13m and alnoat nona abova (aaa Tbbla 2a). An 
axoaliant raviaw of tha known infrarad faaturaa fron an obaarvational 
point of viaw waa praaantad by Harrill (1979) st tha Vtorkahop on 
Thanaodynanioa and Kinatioa of Ouat Fbnaation in tha Bpaoa Nadiua. 

Thia waa augaontad by a raviaw of tha propartiaa of unidantifiad 
infrarad faaturaa (Uillnar at al, 1979) praaantad at tha saaa naating. 
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IMlt Zb 


Hid-Infrartd Sp«otral Ftaitrtt 


Band Gantar (Width) 
(■lorons) 

Saan 

• 

In 

Carriar 

Hhara Obaanrad 

3.1 

A 

Watar loa 

2 

3.3(0. 5)/3.4 

E 

( (^H 7 ) 

2. 5.7.8 

3.4 

A 

( C.M 7 ) 

(D.2 

3.5 

E 

(poly H 2 CO 7) 

3 

6.0 

A 

7 

2 

6. 2(0. 3) 

E 

(watar of hydration 7) 2, 5,7, 8 

6.8 

A 

7 

2 

7. 7(0. 8) 

E 

? 

2, 5.7.8 

8.6(0.3) 

E 

7 

2.5.7.8 

11.2(0.4) 

E 

7 

2. 5. 7.8 

9.7(3.0)/18(3.0) 

AE 

Silloatas 

1.2.3.(5).(6).7.8 

11.2(1.7) 

A£ 

Silicon Carblda 

4.(5) 

Contlnutn 

AE 

C» Ha tals. Natal 

Carbidas 1.4.5.6.7 

* A (Saan in Absorption) 

E (Saan 

In Bilsslon) 


Lagand 


1. Ganaral Intarstallar Madiua 

2. Holaeular Qouds/OMpaot HII Ragiona 

3. OKygan Rich Stars 

4. Carbon Rich Stars 


Planatary Nabulaa 

6. Eruptiva Stars 

7. Hiso. Galaotlo Souroas 

8. Galactic Nuclei 
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Tht aofit widtly h«ld btlUfa art that •■orphoua ailioataa ara 
raaponalhle tot tha faatyraa ohaarvatf at 9*7 aad 18 aloroaa, aillMn 
oarhlda la raaponalbla for tha faatura at 11.2 aloroaa aad taitar ioa 
is raaponalbla for tha faatura at 3.1 alormia. Viaaa ballafa ara by 
no aaana unlvaraal. lha aost eontrovarslal Idantlfloatlon of thla 
group la probably tha Idantlfloatlon of tha 9.7 and 18 aloroo faaturaa 
as dua to ailioataa. Hlllar and Dulay (1978) hava propoaad a aodal 
basad on oonblnatlons of dlatoalo oxldas auoh as S10« M|0« NIO and PsO 
which la plauslbla froa ' apaotrosooplc point of vlaw. Tha aodal la 
also quits intaraatlng in that It attwapta to pro¥lda, aiaultanaoualy. 
a aachaniaa for <tha salaotiva deplation of tha alaaanta and for tha 
proainant extinction faatura obaarvad naar 220na. It doas not attaapt 
to ax plain how such grains fora. This could be a major problaa since 
Day and Donn (1978) hava ahown that co-oondansation of hg and SiO 
produces amorphous silicates rather than separate oxides. This will 
be discu8sa<f in more detail in section III. 

An altarnativa hypothesis, which has received much more notoriety 
than acceptance , is that these features ara produced by organic 
materials (Hoyle and Uickramasinghe, 1977a,b). Huch of the recent 
controversy which surrounds this proposal has centered around tha 
identification of a feature at 3.4vim which could possibly bo due to 
the C-H stretching vibration in hydrocarbons ( (uickramasinghe and 
Alien, VjdO). Tne merits of this model will be discussed further in 
section 111. 
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k» Ml* IM M«ii by Insptotion of Ibblo 2b, thoro oro aaiy 
unosplotiiod footuroo Mhlon ooa bo •oaoolotod Mtth varlota groin 
populofclono. Oboorvationo indloato tbat tho foati^oa at 6.0 and 6*9 
•iorona oould bo duo to a alngit apaolta of gralna ahlla thoaa at 
3.1/3*<it 6.2. 7.7* 8.6 and 11.3 aiorona oould ariao in a aooond 
population (Harrillg 1979). BtMuao tho 6.0 Mid 6.8 nioron foaturoa 
art Mly obaorvod in oonpaot infrarod aourooa in aolooular olouda* tho 
■atonal roaponsiblo for thM la probably daatroyod in aora hoatilo 
onviroMonta. Thia haa lad to tha apooulation that thoaa faaturaa 
■Ight bo duo to to hydroearbona. Spaclfioaliy. tho 6.8 ■ioron foaturo 
oould ariao fron oarbonyl grot^ia in varioua ohoaioal onvironaonta 
whllo tho 6.0 aioron foaturo oould bo duo to tho bonding aodoa of CH^ 
(Willnor ot al. 1979). Tho aoohanlaa uhioh produooa tho aaaoolatad 
poaka at 3.4. 6.2. 7.7* 8.6 and 11.2 aiorona (whloh aro only obaorvod 
in aaiaaion) la atlli qulto eontrovoralal . Wilinor ot al (1979) 
bolievo that a fluoroaoonoo aaohaniaa oould aoocunt for tho 
obaorvationa but note that tho proooaa would nood to bo very 
offloiont. (For oxaaplo* oaoh UV photon in NGC 7027 would nood to 
produce 3 ZR photona.) Dwok ot al (1980) have ahown that thoaa 
faaturaa oould bo produood by tho thonaal aaiaaion of hot duat. The 
ooapoaition of the oonatituont or oonatituonta roaponaiblo for thoao 
foaturoa la atill unknown (Aitkon, 1980). 

Finally, thoro ia a foaturoioaa ooaponont obaorvod in tho 
infrared apactra of nuaaroua aourooa which aight bo duo to tho 
proaonco of aotallic graina, graphite, aaorphoua carbon or a aatal 
carbide. Thia ia oapoeially intoroating ainoo auoh a apoctrua is 
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oharaourittlo of that obatnrad in aovat (May and Haifitld, 1976{ Sato 
tt al, 1978) Mhort It la thomht that ooa atctet aotually ba obaarviaf 
train fonaation. Thia will ba furttoar diaouaaad in aaotion IV. 


b. Viaibla Eatinotim 

Tha viaibla portion of tha intaratallar axtinotion apaotrui ia 
oharaotarixad by a dapandanoa. Thia faot haa long baan known to 
indioata that tha aiaa of tha graina raaponaibla for tha viaibla 
axtinotion ia on tha ordar of tha wavalangth of viaibla light, ainoa 
tha axtinotion of largar graina would ba indapandant of wavalangth, 

•M 

whila that of anallar graina ahould vary aa x . Of Buoh aora 
intaraat ia tha praaanoa of at laaat 39 dlffuaa banda (Tabia 2b) whioh 
ooour batwaan 685nx and ttOni. A thorough r«viaw of tha obaarvationa 
uaa praaantad by Harbig (1975) at lAU Synpoaiua I 31 whila a brief 
overview can ba obtained froo tha work of Hartin (1976) or Mifflaan 
(1977). 

If thaaa faaturaa actually originate in graina - and thia ia 
atill a Matter of debate - than tha earriara nuat ba auoh aaallar than 
average ainoa typical aizad graina ahould ahow an aayBaatrio profile, 
bacauaa the atrangth of thaaa faaturaa doaa not aaao to ba correlated 
with the atrangth of tha far ultraviolet axtinotion but doea ahou sowa 
correlation with tha feature obaarvad at 220nn, tha banda nay ariaa in 
carbon graina. Naaauraaianta of the wavalangth dapandanoa of the 
polarixation have failed to detect any atruotura, thua atrangthaning 
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mi« 2b 

Hm OlffuM Interstellar Benda 


Waveleagthtmi) Central Oepth(S) Full Mldth at Half OeptMA) 


M 2. 80 

16.0 

20.0 

*150.18 

6.5 

3.0 

*172.60 

5.0 

5.0 

*175. t 9 

3.0 

5.6 

*176.30 

5.0 

5.3 

*177.97 

3.5 

1.8 

*188.20 

6.0 

17.0 

536.20 

2.5 

4.4 

5*10. *13 

5.5 

1.0 

5*12.00 

2.0 

11.0 

5*1*1.90 

4.5 

14.0 

5*18.731 

6.0 

4.4 

5**9.38 

5.0 

0.8 

553.50 

2.5 

23.0 

55*1. **6 

3.5 

0.8 

570.512 

7.5 

3.5 

577.83 

6.0 

17.0 

578.041 

37.0 

2.6 

579.496 

22.0 

1.3 

579.703 

22.0 

1.3 

584.41 

3.2 

4.5 

584.979 

9.5 

1.0 

601.09 

4.0 

4.2 

604.20 

2.0 

14.0 

611.30 

5.0 

0.85 

617.71 

7.0 

30.0 

619.595 

14.0 

0.70 

620.306 

16.0 

2.3 

620.649 

16.0 

2.3 

626.977 

17.0 

1.4 

628.391 

38.0 

3.0 

6 il .40 

4.0 

19.0 

635.35 

3.0 

3.1 

637.68 

5.7 

1.5 

637.93 

16.0 

0.86 

642.57 

4.0 

1.1 

659.74 

3.0 

0.6 

661.363 

34.0 

1.1 

666.071 

e e e 

a a e 
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o«M for graphitt oarrloro at laaat for Um bangs at SSJna, S78oi 
aag 628.3m (Hartln, 1978). Nora rtoant nork by UU, York and Sfiotf 
(t980) using obsarvatlona obtalnad froM tha ANS aatallita on a aanpLa 
of 110 hot stars has raaohad tha sma ganaral oonoluslona if ona doas 
not aooapt tha hypothasls that all mall grains ara rasponslbla for 
tha inoraasa in tha far ultraviolat aitinotion. in partioular it 
saaas quits poaaioia that nuoh of tha inoraasa in tha far ultraviolat 
axtinotion is oausad by mall siiioata partiolas uhila a aaparata 
population of anal! graph! ta partiolas oould axist uhioh night oarry 
tha diffusa bands (as surfaoa inpuritias ?). 

Anothar puxxling faatura of tha visibla axtinotion is tha 
ralativaly raoant diaoovary of tha \QS (Vary Broad Struotira). Tho 
VBS is batwaan 50 and lOO na wida, paaks around A56m and shorn a 
nininuB naar 570nn (Hayas and Bax, 1977; Hayaa at al, 1973; Bsx* 

197t). Basad on l^oratory spaotra of a oolloidal nagnatita 
suspansion and on spaotra of nagnatita froa tha orguail Hataorita, 
both of kdiioh ara quits sinilar to tha VBS, Hufflaan (1977) has 
oonoludad that nagnatita grains oould ba prasant in tha intarstallar 
nadiun, Nora raoantly, van Brads and Whittet Cl9d1) hava shown that 
tha VBS appaara to ba a widaapraad faatura of tha ganaral intarstallar 
extinction. Horaovar, thay find that a ohanga in tha shape of these 
features oould ba indioativa of a ohenioal changa on tha surfaca of 
the grains* 

The ranaining faatura in the visibla portion of tha axtinotion 
Qurva is a 'knee* or ohanga in slope at about tOOm. This oould ba 
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tfut to a population of iMrgar gralna (Hufltian. 1977)* Otnor 
oharaotariatioB of thia grain population raatin unknown and of litUa 
apparant intaraat. 


o. ULtraviolat Extinction 

Tha ultraviolat axtinotion apaotrui of aatarial in tha lina of 
aignt to Boat stara atudiad to data ahowa a paak at 217«5na ♦ 2.5na* 
Tnls la alaoat univaraally attributad to tha plaaaa raaonanoa paak in 
aaall graphita partiolaa ( approx iaataly 20 na in radiua). Huoh of tha 
raoant diaouaaion haa oentarad around tha poaaibla affect of an ice 
coating on tha obaarvable proper tiaa of auoh partiolaa (aaa Hacht, 
1981). Oulay. Millar and Williaaa (1979) howavar have praaantad an 
altarnativa aodal to produoa thia feature froa a oixture of diatoaio 
oxide gralna. Thlr Model will be disouaaed in aeotion 111. 

Within the laat year or ao, apeotra have been obtained which 
either do not ahow the peak near 220m (Sitko. Savage and Naade, 1981; 
Snow and Seab, 1980) or which ahow a paak that ia unuaually weak 
(Saab. Snow and Joaeph, 1981; Witt, fiohlin and Stacker, 1981; Witt and 
Cottrell, 1980). Moat regiona which exhibit anoaaloua behavior near 
220na appear to be denae cloud a where the graphite graina reaponaible 
for the featire Bight either be aaleotively daatroyad by aoBe unknown 
BeohaniaB or accrete nantlea to naak it. Thia problem ia diaouaaed in 
acre detail in Chapter S where we ahow that Magnetite graina could 
account for Many of the obaerved faaturea of auoh regiona. We alao 
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•how that in ioast om oaaa, that of ND<MM79» a auhatantlal lot 
•antla ovor a oora of Bagnatita graina oan raproOuoa tha puhliahaO 
ultraalolat apaotrua. 10 tha haat of our hnoidagga thla la tha flrat 
raportad Idanti float ion of tha watar ioa ultraviolat ahaorption adga 
in intaratailar natariaX. Thia ia tha faatura Xiatad at 159ni in 
TabXa 2. 

FoXXowing tha 220 no paak, tha avaraga intarstaXXar axtinotion 
daoraaaaa and raaohaa a nininui hatwaan igonn and 155nn. Thia is 
foXXowtd by an inoraaaa in tha far uXtravioXat axtinotion which 
oontinuas to at Xaast lOOn. k distribution of aaaXX siXioata grains 
is ganaraXXy thought to causa this inoraasa (Hathis, itunpXa and 
Nordsiack, 1977) since it has baan shown that tha strength of tha 
220na feature is not well correlated with it (UU, York and Snow, 1980) 
but that the strength of the 9*7 Micron silicate feature is (Savage 
and Hathis, 1979) 


II. Data froM Scattering and Poli»riaation studies 

Independent neasurenent of the scattering properties of 
Interstellar grains, if the total extinction cross section is already 
known, oan put quite severe liaits on the oonposition and size 
distribution of the grain population. Such Measureaents have been 
stteaptad for the dust in reflection nebulae and the dust responsible 
for the Diffuse Galactic Light (DGL). Unfortunately, such data is 


•xtr«i«ly diffloult to intorprtt dut eithor to unoortalnty in tht 

gtoiotry of ttM aouroo (or stallnr radiation fiald), unoartainty in 

tha dust distribution within tha nabula (or within tha galaxy) and 

unoar taint las inharant in tha thaory of light soattarii^. This last 

is baoausa tha pnasa function usad to intar jarat most obaarvations 

(Hanyay and Graanstain, 1941) dapands on only two paraaatars; a, tha ^ 

albado of the partiela and g, tha avaraga value of tha cosine of tha 

scattering angle. Within certain liaits* alaost any value of one ' 

paranetar can give a satisfactory fit to the observations if a 

suitable adjustment of the other paranetar is made. This is further 

oonpllcated by doubts concerning the applicability of this function 

based both on laboratory maasuremeats (Zarull, 1976; Chylek, Grams and 

Pinnick, 1976) and on Hie theory (see Savage and Hathis. 1979). 

Linear polarisation measurements of dust in oircumstellar shells, 
dense clouds and the interstellar medium yields information on the 
size distributions of these grain populations as well as some 
constraints on the composition of the material which causes the 
polarisation. Heasurements of the circular polarization yield more 
certain compositional constraints on the carrier of the polarization 
as well as information about the probable changes in grain alignment 
along the line of sight. We can thus study either the variation of 
the galactic magnetic field or the meuhanism of grain alignment 
provided that a suitable model for the other is assvjsed. Polarization 
measurements across individual spectral features can be used to 
determine if the feature is contained in the same population of grains 
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rtsponalbl* for the ovorall polar iutlon and thua oan halp to furthar 
dafina and dallnaata apaoiflo grain popuiationa. 

Tha following aaotiona will diaouaa aaaauraflianta of tha Diffuaa 
Galaobio Light, rafleotion nabulaa, linaar and oiroular polarization 
and aaaauraiiants of tha polarisation across individual spaotral 
faaturas. (toh of this inforaation is takan froa tha raviaws of 
Huffhan (1977) and Savaga and Hathis (1979) although a<Mra racant 
obsarvations hava baan addad where these either contradict or extend 
previous interpretations. 


a. Scattering Heasureaents 

There are aany probiens Inherent in the aeasureaent of the OGL. 

In the visible, the DGL is faint coapared to radiation due to direct 
starlight and to the Zodiacal light (scattering of sunlight from 
interplanetary grains - see lAU Syaposiua #90). At 530m, the average 
value of the ratio of the OGL to starlight to Zodiacal Light is one to 
three to six respectively (Roach and Megill, 1961). The ratio of DGL 
to starlight increases to about 1 to 2.5 near the galactic equator 
(Witt. 1968). Unfortunately, there is an uncertainty of at least ?0t 
in the brightness of the average stellar radiation field. The 
uncertainty of the aagnitude of the Zodiacal Light is also large when 
coapared to the intensity of the DGL. Uncertainties in the stellar 
radiation field increase at increased galactic latitudes (Savage and 
Hathis, 1979). 
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In Ura ultravioltt, problaas with th« Zodiioal Liflit 

•r« 1«M Important dm to th« rtlativt waaknaaa of tha aolar apaotrua. 
Unfortunataly, aodala of tha atallar ultravlolat radiation fiald ara 
■uoh laaa raliabia than thoaa in tha viaihla. Furtharaora, alnoa aoat 
of tha ultravlolat flux la provldad by fairly bright atara tdiioh hava 
a rathar Irragular dlatrlbutlon within tha galaxy (Hanry, 1977) auoh 
■odela ara aora diffloult to oonatruot. 

Daapita auoh dlffioultiaa, aoia uaaful Infonaation oan ba 
axtraotad fraa aaaauraaents of tha Intanaity of tha UL as a function 
of galactic latitude, in tha visible, the grains aust ba relatively 
efficient scattarars (a > 0.4); the aost probable value for a is 0.7 i 
0.1 (Savage and Hathis, 1979). XTie value of g is lass certain (and 
cannot be datarainad indepanaantly of a; tha aost probable value of g 
in tha visible is 0.7 i 0.2 (Savage and Hathis, 1979)* in tha 
ultraviolet, Lillie and Witt (19/6) found that 0.6 < i <, 0.9 over the 
spectral range 150na < A < 420na. They also concluded that a(l>300na) 
s 0.7 i 1, a(^220m) * 0.33 ± 0.05 and that a(153mi) > 0.6 ± 0.05. 

This indicates first, that the grains ara strongly forward scattering, 
second, that the extinction peak observed at 220m is at least 
partially due to absorption and third, that tha albedo of interstellar 
grains increases in tha far ultraviolet. Huffs an (1977) finds this 
last conclusion rathar distiarbing since it indicates increased 
scattering efficiency oooibined with decreased absorption at far 
ultraviolet wavelengths. Host solids show increased absorption at 
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tlMM •mrgi«a. Hmm tmuIU havt hoiMvar b««n oonfiraiadi by Hinry 
(1981) Nlw finds a > 0.5 and g > 0.7 In tba far ultravlolat. 

Navarthalaaa* raoant atudiaa of btoa aurfaoa brlghtnaas of tha 
Haropa Nabula (Andriasaa at al. 1977) Indioata that g • 0.25 in tha 
far ultravlolat (Witt, 1977). Siallar oonelualona wara draun froi 
obaarvationa of tha Orion Raflaotlon Mabuloaity (Witt and Lillia, 
1978). Baoauaa it ia hard to find a singla population of grains uhioh 
will hava a forward diraotad (g > 0 7 ) phaaa function in tha viaibla 
and ainultanaoualy produca aora iaotropio aoattaring in tha far 
ultravlolat, Witt (1979) has proposad that thara ia a saoond 
population of grains which cause tha far ultravlolat aoattaring. 

Thasa particles should be extrcaely aaall ao that they scatter 
essentially as Rayleigh particles. Because of tha observed incraasa 
in tha far ultraviolet extinction, the nuabar density of these grains 
would need to be extreaely high coaparad to tha density of the grain 
population responsible for the visible extinction. 

The aost serious problaas associated with the study of reflection 
nebulae are related to the detenaination of the gaoaetry of the systaa 
(Jura, 1977, 1979). Such problaas include the possible variation in 
dust density and grain sixa throughout tha nabula, as well as tha aore 
faailiar problaas associated with the exact dateralnation of the 
relative positions of the cloud, observer and illuiinating star. An 
additional problea in such studies is tha possibility that tha dust in 
thasa nebulae is not reprasentativa of tha general population of 
interstellar grains. Therefore, although there is soae evidence that 
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tht far ultravlolat ptiaat function of intorataliar grains appears to 
ba sort iaotropio than in the visible, this evidence is by no aeans 
conclusive. Because of the inherent uncertainties in the value of g, 
it is isposaible to reach any definite conolusions about the real 
value of a in the far ultraviolet. 


b. Polarization Heasur aments 

The variation of linear polarization as a function of galactic 
coordinates has been well described in the astronomical literature 
(see Mathewson and Ford. 1970; Coyne, (^hrels and Serkowski , 1974; 
Serkowski, Mathewson and Ford, 1975). Variations in the position 
angle of the polarization can be explained by variations in the 
galactic magnetic field and in the aligment efficier^.y of the 
polarizing material. Variation of the strength of polarization as a 
function of wavelength can only be explained with respect to the 
optical properties of the grains themselves. The wavelength at which 
the strength of the polarization reaches a maximun (imax) varies from 
about 410nm in the Cygnus OB-2 association up to about SOOnm in dense 
regions such as Orion, Scorpius and (^hiuchus (Savage and Mathis, 
1979). This is commonly attributed to variation in the average size 
of the grains. Denser regions such as Orion might tend to have larger 
grains due to the accretion of mantles or the effects of coagulation 
(see section IV) and thus Imax would tend to occur at longer 
wavelengths. A region such as the Cygnus OB-2 Association could be 
characterized by a smaller average grain size if many of the larger 
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grains hava bean aliainatad dua to ttia axtraaaly high radiativa flux 
of tha raglon. ilia avar%a valua of imn is approx laataly $50na 
(Savaga and Mi this, 1979)* 

Daspita tha rather large variations in both tha obaarvad strength 
and wavelength dapandenca of interstellar polarization, tha nonsalized 
polarization (P(l)/Paax) can be fit by a function which, for aost 
regions, depends solely on the ratio of l/Xmax (Equation 1) where K 

PU)/Poax s expi-K in^U/imax)} (1) 

usually is set equal to 1.15. A more recent study by Milking et al. 
(1930) which measured the wavelength dependence of the polarization 
from 300m to 1000m indicates that K should actually depend on xmax 
i.e. Ksl.TXmax, due to the observed narrowing of the normalized 
polarization curve in regions characterized by large values of xuax. 
This is thought to reflect a change in the size distribution of the 
grain populations of such regions rather than in the optical constants 
of the average grain material. 

Equation 1 (with K s 1.15) adequately describes most observations 
of regions of mall and large average grain sizes. Iherefore, the 
optical constants of the material responsible for the polarization 
must remain relatively constant frm the blue into the near infrared 
(Martin, 1974). This appears to rule out graphite as the carrier of 
polarization but allows such dielectric materials as ices, silicates 
and silicon carbide. 
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in th« n«ar infrarad U < I600i»)* aquation (1) Ooaa not 
adaquataXy daaerlba ttia dagraa of polar! aation obaarvad toward young 
■olaoular oloud aouroaa. Qulta oftan tha ■aiimai polar! lation at 
thaaa wavalangtha la nuoh highar than axpaotad. In fact, larga 
variationa ara obaarvad within tha aana aolaoular oloud ooaplax. Thia 
Indioataa that tha axoaaa polarisation ariaas oloaa to tha infrarad 
aoiaroa rathar than in tha intarvaning intaratallar aadita. Dyok and 
Lonadala (1981) faal that the noat plausibla axplanation for this 
phanonanon la that tha axoesa polarisation ariaas froa within tha 
cloud Itsalf by radiation which paaaas through a aadiun of aligned 
grains. 

In such sources it is iaposslble to aaasure P(x) at shorter 
wavelengths due to the heavily reddened nature of the oloud. However, 
neasurenents of the angle of polarization within the source can be 
oowpared to the angle aeasured in the visible for nearby OB stars. If 
the angles are correlated, then the excess polarization of the cloud 
is due to alignsent of the grains by the galactic aagnetic field. The 
strength of the field within these clouds however is expected to be 
auch greater than average due to the coapression of the aagnetic field 
lines by the collapsing cloud. Such a correlation has been found 
(Dyck and Lonsdale, 1979) although the evidence is certainly not 
overwlielaing. Tne observations indicate that the grains in both 
sources are aligned by the saae aechanisa - the action of the galactic 
aagnetic field. A oorrolary is that trapped aagnetic field lines can 
play a significant role in the evolution of such dense regions. 


In tht ¥Ulbl«, ■MturMMntt of tiM oirouLnr polar lutloa bava 
boon uaaO to plaoa aoM balpful eonstralnia on tba ocapoaitloo of tba 
polarising aatarial. It la vtil kaoiai ttwt tbt aavalai^th dapandaooa 
of olroiilar polarisation for a dlalootrlo obangaa algn at aoaa 
navalangth, lo. If tha laaglnary part of tha optical oonatant (k) la 
naar 0, than ic • laas, tha wavalangth at whloh tha llnaar 
polarisation raaohaa a aaslsiai (Nartln, 19B0). Aa k Inoraaaaa, xo 
baoooaa largar than xaax. Baoauaa tha Indas of rafraotlon of aatala 
la indapandant of wavalangth, olroulM* polarisation produoad by 
aatalllo partlolaa doaa not naoaaaarlly ohanga algn; althoogh If It 
doaa than xo « Xaax (Savaga and Mathla, 1979 >• Baoauaa tha quantity 
(xo-xaax) la also qulta aanaltlva to tha valua of n. It la poaaxbla 
for Xo to ba approxiaataly aqual to xaax avwi If kdO, If n la alao 
wavalangth dapandant. Thia aituatlon doaa ooour In tha oaaa of 
aagnatlta gralna (Shapiro, 1975). 

Obaarvatlona of tha olroular polarisation In nuiaroua aouroaa 
wara aada by Martin and Angal (1976) • In all oaaaa thay found that xo 
s max. Tharafore tha aatarial raaponaibla for tna polarisation in 
thaaa objaota aust ba a dialaotrio. Slailar atudlaa hava baan oarriad 
out by othar workara (aaa for inatanoa HoHillan and Tapia, 1977). 

Such atudiaa rula out graphita aa tha oarrlar of tha Intaratallar 
polarisation. 

Naaauraaart of tha wavalangth dapandanoa of the polarisation 
aoroaa apaoifio spaotral faaturaa can not only halp to clarify tha 
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int«rr«iationthlpt batwttn th« various faaturta put oan aiao yiald an 
uppar iiait to thair band atrangtha. Obaarvationai aaaauraBanta of 
tha wavalangth dapandanoa of tha polarisation aoroaa tba 10 aioron 
faatura in tna BN aouroa and in tha galaotio oantar hava shown that 
tha voluaa absorption ooaffioiant of tha aatarial rasponaibla for this 
faatura aust ba lasa than 3 x 10^ and 7 x 10^ raapaotivaly (Nartin* 
1978 p.202). This tanda to oonfim tha amorphous natura of thia 
malarial ainoa mora ordarad minarala (auoh aa quartz) tand to hava 
highar volvna absorption ooafficiants (Hartin, 1975* Capps* 1976; 

Capps and Knacka* 1976). 

On tha other hand* polarization maasuramanta across tha diffusa 
bands at a43nm* 5J6m and 628.3m indicate that tha grains which carry 
these features are not aligned (Martin and Angel* 1976). These 
features therefore do not arise In tha same population of grains 
responslole for the silicate feature at 10 microns. If the diffuse 
bands arise In the same population of snail graphite grains thought 
responsible for the faatura at 220m, than there should be no 
pronounced wavalangth dapandanoa across this region. Neasuremants 
have Indeed shown that the variation of PU)/Pmax across the 220m 
feature is adequately predicted by Equation 1 (Gahrals* 1974). In 
this manner* a systematic survey of the wavelength dependence of tae 
polarization across all observed absorption faatwas could ba 
axtaremly oanlflcial. It would be quite Interesting for Instance* if 
It ware discovered that there were two populations of diffuse bands* 
one of which arises in aligned grains* the other In nonallgned 
carriers. 
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in Hodtlt of tht Ch— loti CoapoUtioii oi lntcr»UUar Qr»tn« 


Th«rt ar« at laaat thrta dlffarant approaohaa froa Mhloh grain 
aodals hava baan oonstruetad. Dianiodjinanio nodais Hava baan propoaad 
to pradiot tha oaapoaition of tha graina fornad in tiM ataoaptiaraa of 
stara (Gilaar.. v'C9; Fix« I969««b« 1971) • oirouiatailar ahalla 
(Haekwall, 197D« tba praaolar naoula (Groaaian,l972; Lariaar and 
Groaanan, 1974) and tha ajaota froa auparnovaa (Lattiaar, Sohraaa and 
Groaaaan, 1978; Lattiaar and Groaaian, 1978). Caloulationa baaad on 
nuolaation thaory hava baan parforaad to iM'adiot tha rata of 
condensation and final aiza diatribution of tharnodynaiically 
pradiotad aatarlala in regions of astronoisioal interest. These aodais 
axaapiify a relatively pure theoretical approach. They start with 
basic thaory, naka sons sinplifying assimptions about tha ohraical 
coaposition, tanparatis'a and pressure history of the region oonsldarad 
and thereby calculate the properties of the condensate. To the extent 
that the theoretical framework employed and the assuaptions made 
correspond to reality, the conclusions drawn from such work are valid. 
The thermodynamic approach will be discussed in this section. 
Discussion of models based on nuoleation theory will be postponed 
until the theory is reviewed in Chapter 2. 

A second approach to the |»‘oolem ia purely empirical. In models 
of this sort, an attempt is made to synthesize the observed properties 
of interstellar grains by appropriate combinations of 'well known' 
matariais. Such models include the work of Nathis and co-workers 
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(Mitliis, lwipL« mad lk>r<l»l«clc, 1977 • lwr«aft«r MUI; Mathis, 1979; 
Mathis and Wallaohorst, 1981) who attwpt to obtain quantitative 
nawi MS to observations of the interstellar eatinotion in the range 
none < 1 < lOCOni using Mixtures of grsphite, enststite, olivine, 
silicon carbide, iron and nagnetite grains. Also included are the 
less quantitative Models of Duley and oo-uorkers (i.e. fiuley, Millar 
and Killians, ;979), Hoyle and Uickranasinghe (1977, 1979) and that of 
Hong and Greenberg (1980). Such Models will be described below. 

A third approach is to atteni^ the laboratory production of 
Materials sinilar to those expected to be fomed under astrophysical 
conditions. Such work has been indertaken by Donn and co-workers 
(i.e. Day and Donn, 1978; Donn et al, 1981) who have condensed an 
anorphous nagnesiun silicate fTon the vapor. The infl*ared spectrua of 
this Material resenbles the 10-20 Micron spectra of oxygen rich 
regions. Matrix isolation studies of the interned iates in this 
reaction have been carried out in an attenpt to understand the 
nechanisn by which such species condense (Khanna, Stranz and Donn, 
1980). More refined vapor phase experinents of the nucleation of 
SiO-H^ and Ms-SiO-H^ systens have recently been conpleted (see 
Chapters 5 and 6 of this dissertation) . The experinental approach has 
also been taken by Greenberg and colleagues (Hagen, Allanandola and 
Greenberg, 1979; Allasandola, Greenberg and Homan, 1979) who 
concentrate on the processing of ice Mixtures by ultraviolet 
radiation. Such Materials fora a host of organic coapounds which 
could possibly account for the unidentified infrared SMission features 
discussed earlier and might also help to explain the foraation of sone 


28 



of ttM ooro ooaplox orgmio aolooulos cJiitrvod in donot oloudt. Thla 
•pproMh Mll4 bt dlMuiMd in tiM Moond rovion ohnptor Mhioh will 
onphnalat iaborntory inwontlfatlona ralatad to tha atudy of graina. 


a. Thoniodynaiio Nodtla 

Tha tharnodynanio approach la by far tha noat atraightforward 
■aana of pradioting tha oonpoaitlon of Intaratailar natarlala, 
although a faw rathar inportant initial aaauaptiona nuat obvioualy ba 
•ada. Firat, ona nuat aaauna that tha ayataa undar atudy ia indaad in 
tharnodynanio aquilibriua. Saoond, ona nuat aaauna that tha 
thanaodynaaio data fOr all apaoiaa of inportanoa hava baan inoludad in 
tha oaloulaticns. Third, ona nuat aaauna a particular tanparatura - 
prasaura history for tha systan of intarast. Fourth, ona nuat aasiaa 
an initial oonposition appropriata to tha ragion undar oonsidaration . 
Finally, ona usually assuias that tha aquilibriun approx ination braaka 
down to sons dagraa in ordar that aoaa of tha aarliar oondansataa ara 
prasarvad through latar natanorphosas. Ona now aaraly bagins at an 
appropriately high tanparatura, so that all of tha spaoias ara 
initially in tha gas phasa, and allows tha tanparatura to fall until 
tha first solid oondansas. 

As an axanpla, lat us oonsidar tha oondansation of oorundua 
(Al^O^) lYoa tha vapor. Equation (2) ia a oonsaquanoa of tha usual 

Kaq • P^(AU « P^(0) (2) 
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definition of ttio oquilibritn oonstnnt for the reeotion of e gee end 
•olid. SiBiierly, equale exp(*eO/KT). It ie ueueliy esauied that 
the gaa ia well Mixed ao that by keeping track of the total preaaure 
of the MaJ(W apeoiea (uaually H atoMa) and defining a quantity i auoh 
that it repreaenta the fractional abundance of an elenent in the gaa 
phaae, equation (2) can be reiiritten aa 

expims - AH)/kT) « 

in which the neaaured themed ynaaio quant it iea &H and aS have been 
aubatituted for ac and N^kT haa been uaed for P„. N„ ia defined as 
the nuMber denaity of H atcaa» Mien the left hand aide of Equation 
ii) beconea equal to the right hand aide, Al^o^ grains are said to 
condense fran the vapor. 

When condensation occurs, one usually asaunes that tne reaction 
procedes to ocapletion, i.e. continues until one of the grain 
constituents beconea virtually depleted from the gas phase. The 
composition of the remaining vapor is then recalculated and the system 
is again allowed to cool until the next solid species becomes stable. 
This process ia continued either until the vapor is completely 
condensed or until one decides that the assusption of thermodynamic 
equilibriuB has broken down sufficiently that further computations are 
meaningless. It should be noted that these caloulationa can (and 
often do) include provisions for the reaction of previously condensed 
material with either gaseous or even other solid components within the 
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ayatM. Thua iron fraina, Mhion oondanat at a raXativtly high 
taaparatura, oouid ba eonvtrtad to •agnatita at lowar taaparaturaa, 
oorumtui oouid ba tranaforaad to anorthita, ato. 

Tharaobynaaio aodala oan ba aalb to hava had aaaa aaaaura of 
auooaaa. Gilman (1969) pradiotad that tha oondanaataa from oxygan 
rich M atara ahould contain an appraoiabla fraction of ailicataa. 
Siiicata-iika abaorption at 10 aiorona waa at^aaquantly obaarvad^ 

Hoyla and wiokramaaingha (1962) axaminad tha problam for carbon rich 
stars and pradiotad tha formation of graphita. This was subaaquantly 
obsarved in aarly ultraviolat rookat studiaa (Staohar, 1969) as tha 
prom inant faatura at 220nm. Grossman has had soma suooass in 
formulating a model based upon tha tharaodynamio approach which 
explains the high temperature mineral assemblages found in 
carbonaceous chondrites. Such models therefore serve as excellent 
starting points from which to study condensing systems since they 
provide a frmiework around which existing data can be organized and 
future experimentation planned. 

Unfortunately, some of the basic assuiptions of the thermodynamic 
approach have been questioned (Oonn, 1976) on the grounds that, in a 
system of approximately cosmic composition, kinetic considerations 
tend to preclude the formation of seme of the more complex species 
predicted by thermodynamics. Moreover, Blander and Katz (1967) have 
shown that differing degrees of vapor supersaturation are required 
before species predicted by thermodynamic models would nucleate, if 
indeed they actually were to form from the vapor at all. This would 
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disrupt tilt prsdic^tsd oondsnastios ssqusnos sad oould profoundly si tor 
tiM prodlotsd oosposition of spoeiss fonisd at lonsr twpsrsturss 
sinos tne ooaposition of tho vapor oould bo drastioally dlffaront. 

Hors raoantly, it has baan shown that tha vibratiwtal 
tanperaturas of nolaoulas in tha ralativaly low prassura anvironaants 
auggaatsd as souroaa of grain oondansation ai'a signifies itly lowar 
than tha gas kinatlo taaparatura (Nuth and Donn, 19d1; Chapter 3). It 
has praviously baen shown that tha tanparaturas of snail grains in 
suoh ragions oould also ba significantly lower than that of tha gas 
(Arrhenius and Da, 1973; Ds, 1978). It therefore seens highly 
unlikely that nost of the regions thought to produce grains (see 
section IV) are In themodynanlc equlllbrlui. 


b. Dtplrlcal Approaches 

Observational neasurenents of the extinction and polarization 
properties of Interstellar dust are available from the far Infrared to 
tne extrene ultraviolet (see section 1). If one nakes a few 
assuiptlons as to the types of grains which night fom In various 
astrophyslcal envlronnents then an attenpt can be nade to fit these 
observations by appropriate distributions of suoh grains. The nost 
conprehenslve and quantitative treatnent of this type has been carried 
out by Hathls and co-workers. Each of the other nodels discussed 
below attenpts to explain only a relatively restricted set of 
observations and even these are usually treated only qualitatively. 
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Hathia, iuapla and Nordalaolt (1977) - HM- uaad Mia tiMory to 
oalouiata ttoa abaoluta axtinotion affiolanoy aa a function of 
wavalangth far varloua alxturaa and alaa dlatributiona of propoaad 
grain aatariaia» lhaaa aiaa diatributiona wara oMiatrainad auoh that 
tha total ratio of an individual alaaant to hydrogen oould not axoaad 
ooaaio abundanoaa. lha aixturaa atudiad oonaiatad of graphite* 
olivine, anatatita, iron, ailioon oarbida and Magnetite taken 
individually and in ooMbination (up to three at a tiMa). lha 
oaloulatad extinction aa a function of wavelength waa then ooMpared to 
the interstellar extinction curve by neana of a least squares 
analysis. This analysis was weighted so as to force the fit to 
oonforai to the observed width and position of the 220na feature. 

They found that any Mixture which gave a satisfactory fit to the 
observed extinction had to contain graphite grains. However, graphite 
alone did not yield a good fit. The best fit for a binary Mixture 
contained graphite and olivine. This was followed by graphite plus 
(in order) enstatite, silicon carbide, aignetite and iron. Even the 
worst of tnese ccMbinations (graphite plus iron) however, was not poor 
enough to be excluded by the calculations. Given this fact, it is not 
surprising that the best ternary Mixture contained graphite, olivine 
and enstatite and that enstatite contributed less than 101 to the 
total extinction. 

HRN did not build in a preconceived site distribution, but rather 
allowed the total nueber of grains in a series of "bins" to vary 
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Indvptndvntl y until tht t>«st possiblt fit hnd b««n aohinvnd. Thay 
found howavar that tha final diatribution obtainad for any oonponant 
of a nixtura which gava an aooaptabla fit to tha obaanrationa oould ba 
daaoribad by a powar law diatribution of tha fora n(a) « a*^ whara 3.3 
< q < 3.6. It waa suggaatad that auoh a diatribution of grain airaa 
oould poaaibly raault froa atoohaatio prooaaaaa in tha ragion of grain 
formation. Thia waa aubaaquantly invaatigatad by Biaraann and Harwit 
(I960) who oonoluda that' thia ia indaad tha oasa. 

HRN alao found upper and lower size cutoffs for their 
distribution. For graphite, particle radii varied from about 5nn to 
near lOOOnn while radii of other components remained between 25na and 
250nm. These limits were not rigidly established however since larger 
particles of any composition contribute only 'gray' extinction, while 
smaller particles are in the Rayleigh regime so that the extinction 
per gram is virtually independent of particle size. 

Although HRN calculated the average polarization for their 
mixtures as a function of wavelength, tney found that xmax occurred at 
much too short a wavelength to agree with observation. Hathis (1979) 
later made a more thorough calculation of the wavelengtn dependence of 
the polarization. He assuoed that the sane grain mixture which 
accounted for the extinction measurements could also be used to 
explain the degree of polarization. CDservational evidence cited 
earlier (Gehrels, 1974) indicates that the material responsible for 
the 220nm feature (graphite in the HRN mixture) does not contribute to 
the polarization. Hathis (1979) therefore assumed that only 
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dUltotrio aaUrials suoh as all loon oarblds (n s 2.5) and sllioabes 
(n ■ 1.6) will ba iaportant. Ha arguas that a diatrlbution of 
infinita oyiindars, whoaa radii follow a powar law distribution of tha 
fora n(a) « will yiald tha aquivalant astinotion par unit 

length as a NRN distribution of sf^aras ( n(a) « a**^*^ ) of tha saaa 
aatarial. 

Using thase assinptions he shows that the above distribution of 
cylinders - roughly equivalent to a MRN diatrlbution of elongated 
particles - yields excellent agreement with Serkowski's Law (Eq. l). 
Furthermore, the variation of xmax can be accounted for by a 
reasonable variation of the upper and lower limits of the radii in the 
size distribution used in the calculations. In addition. Hathis 
(1979) finds that the smallest dielectric particles required to fit 
the extinction data are not likely to be aligned since this would tend 
to shift xmax too far into the ultraviolet. He postulates that, as 
larger particles grow, small iron or magnetite particles could be 
incorporated into the grain which would tend to increase the alignment 
efficiency of the larger particles. He therefore demonstrates that a 
HRN mixture not only provx.cs an excellent fit to the average observed 
extinction curve from llOnm < X < lOOOnm. but simultaneously provides 
an equally satisfactory explanation for the observed wavelength 
dependence of the polarization. 

In a later paper. Hathis and Wallenhorst (19d1) extend the miJeV 
into the near infrared and also attempt to fit ooservations of 
anomalous extinction seen towards nighly reddened regions. In the 
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first application t they show that the aaaa indioas of rafraotion and 
aiia distributions uaad to fit tha viaibla-ultravioiat obaarvationa 
also yiaXd satisfactory agraaaant with tha color ratio (I-L)/(B-V) for 
both tha avari^a intaratallar axtinotion and for that of p Ophiuohus 
(a highly raddanad rag ion ) . 

In tha aaoond application, prograssivaly paculiar axtinotions ara 
axMinad as a function of tha raddaning paraaatar (R^ i AyE(B*V)). 
It is shown that as this paraaatar inoraasas, tha raaultant axtinotion 
curves can ba explained by a continuous increase in particle size. It 
is denonstrated th~t first, small silicates are eliminated. As the 
density of the cloud increases, the upper size limit of the silicate 
population increases as well. At the highest values of R^ the small 
graphite component of the mixture begins to disappear. These 
systematic results could also be caused by an increasingly flatter 
particle size distribution as R^ increases, which first effects 
silicates and only begins to perturb the graphite distribution at the 
highest densities. It may be possible to distinguish between these 
alternatives by a careful investigation of the wavelength dependence 
of the polarization as a function of R^ using the techniques of Mathis 
( 1979 ). Such data could aid in the development of models of the 
coagulative growth of grains in denser astrophysical environments. 

Although the model of Mathis and coworkers is quite satisfactory, 
it is not the only model available in the literature. M shall now 
briefly review some of the alternatives. 


Dulty and oowrtetrs hav* advooatad a aodal of intaratallar aolids 
baaad upon tha proportlaa of diatoaio oxldaa (aaa Millar and Oulay, 
1978; Dulty, Millar and Hilliaaa, 1979; Millar and Dulty. 1980). In 
thtat aodtla, tht 10 aioron 'ailioatt* ftaturt ia produotd by SiO 
whilt btit ftaturt at 18 aiorona ia produotd by a ooabination of MgO 
and FeO. Dulty (1980a) olaias that variationa in tht rtfraotlvt indtx 
of the ailioon oxide atrita SiO^^ (1 < x £ 2) aight bt reaponalblt for 
the variation obaerved in xaax of denae olouda, rather than the 
uaually aaained variation in average grain aixe. In addition, Duley 
(1978) haa alao shown that nixed oxide graina ahould be either 
auperparanagnetio or weakly ferroaagnetio and therefore quite eaaily 
aligned by the Davia-Ck'eenatein aechaniai. 

Although Duley and coworkers have used this aodel to account for 
everything froa the interstellar cheaistry of sulfur (Duley, Millar 
and Uillians, I960) to the selective depletion of the eleaents (Duley 
and Millar, 1978; Duley, 1980b) to the origin of the diffuse 
interstellar bands (Duley and McCullough, 1977), few of their 
predictions or calculations are quantitative. Furtheraore, they 
alaost always fail to discuss inconsistencies between their aodel and 
observations. 

As an exaaple, their aodel predicts that the 10 and 20 aicron 
features are produced by the sase population of grains responsible for 
the feature at 220ns. Observational evidence cited earlier suggests 
that these featires arise in seperate grain populations. Their aodel 
also suggests that SiO, NgO and FeO will condense into aaall aixed 


37 


grains without altaration of tha apaotrosoopio propartiaa of tht 
individual oxMaa. Exptriaantal avidanoa (Day and Oonn, 1978) 
auggaata tha oppoaita, avan at oryoganio taaparaturaa (Khanna* Stranti 
and Donn, 1980). Hillar and Oulay (1980) disouaa tha faot that tha 
apaotral oharaotariatioa of altarnativa aodala for tha 10 and 20 
■ioron faaturaa - naaaly tha aaorphoua ailicata of Day and Oonn 
(1978)- Might ohanga at tha taaparaturaa obaarvad in oirouaatallar 
ahella (bIOCOK) yat thay fail to diaousa tha bahavior of thair 
proposed Mixture at auoh taaparaturas. As a final axaapla, Oulay and 
HcCullough (1977) have shown that benxana ions adsorbed on 
protosilicabas could ba rasponsibie for tha diffusa band at A43m. yat 
banzana has navar baan obaarvad in tha intarstallar aaditai. Indaad no 
aromatic hydrocarbon has avar bean obaarvad in any aatrophysicai 
source (Hyars, Thaddeus and Link, 1980; Mann and Williams, 1980). 

In summary, although many of tha ideas of this group, discussed 
in the aforamantionad papers, might have at least some merit, 
sufficient Justification is seldom given for tha assinptions used to 
derive their conclusions. This obviously detracts from the value of 
their work. 

An organic composition has been proposed by Hoyle and couorkars 
as a 'single component* modal of the interstellar dust. A modal based 
on the infrared spectrua of polysaccharides was first proposed to 
explain observations of the infrared features normally attributed to 
both silicates and to water ice (Hoyle and Wickramasinghe, 1977a; 
Hoyle, Olvasen and Wickramasinghe, 1978). These papers did not 
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attMipt to aaplaln tiou auoh a grain oouid ta produoad« nor did thay 
addrtaa tha obvioua problM of how auoh a grain waa to aurviva in tha 
raiativaly hoatila anviroiMant of tha intaratallar aadiui. Suoh a 
Modal pradiota a oorralation batwaan tha 3.07 and 10 aioron faaturaa 
ainoa both ara hypothaaixad to ariaa in tha aaaa population of graina. 
Thia oorralation ia not obaarvad. Although tiia faatura at 10 aiorona 
ia obaarvad in a wide variaty of aouroaa, tha faatura at 3.07 aiorona 
is only obaarvad in fairly danaa olouda (saa Table 2a). 

This aodal now saaas to have bean abandoned in favor of a 
biochaaical origin for intaratallar grains (Hoyle and Uiokraaasingha, 
1979). In this modal, tha faatura at 220na is attributed to 
biochaaical chrcnophoras (Hoyle and Wictaraaaslngha. 1977b, 1978; 
Wickraaasingha , Hoyle and Handy, 1977)* Bacteria ara held rasponaibla 
for tha visible extinction while viruses account for tha far 
ultraviolet observations (sea Fig. 3 of Hoyle and Wickranasingha , 

1979). In a more recent article (Hoyle and Wickraaasingha, 1981), the 
int'rared spectra of E.coli, yeast and E.coli annealed at 300''C are 
ccapared to observations of the unidentified 3.4 micron feature. 
Although the observed spectrus ia reproduced by such materials, the 
authors again fail to adequately discuss the origin and survi'/ytl of 
such 'grains'. It is interesting to note that the authors do discuss 
the fact that the cosmic ray bombardment of such species could produce 
free radicals which could then account for some additional infrared 
features. They fail to note however, that such bombardment should 
easily kill their grain popular icr. 
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To suMorlM, aXthomh aoat oorrospondonoo to tiM oOMrvod 
proportlto of int^ritoUar aatorlalo o«n bo roproduood by • nodol 
boatd OR organlo Mnotituonto, tho aodoX !• far froa aatiafaotory for 
tlM roaaona notad abovo. In addition, it ia knowi ttiat both oxygon 
and carbon rich atara oxpoll oondonaiblo atoaa and aoloouloa into tho 
intoratoUar aodiiai and that thoao oloaonta aro obaorvod to bo 
doplotod froa tho gaa phaao (MMO. A aodol Mhioh faiia to utilixo 
auoh aatoriala auat bo at loaat partially in orror. 

Groonborg haa boon aotivoly ongagod in tho atudy of intoratoUar 

aatoriala for aoro than tiMnty yoara. flooontly, Hong and Groonborg 

(1980) publlahod a unifiod aodol of intoratoUar graina which attoapta 

to oatabliah a quantitative rolationahlp botwoon grain also and 

allgnaont officioncy aubjoot to ooaaio abundanoo oonatrainta. Thoy 

fit tho obaorvational oxtinotion and polarization aoaauroaonta uaing a 

biaodai grain aizo dlatribution. Solativoly largo ooro^aantio graina 

aro noodod to fit tho viaiblo-infrarod oxtinotion and polarization 

data whilo aaaller, baro ailioato (O.OIaioron) and graphite 

(0.013aicron) graina account for tho ultraviolet obaorvationa. 

Although ainglo aizod coroa aro aaaunod, tho aizo dlatribution of the 

coro-aantlo graina ia oxpoctod to follow Equation (A), tdioro a, a^ and 

c 

n(a,a^) ■ n^expl-^Ua-a^l/a^)^) (4) 

n aro roapoctivoly, tho total particle radiua. tho core radiua and a 
normalization oonatant. Tho parameter a^ acta to out off tho aizo 
diatribution for moall particloa. 


40 


Ovtf'alit Umm aifttora do oot dttMpt to obtota roatUt* )£ilok opo 
M qumtlUtivt •• tlioM of Mitlild aod ooworlwrd oltod ibovd. 
and oraaobarg (19B0) ooaolyda tbat aoaa anftaiioaBaat of tAa 
Oavia-Graanatala grain ai*gaiaat •aobaniai ia aaoaaaw*y in ordar to 
adaquataly fit ^aarvationai oonatrainta. Hathia (1979) t on ttw otbar 
hand, finda that ha doaa not avan naad tha full d^raa of aligmant 
produoad hy thia aaana to aoourataiy nodal tha f»lariaation of a wida 
variaty of aouroaa. As Hong and (iraanbarg (1980) do not oonnant on 
any of tha findinga of Hathia and ooworkara, it ia hard to Judga tha 
validity of thair oonoluaiona. Thia nodal doaa pradiot that tha 
graina raaponaibla for tha viaibla axtinotion do hava fairly thick 
nantlaa of prinarily CHO oonpoaition. Tha voluna ratio of oora/nantla 
natariai ia approxinataly 0.1 in thia nodal. It aaana that avan 
highly irradiatad nantlaa ahould produoa aona infrarad faaturaa which 
would ba obaarvabla ovar vary long pathlangtha if auoh nantlaa 
actually account for aa nuoh aa 90S of tha naaa of intaratallar 
graina. Laboratory naaaurananta of tha quantitativa abaorption 
ooaffioianta of highly irradiatad ioaa ara naadad bafora auoh an 
obaarvational taat of thia nodal could ba parfomad. 


IV. Tha Fomation. Growth and Daatruction of Intaratallar Graina. 

In thia aaction I will raviaw tha obaarvational avidanoa ralating 
to tha fonaation, growth and daatruction of intaratallar graina. 1 
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will therefore first consider astrophysioal regions tdiere dust is 
observed to fora. Next I will discuss regions in which grain 
foraation has been hypothesixed. Finally« I will briefly review the 
processes thought responsible for grain growth, aantle foraation and 
grain destruction. The aechanian of grain foraation will not be 
discussed in this section but rather will be deferred mtil Chapter 2. 


a. Observed Regions of Oust Foraation 

Dust has been observed to fora in only two astrophysical 
enviroments, the envelopes of cool stars and the ejecta froa novae. 

1 should note however that there exists soae controversey in the 
interpretation of the observational data used to infer the 
condensation of grains in novae on a rapid tiaescale. Grain foraation 
in circuastellar shells is much less controversial. Several studies 
of this phenoBonon will be suamarized below which illustrate the range 
of stellar types, mass loss rates and observations reported in the 
literature. 

Circumstellar Condensation ; Hagen (197d) studied the c ire um stellar 
gas and dust shells of M giants and supergiants. She combined 
detailed radiative transfer models with high resolution optical 
spectra to derive the column densities of circuastellar material from 
observed line asymmetries. She also used the observed infrared energy 
distributioiic to derive the column density of circuastellar dust from 
the height of the 1U aticron silicate feature using a model which 


inoludta tha effaoU of Mlf abaorptlon by tha oust, from a 
ooapariaon of tha | s to Oust ratio aba oonoiudaO that aiaoat ail of 
tba aatala in tha shalla of stars of M6 or latar ara loo lead into 
grains. 

kithough tha innar radii of tha oircuastaiiar ahaiia ara not wall 
dataniinad, aha oonoiudad that thay Must bagin quite oiosa to tha 
photoaphara (within • 2R^>. Partioia danaitiaa ara oaiouiatad to ba 
in tha ranga 10^ to 10^ atoMa/oo (Praaaura > lO^^^ata). Haaa loaa 
ratea for tha atara in har aanpia ware in tha ranga 10 to 10 solar 
aassas par year. Although tha affect of radiation pressure on tha 
grains nuoieated in the atmoapheres of oool stars has been proposed as 
the nachaniaa which drives mass loss, Hagen (1978) concluded that 
grain formation is the consequence of mass loss rather than the 
reverse. It is unfortunate that gas temperatures within the shells 
could not be determined from her observations. It is worth noting 
however that the two stars in her sample which were observed to 
contain a significant metallic fraction in the gas phase - a Sco and a 
Her - both have hot companions. This would imply that the 
circuasteilar material in these stars is hotter and that either grain 
formation in such regions is inhibited or that grains formed in the 
shell of the cool star are evaporated by the hotter companion. 

Andriesse. Oonn and Viotti (1978) have studied the condensation 

of dust around the high luainosity star n Carinae. They derive a 

'-‘2 

stellar mass loss rate of 7.5 x 10 solar masses per year assuaing 
constant outflow since approximately 1800. Condensation begins when 


gat ooola froi photo aptiario taoparaturaa (• SOtOOOK) to thoaa of 

ordar ^tOOOK. lha daoraaaa In tha vlaual light ourva baglnnlng in 

uhan ooMblnad with tha axpansion valoolty in tha shall of 

10^B/a roughly oorrasponds to tha tiaa at whioh tha initially ajaotad 

anvalope aatarial would hava raaohad tha innar radius of tha obaarvad 

dust shall and gas kinatio taaparaturas would hava baan batwaan 1000K 

and 200(NC. Unfortunataly, partiola dansitias at this radius would 
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hava daoraasad to approxiaataly 10 atona/oo (Prassura ■ 10 atai). 
Nuolaation at such low dansitias is axtraaaiy difficult without tha 
invocation of nassiva dansity anhancaaants (of ordar 10^) dua to 
turbulence. Magnetic conflnaaant, etc. 

lha observed infrared spactrui of tha dust can be Modeled by an 
aaorphous silicate Material, heated oy stellar radiation, which 
re-emits in the Infrared. There is little doubt that this dust must 
have condensed in the stellar outflow. The size of the grains is 
postulated to be ^1 Micron. If this is true, then significant grain 
growth Must have occurred after the initial nuoleation process, 
because tha region of nucleatlon and subsequent grain growth must of 
necessity occur at a radius ouch greater than that of the stellar 
photosphere, a wore careful study of the region interior to the dust 
shell Might be possible. This is in contrast to the situation in cool 
H stars where grain fomation occurs very close to the photosphere. 
Such a study Might detect molecular precursers to the observed solids 
and could possibly yield valuable information on the kinetic and 
vibrational temperatures of these molecules in the region of grain 
foraation. 
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Forrest tt si. (1978) hsve studied the infirsred speotrui of the 
varlsble star OH 26.5 * 0.6 over a period of two years in the 
wavelength range 2 alorons < l < 40 aiorons. The shell sround this 
star is thought to produce aore than 70 magnitudes of visual 
extinction. The underlying source is thought to be an H type variable 
star. The spectrisi of the shell shows prominent id and 20 micron 
absorption features and is therefore thought to consist of silicates. 
The mass loss rate for this star is estimated to be ^10 solar masses 
per year. It is probably an extreme exaaple of mass loss in late 
stars and most likely belongs to an evolutionary stage which occurs 
after that studied by Hagen (1978). The infrared spectrum of this 
source is presented in Chapter 6 along with a discussion of tne 
interpretation of the residuals in the broader 10 and 20 micron peaks 
in terms of laboratory produced amorphous silica, Si20^ and magnesiui 
silicates. 

Hackwell, Genrz and Grasdalen (1978) observed a sudden infrared 
brightening in the Wolf»Rayet star HD 193793 between November 1976 and 
June 1977 at X s 2.3 microns and x s 10 microns. This was followed by 
a decline in intensity at x * 2.3. 3.6 and 4.9 microns and a slow 
brightening at x > 8 microns. This was interpreted as the 
condensation of an optically thin dust shell from material ejected by 
the star in 1976 at a velocity of approximately 2.7 x 10^m/s. The 
infrared speotrua of this object does not show a 10 micron feature and 
the dust is therefore thought to be either iron or graphite. Detailed 
characteristics of the grain speotrua are insufficient to accurately 
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■od«i oondMMtlon in this 'systw although an attaupt has baan aada to 
obtain a aaasura of tha aass loss froa tha atar. This is astlaatad to 
ba 10 solar aassas if tha grains ara iron and v>10 solar aassas if 
thay ara graph! ta. Tha tMparatura of tha photosi^ra of this star is 
- 30.000K. 

Condansation in this systaa is thought to ooour in uavas - unlike 
tha continuous outflows saan in tha souroas diaoussad above. The 
authors believe that tha infrared brightening of this star, observed 
in 1970, was due to an ejection of material which occurred in 1969 
while that observed in 1977 was due to an outburst in 1976. If this 
star is periodic, than another ejection is due in 1982-83 and should 
be followed by an infrared brightening in 1983-84. A search for 
infrared, optical and ultraviolet amission (or absorption) by the 
molecular preoursers to the expected grains - possibly SiC, FeC, C^, 
C^, etc. could prove to be quite informative. 

Koorneef and Sitko (1979) have shown that the wavelength 

dependence of the extinction toward HD 193793 could be explained by a 

population of relatively small graphite grains which originated in the 

stellar wind and were slowed by the interstellar mediia at a distance 
17 

of <^10 cm from the star. This further supports the contention of 
Hackwell, Gehrz and Grasdalen (1978) that grain formation can occur 
periodically in this star. 

In summary, numerous examples of grain formation exist in a 
variety of stellar types. These range from stars with photospheric 
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tnperatures slgnlfioantly in axoaaa of to thoaa with 
taaparaturea of only a faw tiaaa 10^* Son# atara oondanaa ailioata 
duat tmila otnara nuolaata iron or graphita graina. Ua obaarva duat 
foraation in atara tdiioh ejaot aa littla aa 10*^ aolar aaaaaa par yaar 
aa wall aa in aouroaa which night axpall aa much aa 10 aolar naaaaa 
par yaar. Maaa loaa oocura continuoualy in aona atara and aa diacrata 
epiaodea in othera. The duat anvelopa can ba optically thin or can 
produce enough extinction to complataly aaak the underlying atar. No 
one atar can be aaid to typify thoae which produce graina. Each nuat 
be treated aeperately. Yet each environaent in which the nucleation 
proceaa can be denonatrated to occur furthera oir knowledge of the 
proceaa in general, aince each aeperate caae nuat be treated in any 
conprehenalve theory of grain fomation. 

Condenaation in Novae : Obaervationa of duat fomation in novae were 

reported by Hyland and Neugebauer (1970) and Geiael, Kleinman and Low 
(1970) for Nova Serpentia 1970. Viaual obaervationa ahowed a gradual 
decline in the luninoaity of the nova until <t day 60 when a draatic 
reduction in intenaity began. Sinultaneoualy, infrared obaervera 
noted a sharp increase in infrared enission. The nost straightforward 
interpretation of these observations is that, between days 60 and 70, 
an optically thick dust shell fomed froo the aaterial in the 
outflowing nova ejecta. This shell absorbed the stellar visual and 
ultraviolet radiation and reradiated the energy in the infrared. 

This pattern has been repeated for several more recent novae, 
e.g. Nova Vulpeculae 1976 (Ney and Hatfield , 1978) , Nova Serpentia 1978 
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(Sskody, «t al*, 1979; Oahrt* at al.« 19d0a) and Nova Qfinl 1978 
(Gahra at al., 1980b). It should ba notad that not all nova font 
duat. An axaapla of thla latar type la probably Nova Cygni 1975 
(Gallaghar and Nay, 1976) which was followad for v>1 yaar aftar visual 
aaxlauB and lAiloh showad no avldanoa for dust oondansation . 

Savaral aodala hava baan proposed to axplaln (Clayton, 1979; 
Clayton and Hoyla, 1976; Clayton and Ulokranaslngha, 1976) or to 
pradlct (Galli^har, 1977) tha conditions undar which dust will fora In 
novaa. In tnls latter work, Gallagher (1977) contends that If 
lonlxstlon of tha ejected aatarlal beccaas significant before 
conditions suitable for dust oondansation hava developed, then dust 
will not fora. He derives an aaplrlcal relationship for the tlae (In 
days) aftar aaxlaua light at which grain foraatlon will begin 
(Equation 5) as a function of the lualnoslty (L) of tha nova and 

. (320/v)a/L,„)''^ (5) 

the expansion velocity (v) of the ejecta in ka/s. By setting t^ equal 

to t^ (the tlae at which significant Ionization developes In the 

envelope) ha predicts that novae with lualnosltles greater than about 

9 X lO^L . will not produce dust. Nova Cygnl 1975 Is cited as an 
sun 

exaapla of this class of outburst. 

Althoi^h theraal radiation frrz dust grains surrounding novae is 
alaost certain to be responsible for the observed Infrared flux. It Is 
not universally accepted that these grains fora In the nova outburst 
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th»n ynd«r study. It ftss bssn suggsstsd that prs-aatlstlnf fralnst 
sirsady surrounding tha nova bafora tha outburst, oould ba bastad by 
tha Intansa ultravlolat radiation flald uhloh Is knoun to Inoraasa 
with tins sftar manimm llgbt. In psrtloulsr. Bods and Evans (1960) 
nava propoaad aodais for Novsa ;^rpantls 1970 and Aqullaa 1975 In 
whloh tha daoraaaa In vlalbla radiation la a alapla oonsaquanoa of tha 
shift of tna bulk of tha radlativa anargy Into tha ultravlolat as tha 
outer layers of tha nova baoona optically thin and tha hotter surface 
of tha underlying Mhlta dwarf baocxia exposed. Because they postulate 
that tha pra-axlstlng grains have a larger absorption cross section 
for these ultraviolet photons, the emargenca of such Intense 
ultraviolet oaisslon rapidly heats the surrounding aatarlal and thus 
causes tha observed infrared ealssion. 

Bode and Evans suggest that a distribution of 0.01 micron 
graphite particles - left over from previous outbursts - could account 
for the observations. This seems quite plausible both since graphite 
has an increased absorption cross section in the ultraviolet near BOnm 
and 220ru and since such aaall grains would tend to remain within a 
nova system for relatively long periods of time (Bode and Evans, 

I960). A similar interpretation of Nova Cygnus 1976 has recently been 
published (Strickland et al., 1961) in which the pre-existing dust is 
thought to have been heated by the absorption of photons in 
ultraviolet resonance lines rather than from the emerging ultraviolet 
continuuB of the white dwarf. The observed infrared flux requires 
that only about 30 1 of the resonant radiation be absorbed. The 


oaioulcMons of Striokiand at al. (1981) indloata that auoh a high 
affiolanoy la not unraaaonabla , 

Although Strickland at al. (1931) aatlnata that this abaorptlon 
has a ralatlvaly aaall affaot on tha ratios of tha various rasonanoa 
Unas In tha aaargant spaotrua of tha nova, thalr nodal would scan 
nora convincing If thay could show an axanpla of a llna ratio for 
which this was not tha casa. Slnca graphlta grains show enhanoad 
absorption at dOna and 220nn, a rasonanoa llna naar ona of thasa 
wavalangths should show a considerabla daoraasa In Intanslty whan 
comparad to ona batwaan thasa paaks. Although it is inpossibla to 
measure line intensities for astrononical objects at 30m, it may be 
possible to observe this effect near 220m. 


b. Proposed Regions of Grain Formation 

It has been observationally established that grains form in the 
atmospheres of cool stars - and even in the outflow from other stellar 
sources at considerably higher temperatures. Additional regions of 
grain formation - for which no observational evidence exists - have 
also been proposed. Supernovae have been discussed as an important 
source of material enriched in certain short-lived radioisotopes. 
T-Tauri stars and/or primitive solar nebulae could be responsible for 
the establisment of the observed grain size spectrum. Pl{>netary 
nebulae might be important sources of interstellar grains. Each of 
these regions will be briefly discussed in the following paragraphs. 
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Infrared obatrvatlona of planetary nebulae have revealed the 
presence of features usually attributed to dust fomed in both carbon 
and oxygen rich systeas (Aitken et al., 1979)* Therefore aoae 
fraction of interstellar silicate and graphite grains probably fonied 
In planetaries of the appropriate composition at some evolutionary 
stage, unfortunately, the chemical composition of the material 
ejected from a proto planetary nebula is the subject of some debate 
(Zuckeman, 1980). The dynamics of the ejecta and the conditions 
under which the observed grains may have formed are even less certain. 
Natta and Panagia (1981) have observed that the grains appear to be 
concentrated in tne central regions of the nebulae. This would 
indicate either that older grains (those farther from the star) have 
been partially destroyed or that the concentration of grain forming 
material was more abundant in the inner regions of the envelope when 
it was ejected by the parent star. There is other observational 
evidence however, %<hioh indicates that grains may be destroyed in the 
interior regions of planetaries (Harrington, 1981). 

It should be noted that the observed profile of the silicate 
feature in NGC 7027 is not that typified by the signature of the 
Trapezium - the archtypical 'silicate*. Rather it appears to have 
sharper features > such as those observed in mineral spectra - than is 
usually seen in astronomical sources (Jones et al., 1980). This could 
be the result of an annealing process on materials originally of a 
more amorphous natue, possibly even similar to the grains freshly 
nucleated in the atmospheres of cool stars. Of course, the 
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possibility slMiys mists that tbs prossirts in tn« original ajaotad 
anvalopa wars hign anough that a<»a aora oryatallina apaoiaa uy hava 
fonsad. Ona dafinita oonolusion is that tha aaohaniaa by idiioh tha 
grains obsarvad in planatarias nuolaatad and avolvad will raaain 
unknown until auoh aora inforaation bacoaas availabla on tha avolution 
of planatary nabulae thaasalvas. 

T-Tauri stars ara known to axhibit rathar axtansiva dust shalls. 
Thasa probably consist of dust originally contained within tha cloud 
before it had collapsed (Finn and Siaon, 1977; Harvey, Thronaon and 
Gatley, 1979). lhasa stars ara therefore not an iaportant source of 
the total aass of grains in the interstellar aadiua. They aay however 
serve as an extreaely iaportant contributor to tha observed grain size 
spec true. Specifically, Burke and Silk (1976) propose that processes 
such as grain aelting, shattering and nucleation can occur in the 
proto stellar envelope of a newly foraed star. Tnese aaaller grains 
could subsequently be reinjected into the interstellar aediua via 
radiation pressure or in the T>Tauri 'wind'. This aechaniaa aight 
provide a larger nunber of snail grains frca the saae tot«-l aass. 

No convincing observational evidence has ever been presented 
which indicates that grains fora in the ejecta frca supernovae. In 
fact, supernovae shocks are usually considered to be a priae 
destruction aechaniaa of even highly refractory grain cores (Salpeter, 
1977). Anoaalies discovered in the refractory chondrules of 
carbonaceous chondrites do however indicate that scae of the aaterial 
kdiich foraed the solar nebula was enriched in supernova produced 
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iaotopti (CMtron and Druran, 1977). ThU diaoovary liaa proiptad aoaa 

\ 

raaaarohara to oaloulatt tdaraodynaalo oondanaatlon aaquanota for 
auparnovaa ajaota (Lattiaar, SohraM and Qroaaian. 1978; Lattlaar and 
Groaaian, 1978). Suoh a ayataa la, of oouraa, in a atata of aatrana 
diaaquiiibriua (Noth and Donn, 1981; Cnaptar 3) and raaulta which raiy 
on tha attalmant of such aquilihria ara tharafora axtraaaly dubioua. 

A nora likaly aoanario for tha obaarvaa laotoplo anriohaant Mould ba 
tha gradual daplatlon of tna gaa phaaa rafraotorlas onto intaratallar 
grain ooraa whioh had survivad tha auparnova blast uava. Thasa 
ooapoaita graina could hava subsaquantly undargona raprooassing in tha 
prinitiva solar nabula. 


c. Ragions of Hantla Formation and Grain Growth 

As discussad in saotion III, there exist ragions of higher tnan 
normal density which exhibit anomolous ultraviolet extinction. Mathis 
snd Uallenhorst (1981) have shown that a normal distribution of 
typical interstellar dust grains which becomes deficient in small 
particles can reproduce tha observed extinction of some regions, kdille 
others rec^ire that the nuiber density of larger grains increases. 
Although it is possible to preferentially eliminate smaller grains by 
evapor,' ion (McCall, 19r9: bragar, Gahrz and Hackwell, 1981 ), such a 
process can not explain regions of unusually large grains (Mathis and 
Uallenhorst, 1981 ). 
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ThrM proof SMS oould poMlbly l>t rfopontiblt for fuoh 
obtf rvatiOM ; aiaplo grain growth, aantla foniation and grain 
coagulation, it ia poaaibla that all three prooeaaea contribute to 
the fomation of larger graina to acae extent. The proceaa which 
doainatea in any particular region would, of courae, depend upon the 
conditiona in the region. 

Whittet and Bladea (19d0) have concluded that grains In the p 
Ophiuohua cloud do grow by the accretion of volatile mantles since 
the Infrared extinction spectrua of the material in this cloud shows 
no sign of featkw'es at 3.1 or 3.4 microns (water ice and C-H stretch, 
respectively). This is in contrast to the conclusion reached by Cohen 
(1977) who found that abundance constraints, coupled with the increase 
in the average size of the particles necessary to explain the 
polarization observations in such dense clouds, required the 
condensation of CNO rich mantles. She did show however, that the gas 
phase abundances of even fairly volatile metals (i.e. Na and Ca) in 
such clouds was much less than in the average interstellar mediui. 

She concluded that an increase in the mass of absorbing material was 
necessary to explain her observations and therefore ruled out grain 
coagulation as the mechanism of grain growth. 

Unfortunately, Jura (1980) has suggested that a model of the 
extinction towards p Ophiuchus can be based upon coagulation since the 
visual opacity per gram is unusually low towards this star (Bohlin, 
Sav%e and Drake, 1978). This is felt to result from an increase in 
the average size of the graphite component of the 'average' grain 
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population Oy about • factor of thraa. In Jura* a (19S0) aodal 
ailleata graina do not ooagulata, 

Tnia bahavior la In oontraat to tha bahavlor auggaatad by tna 
anaiyala of Nathla and Wallannorat (1991)# l*a. that mail graphlta 
grains only bagln to 'dlsappaar* aftar both tha uppar and louar slca 
llflilts of tha sllloata distribution hava Inoraaaad. Slnoa It has baan 
shown obsarvatlonally that graphlta grains do not contribute to tha 
polarization of such sources, it is difficult to inaglne that tna 
growth of these grains will Increase the value of Xaiax. Jura (1979) 
does not coonent on either of these problens. 

Jura (1979) does however mention tha fact that tha tloasoala for 
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coagulation may be on the order of v>10 s. This tlmasoala la 
comparable to the lifetimes of the clouds thamaalves. If tha grains 
in such clouds are charged (Watson, 1972; Jura, 1976; Oraine, 1979) 
then the timescale for coagulation could increase substantially. Such 
an increase would ir.tan that coagulaMon is unimportant for all but the 
most ancient or compact clouds. 

In summary, grain growth resulting from the adsorption of gas 
phase metals and metal oxides as wall as that ^'esulting from the 
formation of volatile mantles are likely to be important mechanisms by 
which the averse grain size in dense clouds increases. It is 
important to note that processing of volatile mantles by stellar 
ultraviolet radiation might resia. in a significant reduction in tne 
strength of the infrared signature of such materials at j. 1 and J.4 



microns (Sagan and Khare, 1979)* This might indicmta that 'old' 
mantles would be very difficult to observe even when present as a 
quite substantial percentage of the dust, and could thus explain the 
failure of Uhittet and Blades (1980) to detect the expected signature 
of a 'volatile' component. In contrast, grain coagulation is not 
likely to be an important mechanias for the production of larger 
grains in any but the densest and most ancient clouds. 


d. Grain Destruction Hechaniaas 

Several mechanisms have been proposed which could contribute to 
the destruction of interstellar grains. In addition to the possible 
reprocessitig of grains during star formation (Section 111), 
sputtering, grain-grain collisions, photodesorption, evaporation and 
grain surface reactions/ chemical sputtering may all affect the grain 
composition and size distribution. 

Sputtering in both high and low velocity shocks is expected to be 
the most effective mechanisn for the destruction of refractory grains. 
These shocks occur primarily as a result of supernova explosions and 
intercloud collisionr. Grain sputtering can also occir in H 11 
regions, planetary nebulae and the intercloud medium, Cowie (1978) 
has Shown that sputtering can be extremely effective even for 
relatively low velocity shocks (40-50km/s). In his model, destruction 
occurs when a shock travels perpendicular to the magnetic field of the 
cloud. Charged grains are accelerated about these compressed field 
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lin«s and aputtarad by Ha and CNO atoaa. In tilgbar valooity ahooks 
(80-100ta/s) a aignlf leant fraction of allioata or graphite gralna 
oould ba daatroyad. Shull (1978)» In a slnilar study, has raaohad 
essentially the sane oonoluslons. 

Barlow (1978a) has developed eaplrloal fomulae to predict both 
the sputtering threshold energy and the sputtering yield based upon 
extensive review of the experimental data. He concludes that 
supernova produced shockwaves are likely to be the most laiportant 
single mechanism for the destruction of grains. Sputtering by cosmic 
ray particles is not found to be effective. 

Both Shull (1977. 1978a) and Havnes (1980) have studied the 
effect of Intercloud collisions upon the average grain population. 
Shull (1977) considers the results of hydrom^^netic shocks produced in 
low velocity collisions (20>50toB/s) . He finds that between ;j-l0S of 
the dust can be destroyed even in shocks in which the fractional 
ionization of the gas remains low. Destruction results from the 
differential acceleration of the grains within the cloud and 
subsequent grain-grain collisions. Havnes (1980) considers the onset 
of a two stream instability in the shocked region whereby grains less 
than <riD0nn in diameter are brought to rest with resf>ect to the gas on 
a short timescale while those of larger diameter continue at 
approximately the original velocity of the cloud. The net result is 
the formation of two grain populations which stream relative to one 
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anothtr Mi thin tha sm« oloud. Grain-grain oolliaions tharafora oan 
baooaa auffioiantly fraquant and anargatio to oausa aignifioant grain 
daatr Motion. 

Othar procesaes, of little iaportanoa to the average graina in 
the interataliar aediua, can noMaver Oaooaia important Hi thin certain 
astrophysioal anvironaanta. I^aine and Salpeter (1979) have studied 
the fata of dust ambadded uithin a hot plaama (10^-10^K). Although 
sputtering ia important at higher densities, at high temperatures 
(T>10 K) ion field emission can became a significant process by uhioh 
grain material is returned to the gas phase. In dense clouds, 
chemical reactions at grain surfaces - vAiich lead to the production of 
volatile oxides or hydrides - can erode and eventually destroy small 
grains (Barlou, 1978b). Oust in compact (N^>10 atoms/cc) H 11 regions 
can be destroyed by reactive sputtering with H, N and 0 atoms (Draine, 
1979) although the process is not very effective at lower densities. 
Photodesorption is an important destructive process only for weakly 
bound van der Waals mantle materials (Barlow, 1978b) although it does 
help to eject the products produced by grain surface reactions. 
Photodesorption is, however, an effective mechanism for the prevention 
of mantle formation. It may also help to explain the depletion 
patterns of the non-metallic elements in diffuse clouds (Barlow, 

1978C) . 
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V« Su— ary 


Rafractory Interstallar grains are produced as a result of aass 
loss in a wide variety of stellar sources. These grains are then 
modified by processes in the Interstellar medium idilch Include 
sputtering. Ion emission, grain-grain collisions and mantle growth. 

The average resultant grain population can be characterized as a 
fflixtts*e of several components which must Include graphite (or 
Carbyne?) but can also Include "silicates”, silicon carbide. Iron and 
magnetite. The size distribution of the average population Is likely 
to be of the form n(a) « a'** for q « 3.5 and a between 5 and lOOOnm. 

It Is possible to distinguish at least two grain populations from 
measurements of the wavelength dependence of polarization. One is a 
population of small graphite grains which may serve as the carrier of 
the diffuse bands (at least those at 443nm, 57dnm and 628. 3nm.) and 
which are not aligned by the galactic magnetic field. Ihe other Is a 
population of 'silicate' grains which may be associated with small 
iron or magnetite grains and which is responsible for the polarization 
of starlight. A third population of grains may consist of extremely 
small silicates which are not aligned and which are responsible for 
the isotropic scattering observed in the far ultraviolet. A fourth 
population of grains is found in dense clouds and is characterized by 
a shift towards larger average grain size. This shift may be caused 
by the condensation of either refractory or volatile mantles but might 
also be the result of coagulation. 
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Dust, frssniy nuolsatsd in tha outflows of oartain stallar 
souroas, aay ba oharaotarliad as a fifth population sinoa suoh 
notarial raprasanta solids as yat umodifiad by (Mrooessas in Uia 
intarstallar nadiua. lha infrarad spaotrua of ona auoh ragion is 
discussad in Chaptar 6 of this dissartation in relation to the 
laboratory spectra of annealed Mg-SiO aaokas. It is dosonstratad that 
siailar annealing processes nay occur after the initial ejection and 
nucleation of solids in the envelopes of such regions and that these 
nay be observable with available infirared spectroneters. 

It is probably possible to define distinguishing characteristics 
to categorize grain populations found in specific astrononical regions 
such as T-Tauri stars, planetary nebulae or H II regions. Careful 
laboratory studies to detemine the effects of the various processes 
which are thought to be inportant in the interstellar nediua and in 
the so(e*ces thenselves, upon the naterials expected to account for the 
najority of the solids in such regions, night result in a better 
understanding of the relative importance of these effects and of the 
most probable history of the material within the region. A review of 
such experimental efforts will be presented in Qiapter 2. Hy own 
efforts along these lines will be presented in Qvapters A, 5 and 6 of 
this lissertation. 
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ChaptT 2 


ExpTiJitntal and Thporptioal Studip* of Inttrattllar Grain* 


Laboratory studlaa which relate to interstellar grains can be 
roughly divided into two classes. In the first class of experiments 
the primary aim of the research is the production of laboratory 
analogs of suspected interstellar materials. This is usually done in 
order to reproduce certain specific observations - such as tne shape 
or position of an absorption or emission feature - and thus to 
'identify* the component responsible. The second class of experiments 
are those which attempt to define the conditions under which grains 
might form. Currently, most models of grain formation are based on 
the classical theory of homogeneous nucleation. 

This chapter is divided into four main sections. In the first, I 
review experimental studies which were primarily intended to 
'identify* interstellar materials. In the second, I review classical 
nucleation theory and briefly discuss its limitations and possible 
alternative formulations. In the third. 1 discuss the application of 
this theory to the problem of the formation of interstellar grains. 

In the fourth section I review the experimental data applicable to 
such systems and the implications of this data for the models of grain 
formation currently in the literature. 
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I* Laboratory to Inf rittllar Dust 


Until ralatlvaly raoantly tha only "axparlaantal" avldanoa 
ralatad to tha ohaaioal ooapoaition of intaratallar gralna waa 
obtalnad from tha atudy of aataoritaa. Spaoial aaphaaia uaa plaoad 
upon information glaanaO froa tha atudy of oarbonaoaoua ohondritaa, 
ainca thaaa ara thought to ba tha aoat 'priaitiva' olaaa of aataoritaa 
and tharafora tha oioat rapraaantativa of intaratallar aolida (Caaaron, 
1973). 


Unfortuiataly, the temperature and pressure at which aeteoritic 
materials were originally compacted is not well determined. Ihe 
region of the solar nebula in which these processes may have occurred 
as well as tha (devious history of tha matarial before coagulation 
into meteor parent bodies is even more speculative. The processes by 
which these parent bodies are ruptured and the effects of such 
catastrophic events on the structure of the material can only be 
guessed. 1b understate the problem, aeteoritic specimens were not 
produced under well controlled laboratory conditions. 

Within the past 10 years, the study of possible analogs to 
interstellar materials has became more fashionable. Such studies will 
be divided for convenience into four general classes. The first class 
of studies attempts to use common materials as a basis for the 
identification of observed spectral features. Such work often reports 
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the aeasureoient of the epeetrel featiares of a Mineral > aoMetiaes as a 
function of teaperature, grain size or radiation damage. The second 
class of experiaent reports the effects of heat« radiation, or 
electrical discharge upon the spectrua of organic • or at least carbon 
rich - materials. The third type of study aeasures the optical 
properties of glasses or materials produced by the sputtering or rapid 
evaporation/ condensation of refractory materials. The fotarth category 
of experiaent includes those studies which produce smokes. It should 
be noted that these categories are completely arbitrary and that there 
may exist considerable overlap between experiments in separate 
classes. 

a. Mineral Analogs to Interstellar Dust 

The thermodynanic calculations of Gilman (1969) predicted the 
formation of olivine grains in the envelopes of cool stars. An 
absorption feature at 10 microns was sik>sequently detected at low 
resclution in such sources and olivine became the mineral of choice 
with which to model circimstellar dust. High resolution spectra for 8 
< X < 13 microns however do not confirm the presence of olivine in 
ciicuostellar envelopes. The observed feature lacks the sharp fine 
structure characteristic of mineral spectra. The search therefore 
began for a material which had a relatively featureless 10 micron band 
and Uiiich might also be a natural constituent of astrophysical 
env irornents . 
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Day (1974, 1976a,b) dascribtd a laboratory synthaslzad Morphoua 
aagnaalm silicate which has an infrared absorption speotruM 
reseablins ^at of oircuestellar 'silicates'. The aaterial was 
prepared as a gel fron an aqueous solution and its infrared spactrua 
recorded as a function of annealing teaperature. Annealing in a 
vaouuB transforas the spec true of this aaterial froa one idiich 
reseables that of interstellar dust to one identical to forsterite 
(the aagnesiuB end meaber of the olivine solid solution series). It 
is interesting to note that when Day (1974) heated a sample of the 
matrix of the Murchinson Heteorite (a type II carbonaceous chondrite) , 
the infrared spectrua of which also reseables that of the interstellar 
feature, he obtained a spectrun which exhibited many of the prominent 
olivine signatures. 

Zaiko^ki, Knacke and porco (1975) have suggested a 
phyllosilicate model of interstellar grains based primarily on studies 
of type I carbonaceous chondrites. Their analyses C( ifirm that 
phyllosilicates are the dominant silicate minerals of such meteorites 
and that the infrared spectra of several phyllosilicates (i.e. 
montmorillonite or chlorite) resemble, to first approximation, the 
spectra of astronomical sources. 

Unfortunately, associated peaks near 3 and 6 microns, produced by 
hydroxyl groups within such minerals, are not correlated with the band 
strength of the 10 micron feature as would be predicted by such a 
model. In addition, the direct condensation of such complex species 
from the vapor . needed to explain clrcuostellar ob*’ ;rvations - 
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should prove to b« extrsasly difficult. However it is easy to iassine 
the forastion of such aaterisls frea pre-existing axorphous silicates 
in regions of higher pressure such as the priaitive solar nebula. 

Knacks and Kritchaer (198O), in a study of the infrared spectra 
of carbonaceous chondrites, found that hydrated phyllosilicates aight 
account for a long wavelength wing found in the 3. 07 aicron water ice 
absorption featire. The hydration of axorphous silicates, on a axall 
scale. Bight be possible given the degree of reprocessing expected for 
the average interstellar grain. However, alternate explanations for 
the observed shape of this feature exist and will be discussed in the 
next subsection. 

Because high teaperature condensates such as olivine and 
enstatite are predicted to form in circuastellar environments, it has 
been suggested that the irradiation of such materials by energetic 
particles could alter their spectra sufficiently to match that 
observed in the interstellar mediin (Uickramasinghe , 1971). Pybwad 
(1971) investigated the effects the irradiation of silicate minerals 
by KeV energy protons while Day (1977) used 1-2 MeV protons. Neither 
found any significant alteration of the spectra even for radiation 
doses which were factors of 10-100 higher than those expected for 
typical grains. 

KrStchner and Huffman (1979) have measured the optical constants 
cf an extremely disordered olivine sampl' produced by the exposure of 
a polished, natural olivine to a dose of 5 x 10^^ (1.5 MeV Neon ions) 
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oa TtM dlalMtrio fuDOtions fbr tti« 8 aiorons < x 

< 30 alorons have b««n datarainad froa an aaaXyais of tha apaoular 
raflaotanoa of tha aaapla. Tha aatlnotlon produoad hy kylaigh 
partlolaa of this ooapositlon wara ealoulattd uslnt Mia lhaory. 
Extinotion aaasuraaants for an aiorphous oiivlna aaoka wara also takan 
and found to agraa with oalouXations. Bacauaa tha radiation dosaga 
usad to produea the aaorphous olivine was auch higt^ than that 
axpeotad (Mrar tha lifatiaa of typical interstellar grains, nothing 
oould be concluded about the foraation of such aaorphous aatarials, 
Krttchaar and Huflkan (1979) did conclude however that the aatarial 
responsible for the 10 and 20 aicron features oould be an aaorphous 
'v)iivine* . 

X have already noted that tha best aodel for the visible and 
ultraviolet extinction and polarization aeasureaents (Mathis and 
coworkers- see chapter 1) is based upon the optical properties of wail 
known terrestrial aatarials. Such infomation has been available in 
the literature for nany years for most of the less complex and more 
prosinent mineral species. 


b. (Organic Pcljners and Ices 

As discussed in section III of Qiapter 1, Hoyle and coworkers 
have suggested that the observed infrared spectra of astronomical 
sources could be explained in terms of models based upon the 
laboratory spectra of polysaccharides (Hoyle and Wickramasingha, 1977) 
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or bioohaaioal oaapoimds (Hoyle end Wiokraia 8 inghe« 1981). Iheir 
visiOle^ltraviolet extinction eodel is based upon laboratory 
aeasuraaents of the optical properties of viruses, bacteria and 
biocheaical chroaophores (Hoyle and Wickraaasinghe , 1979). Much of 
the current controversy surrounding this aodel revolves about the 
recent detection of a 3 .S aicron feature in the spectrua of IRS-7 
(Wickraaasinghe and Allen, 1980). liMile Ouley and Uilliaas (1979, 
1981) interpret the feature as the signature of C-H functional groups 
foraed on the surface of interstellar graphite grains, Hoyle (1981) 
prefers to interpret this feature as indicative of interstellar 
microbes. 

Little more need be said about such models except to point out 
that they are obvious examples of what happens when the 
'identification* of a grain material rests solely upon the shape of 
its spectrin. Although this is indeed an important constraint, one 
must also consider such factors as the relative stability of the 
material in the regions in which it has been observed and the 
mechanism by which such materials might form. 

Studies of two other types of carbon based polymers have been 
performed in the laboratory. Khare and Sagen (1971, 1973, 1979 and 
references therein) have produced a relatively stable organic polymer 
by electric dischargii through, and ultraviolet irradiation of, 
mixtures of fairly siiople starting materials (i.e. CH^, NH^, H^O, H^S, 
etc.) The properties of this material vary from run to run but may 
nevertheless represent a small fraction of interstellar material. 


This coaponsnt could bo produosd as mantled grains approach the adgas 
of interstellar clouds and thus become exposed to the interstellar 
radiation field. 

Carbynes are solid allotropes of carbon formed from alternating 
triple and single C-C bonds, bhittaker and Wolten (1972) have 
synthesized carbynes at temperatures in excess of 2600K in the 
laboratory and Whittaker (1978t 1979) subsequently found several types 
of carbynes in samples of six terrestrial graphites. Carbynes have 
also been identified in meteorites (Hayatsu et ai., 1980) where they 
have been shown tc be carriers of primordial noble gas (ttiittaker et 
al.. 1980). Hayatsu et al. (1980) have demonstrated the formation of 
carbynes from the disproportionation of CO in the presence of a 
enromite catalyst at temperatures of only 3CK)-400C. Webster (1980) 
has suggested that carbynes may be the primary condensate in the 
atmospheres of carbon stars. For this reason, it is important to 
obtain the ultraviolet extinction spectra of these materials. Such 
experiments are currently in the planning stage (Hecht, 1981) • 

A study of the condensation of water at the low temperatures 
expected in dense interstellar clouds indicates that an amorphous ice 
is formed (Leger et al.,1979). Ihe spectrin of this amorphous solid 
more closely resembles the observed shape of the 3. 1 micron feature 
than does a more crystalline species. This study also indicates that 
the band strength of the 3,1 micron feature in amorphous ice is only 
v2/3 of that in crystalline ice. 
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Gr««nb*rg and ooMorkara (Hagan, Allaaandola and Graanbarg, 1979; 
AUaaandola, Graanbarg and Noraan, 1979) Hava raoantly daaoribad 
axperimantal atudias of tha low taaparatura (<i>12-42K) ptotolyaia of 
alaple CO, CH||, H^O and NH^ aixturaa. The purpoaa of tha expariaanta 
was to Inveatlgata the produota fonaad by auoh prooesaaa in tha 
■antlas of Intaratellar grains by tha action of ultraviolet radiation. 
The observed products were analysed as a function of the length of 
photolysis by infrared, visible and ultraviolet spectroscopy. In 
addition, the nass spectra of selected residues were obtained after 
tt^ system had been brought up to room temperature. It was found that 
some of these residues were stable to temperatures as high as 
550-600K. 

Hagan et al . also found that a large nunber of molecules and 
radicals observed in tha interstellar mediua can be produced by such 
processes. It was also shown that grain mantles which contain such 
species also contain a large quantity of stored chemical energy. This 
energy could be released in relatively low velocity grain-grain 
collisions and result in the evaporation of the mantle and the release 
of a significant nuaber of complex organic molecules. Such processes 
might be the primary mechanism by which the more complex organic 
species, which are observed in dense clouds, are produced. 

Finally, it was shown that ultraviolet irradiation of CO/H^O 
mixtures results in a feature at 3. 1 microns which closely resembles 
the shape of the observed band attributed to water ice in the spectrum 
of the KL Nebula. It is not known (by this author} if this results 
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from a braatalown of bht H2O lot struoturo or froa a ohaaioal raaotion 
with CO* Studiaa of tha affaot of tha irradiation of othar aixtoraa 
upon tha ahapa of thia faatura ara raportad to ba in prograaa (Hi^an, 
Aliaaandola and Graanbarg, 1^79). 


o. Sputtering and Laaar/FIaah Evaporated Condensataa. 

Lord (196b) vaporized poudara of iron oxide, aagneaiuB oxide, 
aluainui oxide and dunite with a flaah diacharga lamp but waa unable 
to either cheaically or structurally Identify hia condenaation 
products which foraed es airrors as well as whiskers and spherules on 
the glass walls of his system. Because he was unable to control tha 
ccaposition of the vapor, could not achieve a homogeneous vapo* base 
and could not control the rate of condensation, he abandoned the 
project. 

Blander et al. (1969) used a laser to vaporize basalt targets and 
found (not surprisingly) that the recondensed material was enriched in 
volatiles while the residual basalt became enriched in refractory 
oxides. Heyer (I97I) controllably condensed a vapor produced by the 
sputtering of silicate targets in an argon glow discharge onto 
substrates at known temperatures. Silicon, magnesium and iron were 
found to be depleted relative to calcium and aluainun at substrate 
temperatures above 50OC while the Ca/Al ratio was found to remain 
unchanged over a wide temperature span (100C < T < 850C) . He 
therefore proposed that although the bulk chemistry of chondritic 
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Mttoritts Might b« rtlatttd to a high taaparatura oondansation prootaa 
in a ralativaly hoaoganaous aolar nebula, individual chondrites are 
unlitceiy to be priaordial gas phase condensates because of their 
variable Ca/Al ratios. 

Stephens and Kothari (1978) have used a pulsed laser to evaporate 
sanples of astrophysical interest such as olivine, pyroxene, aluaina, 
carbon, etc. Ihe samples recondensed in the laser pluae, usually at 
pressures in excess of one atmosphere in an Ar, 0^ or buffer gas. 
Silicates vaporized in as well as caloiin carbonate and carbon 
vaporized in Ar or recondensed into glassy spheres 20»30nm in 
diameter. Aluaina and silicon carbide recondensed into crystalline 
grains. It is interesting to note that the condensates formed in 
their system always formed chainlike necklace structures. 

Stephens and Russell (1979) used the system described above to 
produce samples of glassy olivine and enstatite from natural samples. 
The infrared spectra (4 microns < x < 14 microns) of these amorphous 
glasses were then compared to samples of the same composition which 
had been ground to similar dimensions. These spectra were also 
compared to the observed spectrum of the Trapezium. It is shown that 
although the spectra of the amorphous glasses are narrower than th«f 
observed profile of the Trapezium, variations in the composition of 
the grains or gradients in the source temperature of the dust cloud 
could broaden the laboratory profile to match the observations. This 
is not true of the ground samples which are shown to exhibit 
considerable structure Ihin the 10 micron band. 
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Stephens (19S0) used the ssae spperstus to aessure the extinction 
of small, vapor condensed silicate, carbon and silicon carbide 
particles in the range 130nm < i < 800m. Although the silicate and 
carbon particles produced in this study were amorphous, silicon 
carbide condensed in the low temperature ( B) crystal form. The most 
interesting finding of this study was that amorphous or glassy carbon 
particles of mean radii I3nm and 6na produced extinction peaks at 250 
and 235 nm respectively. Graf^itic carbon (or possibly carbynes) must 
therefore be responsible for the observed extinction peak at 220nm. 

Day (1979) has reactively sputtered cathodes of I^Si and Mg2Si 
compositions in an Ar/O^ atmosphere to produce amorphous films vOf 
^^SiO^ and Hg2^0|^ composition respectively. He measured the optical 
constants of these films in the infrared (7 microns < X < 33 microns) 
and used these constants to calculate the extinction of Rayleigh 
spheres of the same nominal composition. Both compounds exhibit broad 
extinction peaks centered near 10 and 19 microns which compare 
favorably with those found in astronomical sources. 

A similar study was done on the properties of amorphous iron 
silicates produced by the reactive sputtering ol FeSi and Fe2i>i 
cathodes in Ar/O^ atmospheres (Day, 1981). Although the peak 
positions for the FeBiO^ and Fe28iu^ films are essentially the same as 
in their magnesim counterparts, the calculated extinction 
efficiencies are 30S to 75t lower. Measurements of the far infrared 
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tranaplission of both aaorphoua nagnaalua and iron allloataa indioataa 
that tha axtinotion falla off aa 


d. anokts 

Kaaijo at al. (1975) produoad fina particlaa of carbon, iron, 
iron oxide and silica, tha naan diamatar of which rart((ad batwaan 
4-100ru. Iron and iron oxide condensates exhibited aone crystal 
structure while carbon and silica anokes were amorphous. In all cases 
however, the small particles were observed to form chainlika 
'necklace* structures similar to those found in the laser condensed 
aookes of Stephens and Kothari (I97d). 

Day (1975) produced v<50<-60nm diameter forsterite spheres by 
striking an arc between hollow graphite electrodes containing natural 
fosterite. This method was later used in the study of Kritchmer and 
Hufilnan (1979) mentioned previously. Although the initial condensates 
were glassy, annealing at 800C converted the particles to crystalline 
olivine. The optical constants of this material were then measured in 
the range (9 microns < i < 25 microns). These spheres were probably 
not true aookes - vapor phase condensates - but rather were more 
likely soall molten fragments thrown from the electrodes and rapidly 
quenched by the surrounding gas (Day, 1976). 

A similar technique - evaporation in an electric arc - has been 
used by Lefevre (1970) to produce smokes of iron, carbon, silicon 
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oarbidt aiid silloa, Qrain aiiaa In thla atudy wara In tha rang# ( 10 m 
< r < 200na) . Cptioal abaorption maabrasanta of tbaaa particlaa wara 
aada and omparad with Hia oaloulationa for infinita oyllndara in tha 
ranga 360na < X < 700na. Good agraanant waa obtalnad only for allioa 
- poasibly indioativa of tha naad for battar optical oonatanta for tha 
othar aatariala. in all oaaas the partiolaa again foniad ohainlika 
'naoklaoa' atruoturaa aiailar to thoaa aantionad pravioualy. 

Day and Donn (1978a,b) deaoribe the production of Si^O^ (and 
aaorphous oagnaaiui silicate) aaoka by tha evaporation of SiO (plus 
nagnasiuB) into Ar, or Ar/H^O atooapheres at low pressures (P < 10 
torr) . In this study the condensate was shown to condense directly 
fron the gas phase and was extremely sensitive to the temperature of 
the ambient gas; i.e. no condensate formed at T s 500C at SiO-Mg 
evaporation rates sufficient to produce copious quantities of smoke at 
room temperature. The infrared spactriai of the amorphous magnasiin 
silicate showed some similarity to the observed 10 and 20 micron 
feature of astronomical sources and annealed (at 1250C) to fora 
crystalline forsterite. 

Althoii^h the experimental conditions under which particle 
formation occurred in this system were not well defined, these results 
were significant since they demonstrated the non-equilibriua nature of 
the condensation process in such refractory systems. Ihis was evident 
both from the amorphous nature of the condensates and from the drastic 
inhibition of particle formation at &00C in a system predicted 
thermodynamically to condense at .riiiooC (Larimer and Grossman, 197k) 



from a vapor lowar in oonoantratlon than that atudiad. lha 
oondanaation of and ayatoaa wart invaatigatad «ora 

thoroMghiy for thia diaaartatlon (aaa Chaptara 5 & 6). All of tha 
qualitativa findinga of Day and Donn (1978a, b) iMra oonfinaad and put 
on a auoh aora quantitativa baaia. In addition, tha infrarad apaotra 
of thaaa oondanaataa wara ahoun to baar a aignifieant raaaablanoa to 
the 8-25 aioron spaotrim of OH 26.5 "*■ 0.6 aa tha aanplaa wara annaalad 
in vacuum at lOOUC for timaa ranging up to 167 houra. 


II. Huoleation Thaory 

the partial preaaure of a chemical species, in equillbriun with a 
liquid or solid of similar composition, is determined uniquely by the 
temperature of the system. If, however, no solid or liquid phase 
exists, then it is possible for the vapor to become supersaturated as 
the temperature of the system decreases (assuae that interactions with 
the walls are insignificant). Nucleation theory attempts to predict 
the rate of condensation as a function of the degree of 
supersaturation, ambient temperature and chemical composition of the 
cooling system, one of the earliest solutions to this problem is 
based upon thermodynamic principles and assuaes that the condensing 
clusters can be approximated by liquid drops (decker and DSring, 1935; 
Volmer and Weber, 1926). This approach is known as the becker-Dbring 
or Qassical theory of homogeneous nucleation (see Reiss, 1977, 
Abraham, 1974). 


Thift Motion Mill to dlvidod into thrM ouooootions. In ttio 
first* 1 Mill briofiy skotoh tho basic dorivation of olassioal 
nuoltation thoory. in tho aooond, l will ouUino soot of tho 
•odifioationa and altarnativos to olassioal thoory whioh havo boon 
suggoatod ovor tho yoara. In tho third* 1 will roviow an application 
of olassioal thoory (baaod upon a stoady atato approx iaat ion) in whioh 
tho tiao dopondont nucloation rato for a cooling gas is oaloulatod as 
a function of the cooling trajectory and tno inhoront proportios of 
tho oondonsato. 


a. The Classical Theory of Homogeneous Nucleation 

If we assuae that the bulk of the cooling* supersaturated gas is 
in the form of monomers and that these react to form a quasi steady 
state distribution of n-«ers, then the equilibrium number density of 
suoh clusters will be given by 

Mjj * k^exp(-sG,/kT) (1) 

where is the number density of monomers, is the number oensity 
of n>mers* &G* is the free energy of formation of an n-mor from 
monomers and T is the temperature of the system. 

If one assumes that the clusters are liquid drops* then the free 
energy of foraation of these drops from monomeric vapor will be the 
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sia of tMO toms, on« dopsndsnt on ttos voluss and tha othar on tha 
surfaea araa of tha olustar. Tha fraa anargy of this systaa is given 

SG. s <4ikTinS ■¥ 3 a (2) 

■ n 

by Equation (2). where S is the super saturation ratio, s^ is the 
surface araa of a cluster and o is the surface tension of the liquid. 
Since the cluster is assuaed to be spherical, its radius, r, and 

, 1/3 2 

surface area are given by r s [(3na)/(4*p} j and s^ x Swr 
respectively. These can be calculated if one knows the density of the 
cluster, p, and the mass of a aonoier, m. Equation (2) can therefore 
be written as 


aG, X -(4«r^)/(3v)kTtnS ♦ 4»r^o (3) 

where v is the voluae per aoncaer in the cluster. 

The aaxiaui in aG,, the activation energy barrier, aG,*, in the 
formation of clusters from monoaers and r , the radius of the 
•critical' cluster at which aG, is observed, are easily calculated by 
taking the derivative of Equation (3) with respect to r. These 
quanti'^ies are given by Equations (4) and (5) respectively. The 

aG,* = (l6io^v^)/(3[kUnS3^) (M) 

r„ = C2ov)/(kTtnS) (5) 

c 
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•^Ulbriun aunbtr of oritioal »lMd oluatoro it Atund t»y tubttiiuUAi 
Equation (t) into Equation (l). 

lha attndy atata rata of fbnnation of oritleal aiaad oluatara, J, 
ia tlia produot of a nuibar of faotora. Ihaaa inoluda ttaa aonoaar 

iaplngaaant rata par unit aurfaoa araa of tha oluatar, p, tha aurfaoa 

a 

araa of tha oritioai siaad oluatar, a , tha ntabar danaity of oritioal 
oluatara, , and tha rata at uhioh thaaa oritioal oluatara 
ra-avaporata . For aiaplioity, tha atiokii^ ooaffioiant of a aononar 
onto a oluatar ia uaually aaauaad to ha unity, lha nuoiaation rata ia 
tharafora givan by 

J s (2o/«a)^^^vN^2e*p[-(l6»oV)/(3lkT}^UnS>2)3 (6) 

whara it haa baan aaauiad that Z, tha rata at which ra-avaporation 
ocoura, ia givan by Equation (7) (Raiaa. 1977) • 

2 a (AC,*/t3ikTn^^J)'^^ (7) 

It can ba aean froa inapaction of Equation (6) that tha 
nuoiaation rata ia axtraaaly aanaitiva to tha aabiant taoparatura, tha 
oonoantration of aonooiara (both and S) and to tha aurfaca anargy of 
tha oluatar, o. Nuoiaation thaory aaauaaa that tha aurfaoa anargy of 
a oluatar ia aqual to tha aurfaoa anargy of tha bulk liquid (or aolid) 
and that v, tha voluaa of a nononar in tha oluatar, ia alao oqual to 
ita bulk value. Criticiau of thaaa aaaunptiona, aa well aa objectiona 
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ooaottrning tht proper foni of AC, • hovo Imm tho buis of virtually 
all of tho oontrovoray in this fiold. 

Ihoaa oojaetions, as util as aoat of the suggested aodifioations, 
uill be discussed in the next subseotions* Houever, it aust be noted 
that experiaental evidence obtained on systeas Which *equire 
relatively low supersatwations in order to nucleate, and which 
condense as liquids, are found to agree quite well with the 
predictions of classical nucleation theory (Adaason, 1976 ). 


b. Hod if i nations and Alternatives to Qassical Nucleation Theory 

A long standing debate in the literature concerns the 'proper' 
fora of aG,*. Lothe and pound (1962) feel that Equation (4) does not 
adequately account for the statistical aechanical contributions from 
the translational and rotational degrees of freedom of the aonoaers as 
they coabine to form n-aers. Their alternative fomulation, when 
applied to the nucleation of %«ater vapor, indicates that the 
nucleation rate should be 10 to 10 ' times faster than the rate 
predicted by classical theory (Lothe and Pound, 1962). 

Several authors (Neiss, 1977; blander and Katz, 1972) have argued 
that such factors are already accounted for by the liquid drop 
approximation. More important. Blander and Katz (1972) have shown 

ft 

that there is a mistake in the original derivation of AG, ar.d that 
Equation (3} should contain an additional term on the right hand side 
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(♦ KTuiS). This in tn additlonai factor of in the 

prt-cxponcntiai factor of fiquation (6). Ihia factor ia iK>t iaportant 
for the nuclcation of aost syataaB atudind to data ainct the required 
auperaaturationa have uaually been reiatively aeall. Howaver, for the 
nucleation of aore refractory aateriaie auoh aa SiO, FdO» etc., S can 

c 

he on the order of 10 and thia additional factor can therefore hecoae 
al^nificMt. 

Another aajor problea in the claaaioal foraulntion of nucleation 

theory Is the assuaption that quantities auoh as surface tension, 

density and &G, for bulk materials, can be applied 'u> aaall clusters 

(Nishioka and Pound, 1977; Tolaan, 1946. 1949; Abraham, 1974). Such 

17 

assumptions could lead to errors as large as 10 (Oriani and 

Sundquist, 1963) for the nucleation rate of water vapor and night 

17 

actually cancel the predicted increase of rio mentioned earlier 
(Lothe and Pound, 1962). Because of such uncertainties, Bauer and 
FTurip (1977) have proposed a model for the nucleation of highly 
supersaturated systems based upon a kinetic approach. They find that 
the rate determining 'step* occurs at a size greater than that of the 
critical cluster and that the overall nucleation rate is determined cy 
the rate at which AG«(r) decreases beyond the critical cluster. Their 
model predicts that refractory compounds tend to nucleate slower than 
the rate predicted by classical theory. 

An alternative, kinetic approach to nucleation, which utilizes 
quantities which can in principle be measured in an equilibrium 
system, has been formulated by Katz and Donohue (1979). It Is 



lattrttting to nott that tholr proaloto4 niwltotioa rott (atotfilni oa 
Idoal t**) Moatloal to that praOiotod by olaaaloal tiMory. Of 
oooraa auolaatloa, by 4aflaitloa« oaa atvar ooour ia aa idaal taa. 
Navarthalaaa, thalr fOraaUaa doaa offar tba poaalbUity of ualoi 
•aaaurad oluatar ooaoaatratioaa , obtainad uadar aquUlbriua 
ooaditioaa, to pradiot tha auoiaatloa rata la a auparaaturatad syataa. 
Whaa oaa raaaabars that thaoratioal rataa oaa vary by as auoh aa 10^** 
dapaadiiiA upoa tdUoh tama ara daaaad iaportaat, auoh a rallaaoa upoa 
an axpariaantaliy baaad prad lotion la Indaad an aaaat* This foraallaa 
has raoantly baan axtandad to traat tha oaaa of auoiaatloa for syataaa 
in which tha oondanalbia aonoaar la produoad by tha ohoaloal raaotlon 
of originally nonoondanaibla gas (or liquid) phasa oonponanta (Kata 
and Donohua, 19d2) • 

Finally, Giliaspio (1981) has shown that hoaoganaoua nucltrtion 
can t>a viawcd in tarns of a stoohastio prooass. In this nodal, tha 
aquiiibriua, pra-nuoiaation olustar distribution is dataminad by tha 
probability that an (n-D-nar will gain a nononar to baoona an a-nar, 
varsas tha oonbinad probabilitias that an n-nar will aithar losa or 
gain a nononar. It is shown that, with nuiarous raasonabla 
assunptions, this nodal raduoas to an axprassion for tha nuolaation 
rata idantioal to tnat pradiutad by olassioal thaory. In addition, 
tha nodal allows ona to caloulata tha 'tina lag* bafora the onset of 
nuolaation due to tha tina required to establish a population of 
oritioal nuclei. Again, it is in principle possible to naas^a tha 
various quantities necessary to pradiot tha nuolaation rata of a 
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partioulwr systai of intorott. Ihlt la o^Meially trua of a ayatw in 
which the alaa of the oritioal nuclaua la aipaotad to ha aaall. 


o. Tiaa Oapandant Nuoiaatlon Thaory 

Draina and Sal pa tar (1977) hava davlaad a nodal, baaad upon 
olasaical nuoiaatlon thaory, which can ba uaad to pradlot tha 
nuoiaatlon rate, tha avaraga slsa of tha final clustars and tha 
diaper alon about thia naan for a cooling, auparaaturatad vapor, lhaaa 
quantitiaa are calculated in tama of two dlnanalonlaaa parnaatara, e 
and rt* *3 wall aa other rather atrangaly dlnanalonad quantitiaa. 

Theaa additional quantitiaa are uaad to daacriba tha cooling 
trajectory of tha gaa and tha thamal >ropartlaa of tha oondanaing 
cluatera. e Is defined in auch a way that k 6 || Is aasantlally equal to 
tha aurfaca free energy per aurfaoa alta for a oluatar of N nononara. 
n ia defined aa tha nuabar of tinea an avaraga nononer can be expected 
to collide with (and atick to) an average aurfaca aite during one 
auperaaturation ratio e-folding tine on the cooling trajectory of the 
gaa. 


Ihe derivation of thia nodal ia extrenely ccnplex - and ia 
further obacured by the uae of a hoat of paranetera apeclfically 
defined in tema of auaa and ratioa of nore faniliar quantitiaa. Thia 
model alao relies heavily upon classical nucleatlon theory and is 
therefore subject to most of the sane problens mentioned in subsection 
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Drain# ana Salpatar (1977) find that tha 'orltioal* 
suparaaturatlon at whioh nuolaation ocoura from tha vapor ia givan by 

inS^ « a,“Uu(b,/B,')(tnA/8ti»n)-(8inn)”^) (8) 

where B, i (2T/e)^^^(tnn) B,' * B«-3c/2* * « 
t(Tf-1)"'-a/hT]"\<Un(e/T)/dtnTJ, B s icT(din{n^^)/dtnT) , y s 
1<»(dinT/din[nl) , n ia the total denalty of oondenaible naterial, 
a* (2^^/3)(i/3)'^^(B,'/B,)‘*(Ve,)^(inn)^^^/U<f^). X s 
c^(1-X)~^n~\ X la the uaaa fraction of condenaible Material, la 
the concentration of aoncmera and <f^> • 1-exp(B,)> 

Of theae quantities, B can be determined experlnentally; y la 
deteralned by the asauxed cooling trajectory and 'Ideality* of the 
gas. The dependence of e on T Is virtually anyone's guess as la the 
dependence of e on N* The concentration of nonoaers la tiiae dependent 
and n and X are arbitrary Input parameters. Many of these terms are. 
of course, redundant - not to mention confusing. Nevertheless, the 
average size of the final cluster distribution is given by 

<Nf^nal^ - (2/t)^^^(3e/2®)(B,/B,')^<f^>(0/T,)^(n/inn)-^ (9) 

and is thus quite sensitive to both e and n since Bt* Bt* t and <f^> 
also each depend on these quantities. 
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Tht r«lativ« dlsptrston in the aean cluster siss dlstrlhutloo •> 
in tsnis of the nusbtr of aonosers inoorporsted into tho cluster - is 
given by 

s 0.28(e/T)^^^(tnn)’^^^ (10) 

where the fsctor 0.26 was aapirically chosen to aske the equation 
agree with the results of nuaerical calculations. The application of 
this theory to the condensation of grains in the ataospherea of stars 

j 

will be discussed in section 111 below. j 

1 


111. Astro physical Applications of Nucleation Theory | 

1 

J 

i 

f 

# 

Nucleation theory has been used to calculate the conditions under | 

which graphite grains aight condense in the atmospheres of cool. I 

I 

carbon rich stars (Donn et al., 1968; Fix, 1969). Sal peter (197Sa) 
used a fora of nucleation theory, modified to account for the 
temperature increase of the clusters because of the high latent heat 
of condensation of crystalline species, to treat the condensation 
process in both oxygen and carbon rich stars. He then predicted the 
conditions under which either dust alone, or a mixture of both gas and 
dust Bight be ejected from the atmospheres of s^ch stairs (Salpeter, 

1974 b) . 
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An smary of the problMS enootmtered in the 

epplioation of nuoleation theory to oondensetion prooessee in stellar 
stoospheres is given by Tsbek at al.(1975). lYiese authors 
specifically discuss the effects caused by the assuaption of bulk 
parameters (gG, y, crystal structure) in such calculations. These 
authors also discuss the result of using the Lothe-Pound formula 
rather than the classical model. They show that the nucleation of 
carbon grains is extremely sensitive to the value adopted for y* They 
also note that the value for the surface free energy, y, which might 
be appropriate for graphite clusters, is a matter of considerable 
uncertainty. 

Yamamoto and Hasegawa (1977) ohereafter YH- independently 
formwiated a time dependent model of nucleation to calculate the size 
distribution of the grains formed in astrophysical systems. In 
general, their formalism is quite similar to that of Draine and 
Sal peter (1977) - hereafter 06 - in that nucleation is found to depend 
sensitively on two parameters, y and a, which are analogous to e and n 
of 06. The application of YH theory is made much easier due to the 
inclusion of some extremely useful tables; one of which tabulates r . 

m 

the radius of the finax clusters - as a function of y and log a is 
particularly valuable. 

YH apply their model to calculate the condensation sequence of a 
gas of solar composition. They obtain qualitative verification of the 
conclusions derived by blander and Katz (1967), i.e. that nuci.eation 
barriers can effectively disrupt the thermodynamic condensation 


92 



Mqutnot propostd by Laria«r (1967) or Grosnan (1972). Of apaolal 
inttreat. 7H art oaraful to note aoat of the poaaibie probleaa whloh 
oQuld alter their oonoluaiona. Such probleaa include unoertaintiea in 
effeota of the oondenaation of apeoiea fonaed by oheoiioal reaotiona 
in the aas phaae end poaaibie effeota of heterogeneous nuoleation on 
the aurfecea of 'pre-existing' grains. In the case of nuoleation in 
stellar ataospherea, these would be the grains which condense first 
and could therefore lower the condensation barrier for less refractory 
materials. 

Oeguchi (1980) has used tne combined YH-Di nuoleation model to 
study the effect of grain formation on the gas phase chemical 
composition of cool clrcuasteilar envelopes. He shows that although 
the YH-D2S model is based upon classical nuoleation theory, inclusion 
of the correction factor proposed by Lothe and Pound ( 1962 ; only 
changes the total number of predicted grain nuclei by a factor between 
8 and 64. IXie to the other uncertainties inherent in the analysis, 
this disagreement is unimportant. 

Deguchi (i960) finds that -/ 90 -yyt of the metals in tne ejected 
gas condense into grains when the mass loss rate is eo’^al to 10 
solar masses per year. More than 991 of the refractories condense at 
mass loss rates less than j x 10”^ solar masses per year. He also 
calculates tnat nuoleation occurs at v4|-6 times tne stellar radius. 
Uiven the uncertainties in the model, these conclusions are in 
reasonable agreement with the calculations of Hagen (1978). However, 
this may be fortuitous since Hagen (1978) concludes that grain 


9j 


foraiatlon ooours «s a oonaaquanoa of atallar aaaa loaa whila Daguolii 
(1980) finds that nuoiaation drivas it in stars haaviar than 1.5 solar 
aassas. 

Drains (1979) has appliad tha tima dapandant nuolaation nodal of 
OS to tha oondansation of graphita in carbon rich pro to pi ana tar y 
nabulaa and novae, altiiough tna theory successfully predicts the 
formation of graphite in such sources, it cannot fona silicate grains 
in the outflow from normal giants and supargiants without the 
invocation of density enhanoeiaents on the order of factors of 100 over 
the densities derived using the assmaption of a steady flow. 

In a later paper (Oraine, 1981), condensation in such outflows is 
snown to be possible due to the fact that both the 810 monomer and 
small (10 atom) 'olivine* grains can be characterized by temperatures 
considerably less than that characteristic of the translational 
velocity of the gas. Even with these approximations however, the size 
of the critical clusters tnought to nucleate in such environments is 
barely larger than that for which the D8 model begins to fail. In 
addition, one can only speculate about the properties of a 10 atom 
'olivine' cluster since olivine 'monomers' already require at least 7 
atoms. Condensation processes in such systems must surely oe 
influenced oy the reaction rate oi' HiO (or f-eO) ♦ 810 molecules to 
form critical clusters in the gas phase. This problem is discussed in 
somewhat more detail in Ctiapter 6. 
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L«f«vr« ( 1979 ) hat ahotn that tha growth of graina whloh nuolaata 
in stailar aouroaa la vary aanaltiva to tha taaparatura of tha 
cluatars and that thla taaparatura ia, in turn, likaly to ha dlffarant 
than that of tha aahlant gaa. Thia raauit la not dapandant upon tha 
foraatlon rata of orltloai oluatara but rathar oonaidara only tha 
aubaaquant growth of thaaa oluatara Into aacroaooplo graina. Ha 
conoludaa that tha auparaaturation ratio uaad in thaoriaa appiiad to 
aatrophyaioai ragiona ahould ba dafinad in taraa of tha aquilibriuM 
taaparatura of the critical cluater rathar than that of the aabiant 
gaa. Thia. of oouraa. requiraa that one be able to calculate tha 
taiaperaturea of aonoaera, dlfflera,...n-aera....&«all graina. Thia 
problaa is discussed in aore detail in Chapter 3. 


IV. Nucleation bcperiaents at High Supersatirations 

Very little experiaental work has been done on the condensation 
of refractory materials , although the nucleation of aore volatile 
species - i.e. water - has been studied for decades. Such work formed 
the experiaental basis upon which classical nucleation theory was 
ouilt. It is only aore recently that the vapor phase nucleation of 
refractory species becaae interesting. This is due ir. large measure 
to the importance of such processes in the formation of particulates 
in tne combustion of coal and other fuels, and the desire to control 
pollution from such sources. Refractory aerosols are also produced by 
the recondensation of meteoric materials in the upper atmosphere and 
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in niMinroua induatrlal prooasMs. It atoiiid oovioua fron t^.• 
pravioua Motions that thoM prooaaMa art also important in 
aatrophyaioa , 

Oorftid and Hudaon (1973) attaoiptad to «xpt>*iBcn tally iotarmine 
tha dagraa of graphita oondanMtion Mhich night ba axpaotad in tha 
axpanding anvalopas of cool, carbon rich atara. This waa dona by 
allowing hot (^OOOK) nixturaa of nethana and hydrogan to axpand 
through a auiall orifioa into a vacuin. lha final CH^/H^ ratio was 
daterninoJ by nasa apactrosoopy and oonparad with tha initial ratio. 

A decraaaa waa tawan to be indicative of the foraation of grai^ite in 
the expansion. It should be noted that the foraation of graphite 
could not be directly obaerved in this systra. me expansion waa 
scaled so that the number of temiolecular collisions in the experiment 
waa equal to the nuaber expected over the timescale of the stellar 
outflow. Ihey concluded that little or no graphite condensation could 
be expected in such sources and obtained an upper limit of 10H of the 
available carbon condensed into grains. 

mese results might be due to the short timescales over which the 
reactions were followed in the laboratory as opposed to those 
available in stellar sources, mey might also be explained by the 
fact that, in this experiment, carbon atoms were not readily available 
but had to be produced by the dissociation of CH^. However, the 
dissociation of methane should become increasingly difficult as tne 
gas expands due to the decrease in tne average energy available to 
break the C-»H bonds as the gas cools. An equilibrlua carbon atom 
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oonoMtration mb probaoly n«v«r approBOhed in this systM snd thair 
oonciualons are therefore at least suspect. 

Bauer aiKl eoMorkers have studied the nuoleation of supersaturated 
aetai vapors as a function of ambient temperature in a series of shook 
t^ experiments. Freund and Bsuer (1977) experimentally determined 
the binding energy per atom for the reaction nFe ♦ Fe^. They found 
that the expression sE^/h s aE^(l-n adequately represented their 

data. is the siAllmatlon energy of the bulk metal. This result 
is interesting in that It demonstrates that aoaoaers are bound much 
less strongly to small clusters than to the bu k metal. 

Frurlp and bauer (1977a) determined the temperature dependence of 
the critical supersaturation ratio for iron, lead and bismuth. They 
showed that, even allowing for a substantial degree of experimental 
error, neither the classical nor the Lothe-Pound theory could be used 
to predict their results, m addition, due to the probable dependence 
of the surface tension on temperature, they questioned if such 
measurements oovld reasonably be used to distinguish between these 
theories. 

Frurlp and Bauer (1977b) developed a tecrmique by which an 
independent determination of the nucleation flux could be made. They 
then used this technique and tneir previous results on the temperature 
dependence of the condensation of lead, bianuth ana Iron to develope a 
kinetic theory of nucleation (Bauer and Frurip, 1977). This theory 
completely abandoned the liquid drop approximation and the concept of 


97 


a 'oonatralnad at^Uibrlui diatribution* , un for tuna taly, tbla thaory 
la axtroMaiy uiffloult (if not Inpoaalbit) to apply to Most ayataaa of 
astropbyaioai Intarast ainoa it raqulraa tha input of datailad 
raaotion rata oonatanta for all apaciaa oV iaportanoa* Suoh datailad 
infomation ia ob>^ioualy not availabla at tha praaant ainoa avan tha 
chaaioal ooapoaitiun of aooia ^rain fonain^ ragiona ia atill a aattar 
of dabata (2uckanian« 1980). 

Racantly, :$taphena and Biuar (1981) hava raportad aaaaurauianta of 
tha critical praaaura for tha onaat of avalancha nuclaation aa a 
function of acnbiant tacaparatura for tha condanaation of Fa, 81, Fa/8i, 
FaO^ and *11 caaaa thara was substantial disagraaiaant 

between tha maasuraoiants and pradiotions basad on either tha classical 
or the Lotha-Pdund formulations of nuclaation theory. 

8upar saturations observed in these shook tube experiments ranged from 
v>iu for iron to several orders of magnitude for tha 8iU^ systaia. 

A comparison of the snock tvkie data of Stephens and bauer (1981), 
ootainad at temperatures between 12!>0K and U200K, and that obtained in 
a quasi steady state system at temperatures between 7!>0i( and lUJQK is 
presented in Qiapter This chapter reports measuraaents of the 
,emperature dependence of tna onset of avalanctie nucleation in an 
8iQ-H^ system as a function of teiaperature in a system similar to that 
used in the work of i;uy and Dunn (I9^8a,b). When the magnitutle of the 
possible error in both experiments is canbined with the fact that the 
tetaperature ranges of the measurements do not overlap, a detailed 
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ooipitrlaon of thcM two oota of roaulta ■awiinAloaa, Howtvor, tht 
ovorail airaaaiant botwaon tbtao difforant aatlioda appaara to ba quita 
good. 


V. buAiaary 

It baa bean Ueisonatratad (Kritcbmar and Huffaan, 1977: 0<*yi 19 / 9 « 
1961) that an amorpboua iron or ni^naaiua sliioata can produce 
infrared extinction (6 aicrona < X < niorona) which is quite 
similar to that observed in obscured, oxygen rich, stellar sources. 

Ihe spectra of more crystalline materials Joes not reproduce such 
Observations. Host theories of nucleation rely upon the attainment of 
equilibrium conditions (at least before the onset of nucleatxon) and 
predict the fonaation of a thermodynamically stable (in most cases a 
cr ystal line) sol id . 

The results of nucleation experiments on retVactory materials, 
similar to those whicn might condense in astronomical environments, 
deviate considerably from the predictions of both the classical and 
Lothe-Pound formulations of nucleation theory. It must be emphasized 
that tnese experiments were performed in terrestrial laboratories 
Where the Kinetic temperature of the gas, the vibrational temperature 
oi the clusters and the surface temperatures of tne grains were 
virtually identical. This is not likely to be the case in most of the 
regions currently believed to be sources of fresh interstellar grains. 
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This point it ditousatd in tort dtttil in Ghtpttr 3 


Exptrittntai rttultt tor tho avtltnoht oondtnattion of SiO-H^ and 
Hg-SiO-H^ ayataaa as a funotion of tho aabiont gaa taaporaturo art 
roportod in Qiaptora 5 and 6 rospootivoiy. 
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Chapter 3« VlbretionaX PlaeQuiUhrluM in Retiona of Grein Foritlon 


It is well known that in diffuse interstellar olouds the internal 
energy of atons and noleoules shows extreae deviations fron thenial 
equilibriuB. Donn (1975. 1979> pointed out that significant 
deviations for molecular vibrational energy levels would occur in the 
much denser regions where grain formation has been postulated. These 
regions include novae and supernovae ejecta, the shells surrounding 
cool stars, planetary nebulae and several models of the primordial 
solar nebula. An approximate calculation by Thompson ( 1973) indicated 
that nonequilibrium vibrational distributions are possible in the 
atmosptieres of cool supergiants. Unfortunately, the criterion adopted 
by Thompson (1973) for the onset of non-LTE in astrophysical systems 
lacks precision. Further, it does not permit the quantitative 
analysis necessary to facilitate improved calculations of molecular 
abundances or to gain a true understanding of the nucleation process 
in such systems. In addition, his analysis suffers from the adoption 
of incorrect relaxation efficiencies, (ieballe and Iriollman 

( 1972 ) have observed an apparent non-equilibriun temperature 
distribution in the CO spectra of cool stars. Hore thorough model 
calculations by Carbon et al. (1976) have shown that such effects 
might indeed be expected. 

In this ctiapter we first calculate the vibrational distribution 
and corresponding vibrational temperatures for the diatomic molecules 
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oarbon nonox Ida and ail loon nonox Ida ovar a wlda ranga of gas klnatio 
tanparaturaa and praaauras for aavaral valuaa of a diluted blaok body 
radiation fiald. Wa than astinata the vibrational tanparature of an 
idealized polyatonio molecule under ainilar conditions. Carbon 
nonoxide vms chosen because of its abundance in nost regions of 
interest and silicon monoxide because it plays a major role in the 
condensation of nunerous predicted grain materials. Most of the 
necessary nolecular data are available for both. 

The factors reponsible for the disequilibrim at lou pressures 
are (1) low collision frequency, (2) general inefficiency of 
translational-vibrational energy exchange during collision 
(Kondrat'ev, 1964; Cottrell and McCoubrey, 1901) and (i) radiative 
lifetimes for vibrational transitions, which are generally in the 
range 10"^ - 10"^ s (Penner, 1959). The net effect is to make 
radiative decay comparable to or faster than collislonal excitation. 
Consequently, at low pressures the vibrational energy distribution 
will be seriously perturbed from an equilibrium tk>lt 2 iaann distribution 
by the radiative depletion of higher levels. 

This effect cah have a sighificaht ihfluence oh a humber of 
physico-chemical processes in low density cosmic clouds as is shown by 
the quantitive calculations of Dalgarno and Roberge (1979). This 
effect is of particular significance for molecular abundance 
distributions and condensation processes. It raises serious questions 
about analyses of these phenomenu using the assumption of thermal 
equilibriuo. 


lUb 


II, Method of AnalyaU 


a. Model 

As noted by Thompson (1973)* at tempei atures below 
electronic transitions are rare and do not significantly contribute to 
trie vibrational population of CO. Furthermore, condensation of grains 
is expected to occur at temperatures below 2000K in regions idiere the 
radiative energy density has been considerably diluted. At such 
temperatures, the bulk of both CO and SiU occupy low vibrational 
states even at equilibrium. We treat both systems as multilevel 
molecules and solve for the populations of the lowest twelve levels at 
which point more than 99.99999)^ of the population has usually been 
accounted for. 


A "vibrational temperature" can be defined by the ratio of the 
populations of two states sucn that 


N 

t 


-hcu 
exp — — 
kT 

V 


( 1 ) 


Where and are the populations of the upper and lower states, 
is the transition energy in cm'^ h, c, and k have their usual 
meanings and is the viorational temperature. In equation (1) tne 
statistical weights have been set equal to one as for the case of a 
harmonic oscillator. Tnis is a good approximation for the low 
vibrational levels vmich are Important here. 
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An txaot oaloulatlon of tht population of vibratimal lavala 
raquiraa knouladgt of tha radiation fiald, oollisional oroaa aaotiona 
for tha axoltatlon of individual lavala by all colliding apaoiaa, aa 
wail aa tha radiativa tranaition prooabiiitiaa froa oocupiad lavala. 

In addition, tha ohaaioal oaapoaition, taaparatura and praaaura of tha 
region auat be known. Baoauaa of thaaa raatrictiona an axaot 
calculation for any particular region cannot be carried out. However, 
a aolutlon utilixing a ''iw raaaonable aaauaptions and aoaawhat 
Idealixad oonditiona can be expected to yield raaulta sufficiently 
accurate for the treatment of problems in which the vibrational 
populations beceme important. 

Vibrational equilibriua for harmonic oscilletors has received a 
very general treatment by Rubin and Schuler (1956a, b; 1957). Iheir 
analysis is a good approximation for the low vibrational levels of 
molecules such as CO and SiO. We adopt their n level system and 
introduce diluted black body radiation through a dilution constant, W. 
Ue then determine the steady state vibrational distribution rather 
than the relaxation rate which was of interest in the much higher 
pressure systems of Rubin and Schuler and in which equilibriua is 
ultimately attained. 

For n level systems we obtain the following set of differential 
equations fur the steady state population of a given level, with the 
approximation that oollisional and radiative transitions occur only 
between adjacent levels and all energy level spacings are the same: 


10S 


( 2 ) 




^^10 ■*“ *10 * ^ 10 ^ 


^^01 * '^01^ * ® 


* *X-1^S-1,x'*’^X-l,X^ * **x^1 ^^'x^-l.x * ^X*1,X ■*■ *X+1,X^ 


" **X^S,X+1 '*“ ^X,X>fl * ^X,X-1 ■*” ^x,x-l '*’ *x,x-l^ * ^ 


t4ier« is the population of level x and and represent 

the probability of collisional, spontaneous radiative and induced 
radiative transitions from level x to level y, respectively. 


We use the following relationships (Rubin and Schuler, 1957). 


J-J*1 • “=>10 "P 

^♦l.j 

s (j ♦ 1) C 

J^l.j 

Vi.j 

s B 

= olexp(|y -1)J”' s aQ 

a s dv 

h«^ • * , 
J .J+ 1 


1.0 


Vl.J * *1.0 


where A c^nd b are the Einstein coefficients, and is the black body 
radiation density at the energy of the transition of interest. 
Equations (2) and (3) are solved for population ratios and yield: 
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li„ (x)(C,„ •xp(-hCtt/kT) ♦H0 A^q) 4- I x^l )( C,„+(U WQ)A,n) 


'10 
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(x+1)(C^Q exp(- how/kT) ♦ q) ♦ (xHC^^ ♦ (1 ■*• WQ)A^q) 


( 5 ) 


where T is both the kinetic and color temperatu*e of the system. This 
approximation for T will be considered further In section V. 

b. Molecular Constants 


The Einstein A's for SIO were taken from the work of Hedeiund 
find Lambert (1972). For CO, A was estimated using the 

j+ 1 , j 

approximation 

A . s (J-t- 1)A_ where A.,, was calculated from Cnakarian (19Y0): 

J f I Iw 1C* 

^10 * zJ (TM^) (6) 

-1 ^ 

Where I'M is the transition moment (for Co, TH s 1.04 x 10 ) and Z is 

a constant which depends on the particular vibration ~ rotation 
transition of interest. For the P branch Z a (J^ ♦ 1)/(2J^ ♦ 1) 
whereas for the R branch Z a Jy/(2Jy ♦ 1) where is the rotational 
quantun nuaber. :>lnce we are concerned only with the total nuaber of 
transitions per vibrational level we sun the K and P branch transition 
rates and Z tneretore equals 1. Equation (6) was also used to 
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oaloulate values of for our abdel polyatoalc aoleoula using 
reasonable values for m and IM. 

for* CU, SiO and the aodel polyatoaic aolecule was initially 
calculated frua the work of Millikan and idiite (1963) as adopted by 
Ihoapson (1973) for a haraonic oscillator (see Herzfield and Litovitz, 
1959). 



where P is the total pressure in atmospheres, |i is the reduced aass of 
the colliding system in units of molecular weight, e is the energy of 
the transition in units of and T is the kinetic temperature. The 

expression from which this formula is derived, an empirical relation 
fitted to a large number of experimental relaxation measurements >jnder 
a very wide range of parameters (i.e., T, P, it. 9) is good to better 
than a factor of 2 in most cases (Millikan and White, 1963). 

The formalism of Millikan and White 096j) does not apply to the 
relaxation of diatomic species by atoms with open electronic shells. 
These atoRMOolecule collisions have been shown to be much more 
effective than previously predicted In transferring energy from the 
translational to vibrational degrees of freedom (von hosenberg, Taylor 
and Teare, 1971; West. Weston and Flynn, 1977; MacDonald and Moore, 
>97d). Nikitin (1974) has proposed that Stark splitti.ng of degenerate 
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•iMtronio spln-ort»it ooMpoiMfita oan laad to roaonant oolltaiona ami 
thus to aora affiolant translational to vibrational energy exchange. 

Mon-ad iabadie electronic energy transfer will alao be iaportant 
in colliding systeas where a stable reactive coaplex can be fonaed 
(see the review by anith, 1976). In this case, sv can be greater than 
♦ 1. Carbon et al. (1976) have desionstrated that these energy 
exchai^e aechaniaBS can have a significant effect on the calculated 
deviation of CO froa local thenaodynaaic equilibriua in the 
ataospheres of cool stars* These effects will be discussed in section 
V. 


3 

Because for H atoa-CO collisions is at aost a factor of 10 
higher than predicted by Hilllkan and W)ite (1963) (see Figw*e 1 of 
Carbon et al. (1976)) we have for simplicity used a collisional rate 
(C^q) which is a factor of 10^ higher than that calculated from their 
equations for Uils system. Similarly, because Fe atoms also very 
efficiently excite CO vibrational levels (von Rosenberg and Wray. 

1971 ) and because this discrepancy is also of order 10^. for iron we 
have again used a collisional rate which is 10^ times that calculated 
from the work of Millikan and White (1963). These approximations 
should underestimate the disequilibriia of the system. C^^ for 0 
atom-Co collisions was derived from the work of Center (1973). 
for H atom-SiO and CO-CO collisions were used as calculated from the 
equations of Hilllkan and White (1963). Although the use of these C^^ 
values probably overestimates the degree of disequilibriun for H 
atom-310 collisions, they should be correct for the CO-CO system. 
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The collisional relaxation rate. is proportional to the 
product of the pressure and the collisional efficiency. An increased 
efficiency coabined with a proportionate decrease in pressure will 
leave the systea laichanged. Therefore in the Figures the effect of a 
different choice of collisional efficiency is to aultiply the pressua*e 
scale by a factor equal to the ratio of the original efficiency to the 
new efficiency. 

Using equations (1), (4), and (3). the vibrational tenperature of 
a particular excited state with respect to the ground state can be 
calculated . 

III. Disequilibrium of Piatomic Molecules 

Figure 1 is a plot of the vibrational temperature of CO as a 
function of H-atom pressure at a kinetic temperature of 2000K for 
various values of the radiative dilution coefficient W, (W s 1 for a 
black body radiation field. W = 0 in the absence of radiation). It is 
obvious from the figure that althoi«h the radiation field provides a 
lower limit to the vibrational temperature it has little efffect on 
the pressure at which the system begins to depart from equilibriun. 

It snould be noted that there is little difference between the degree 
ot disequilibriuD for W equal to zero and w equal to a few tenths 
until the vibrational temperature is well below tne kinetic 
temperatire . 
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Figure 1. Vibrational tempcratawe of CO as a function of total 


pressure in an H-aton gas at a kinetic tesperature of 2000 K for 
various values of the radiative dilution coefficient, W. (W s o in 
the absence of a radiation field.) 









It should be noted that for ^undent noleoular species suoh as CO 
the optical depth at the wavelength of the fundanental frequency could 
be considerably higher tnan average due to the trapping of resonance 
radiation by the CO molecules. In most cases of Interest however, the 
density of molecular species should be low enough that this effect 
does not become important. 

Another Interesting feature is that under any particular set of 
conditions, all vibrational level populations can be represented by a 
Boltzman distribution at a characteristic vibrational temperature. 

This results from the fact that our model for diatomic molecules has 
an analytical solution which is independent of the vibrational level; 



KT ^ 0^*1 0 * 

^ ^ Kcs “r:;"-7A“~r5riJprteS7kT7 


Unfortunately, this equation requires that both the Einstein A 
coefficients and the coilisionai efficiencies for the higher 
vibrational levels be given by the idealized formulae of Kubin and 
Schuler (1^7). This precludes the use of more accurate results, sue 
as the calculation of the Einstein A coefficients for SiU by Hedeluni 
and Lambert (197^). Although equation (4) assumes the same idealized 
dependence of the C and A coefficients for higher vibrational levels 
as does (d), it is a trivial task to separate the variables in such a 
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way that aquation (4) could aoooaodata battar data aa it bacaaas 
availabla. This is not trua of aquation (d)> 

Thompson (197i) stated that non-LTE affects become important whan 
the rate of downward collisional transitions equals that of 
spontaneous radiative transitions. We find from equations (1). (4) 
and (5) that when these races are equal at a kinetic temperature of 
2000 K and a dilution coefficient, W, of 0.5, the vibrational 
temperature of CO is 1705 K. In the absence of a radiation field (W s 
0) the above conditions yield CO vibration temperatures of 1380 K. 

When Thompson's criterion is satisfied the system is already far 
removed from equilibriin. Further, his criterion provides only a 
qualitative measure of disequilibrium and cannot be used to determine 
the consequences of the disequilibrium upon other phenomena, for 
example, molecular equilibrium or particle condensation. 

Figure 2 is a plot of the average vibrational temperature of CO 
as a function of H-atom pressure for various values of the gas kinetic 
temperature in the absence of radiation. Figure 3 is the analogous 
plot for ^10. Because of the use of classical collisional 
efficiencies for H-BiO, the BiO pressure scales will be too high. If 
this system benaves as does H>C0, the correction factor may be as high 
as 1000. As there are no experimental measurements for 
collisions, the classical results were retained. 
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Figure 2. CO vibrational teaperature aa a function of total pressure 
in an H atom gas. Each line represents a calculation at the gas 
kinetic temperature of the intercept, ri x 0. 






Figure 3. SIO vibrational temperature as a function of total presswe 
in an H-atom gas. Each line represents a calculation at the gas 
kinetic temperature of the intercept, U s 0. 



3000 






IV. Piatquillbriw of Polyitoalc HoltoulM 


Th« axtraae vibrational diaaquilibriua found for diatonic 
moleculaa in clouda and ita oonaaquenoea for diaaooiation (Daigarno 
and Robarge, 1979) oalla for an invaatigation of tha vibrational 
behavior of polyatonic nolaoulaa. lhaaa sore oonplax apaoiaa with 
nany vibrational nodea cannot ba analysed as accw*ataly as can 
diatcaic aolaculas (HoDonald, 1979). Lnbart (1962, 1977) reviewed 
intranolecular vibrational energy conversion. Generally, the 
relaxation process can be characterized by a single rate, indicating 
rapid colllsional internal energy transfer, if the energy of the more 
energetic fundamental is less than twice that of the lesser. 

Otherwise each mode has a separate relaxation rate. 

In order to carry out an analysis for polyatomic molecules 
similar to that for diatomics we used a hypothetical molecule with two 
fundamental modes. In our initial calculations the energy of mode 2 
was assuned to equal three times the energy of mode 1. We then assume 
that collislonally induced energy exchange between modes is relatively 
rapid compared to other relaxation processes (Figure 4). This is 
based in part on Lambert's discussion cited earlier and in part on the 
effectiveness of internal energy transfer when long lived complexes 
are formed (McDonald, 1979). Ihe latter appears possible in 
atom-molecule collisions. We also assune that our hypothetical 
molecules are not excited in both vibrational modes simultaneously. 

We therefore ignore all transitions which would either originate or 
end in such a state. Ihls greatly simplifies the Kinetic equations 
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Figure 4. Schmatio representation of transfer processes in 
nypottietical dual aode poiyatoaic nodeis witn a 3:1 and 3^2 ratio of 
vibrational fundMaental energies. 


m 



MODEL POLYATOMIC MOLECULES 



Rate of A » B > C Rate of A » D > B > C 
No off resonance Off resonance transfer 

transfer allowed allowed (D) 


B and C contain both collisional and radiative terms 
A and D contain only collisional terms 


uhich deaerlb* Uit steady state population distribution and allows us 
to crest separately the high and low energy sodas in a sanner which is 
analogous to our treatsent of diatonic solecules in equation (3) by 
the addition of terms of tiM fora M (G, „) where G, „ is the rate 

of transfer oetween resonant levels of sodas 1 and 2 and the Kronecker 
delta allows transfer between nodes only at resonance, luc validity 
of this assuiption will be discussed at the end of this section after 
the model and its results have been presented. 


Because of rapid intersode relaxation and the occurrence of 
several high energy nodes in polyatosic solecules (Herzberg, 1945) oiv 
two node solecules can be taken as an approximate representation of 
real, mor.> complex polyatomics. To calculate the vibrational 
distribution we allow the following set of relaxation processes: (1) 
radiative excitation and de-excitation involving adjacent levels 
within a mode; (2) similar collisional processes and (3) rapid 
internal energy transfer between modes near resonance. Our initial 
treatment aasuaes that (3) is much sore rapid than either (1) or (2). 
We therefore obtain the steady state populations in mode 1 from 
equation (9) 


31 * = Nx+/\+l,x"'S+l.x^^x+1,x^ **x-1^^x-l,x"’Vl.x^ VS,x*3x.y^ 

*^x^*x,x-l'*^x,x-1*^x,x-l'‘’\,x-fl'^x,x-*-l'*^x,y*3x,y^ " 


where the suoscript x indicates vibrational levels of the low energy 
mode, and the subscript y indicates levels of the high energy mode. 
:>imilarly equation (10) was used to obtain the steady state 


populations in mode 2. 



~ **y^\,y-i'^y.y-l‘*^^y,y-i‘^y,y+i 


< Vi.yS-i.y> " “x'^^x.y^ax.y^ 


'*‘^y.y+i'^y.x*3x,y^ ® ^ 


( 10 ) 


In ttiis case, every level in the high energy node will be in resonance 
with a low energy mode. 


Figure 5 is a typical plot of the vibrational temperature of our 

hypothetical dual mode polyatomic molecule as a function of H atom 

pressure. The plot was generated in the same way as were those for 

the diatomic molecules rsing the following parameters; = 2000 

K, W = 0.1; mode 1, u s 1600 cm“\ TM r 10“^^; mode 2, u = 600 Qn”\ 

— 2C —1 

IM = 10 cm ; molecular weight of polyatomic molecule = 60. Each 
of the upper four curves represents the vibrational temperature of the 
low energy fundamental when the rate of intermode transfer (Process A 
in Figure 4) is set at 10”^, 10“^, 10”^ and 1 times the collision rate 
(C). Similarly each of the lower curves represents the temperature of 
the higher energy mode at these same intermode transfer rates. 


In contrast to our diatomic calculations, temperatures of various 
levels in a single mode nonnally differ by up to about 30 % from the 
average. As coupling between the modes is increased this discrepancy 
increases until in the highly coupled systems the lowest levels of the 
iow energy mode are at considerably higher temperature than are the 
upper levels. Tliese upper levels have temperatures comparable to 
those of the high energy mode and are therefore far from equiiibriun. 
It s.ould be noted that variation of the molecular parameters used to 
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Figure 3. Vibrational tmperature of 3 *’ 1 polyatoaic model as a 
function of total pressure and intermode coupling efficiency. Upper 
curves are the temperatures of the low-energy fundamental, lower 
curves are the temperatures of the high energy fundamental. 
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gtn«rat« Flgurt 5 not ohingt tht gtntral trends. Pwr in»t«no«, 
dtortaslng tht diitrgy or tranaltion noaont fbr olthor sod* still 
shifts ths ourvos to lowsr prsssurss. Siailsrly, inorsssing W would 
sst s loMtr liait to ths vibrstionsl tsapsrsturs ss in ths osss of CO 
or SiO» without sffsoting ths prssst«rs st which non-LTE offsets bsgin. 
Howsvsr, variation of ths collisions! sffioisney can no longer bs 
oaapsnsatsd for by a siapls shift of ths prsssurs seals sines this 
ohangs now offsets ths rslativs iaportanos of proosssss (2) and (3). 

It Wight bs argusd that our aodsl allows too fsw in traiol ocular 
transitions to glvs an aoourate representation of the upper levels ot 
the aolecule. However, the addition of aore intraRolecular 
transitions should tend to increase the rate of depopulation of the 
upper levels by allowing quicker access to the higher energy 
transitions wliich first depart from thermal equilibriua. In this 
respect our aodel appears biased towards thermal equilibriua, yet 
significant non-LTE effects are observed. 

In order to test the above hypothesis and so that we night have 
data on the behavior of systems where the ratio of the fundamentals is 
less than two to one, we constructed a model in which every third 
level of the lower energy mode was in resonance with every second 
level of the higher energy fundamental (Figure g) , We then calculated 
the vibrational temperatures for two cases. 

In the first case we allowed only resonant intermode transfer. 

In the second we allowed both resonant tranfer and non-resonant 


intcniod* transfer. The rata of noiwraaonant tranfar was oaloulatad 
using tha tbrauia of Hlllikan and Uhita (1963) and tiia approiriata 
anargy dlffaranoa batwaan iavala. In both oases wa obsarva 
significant departure fro« LTE as the pressure and radiative dilution 
coefficient decraasa. However, when off resonant transfer was allowed 
we find that vibrational temperatures in the transition region were 
significantly 01511) lower than when these transitions were forbidden. 
As the pressure continues to decrease however, the vibrational 
temperatures of these two cases converge. It should also be noted 
that the onset of non-LTE effects occurs at virtually the same 
pressure in both calculations, only the degree of disequilibriun in 
the transition region is effected. 

We now discuss our assumption that no molecule is simultaneously 
excited in both vibrational modes. Only about twenty percent of an 
equilibriun distribution of our 3:1 model polyatomic molecules 
(u»,| = lb 00 cro“\u 2 * 600 °'"""^) a temperature of 2000 K will be excited in 
both vibrational modes. As the vibrational temperature of tne 
molecule decreases with decreasing pressure and radiation density this 
fraction rapidly approaches zero. Furthermore, inclusion of all of 
the transitions which were excluded by this assumption would simply 
increase the coupling between vibrational modes. We have shown in 
Figure 5 that as the coupling between modes increases, the vibrational 
temperatures of the modes approach one another. Inclusion of these 
transitions would therefore not tend to keep any mode in thermal 
equilibriun. Since we are not attempting to calculate exact 
vibrational temperatures for any specific molecule, but instead are 



interested in discovering the general trend of polyatomic molecules in 
regions of low pressure and radiation density, we feel that this 
assuaption has little effect on our conclusions, he emphasize again 
that these model polyatomic molecules are used only as an aid to 
understand the effects of interaction between coupled modes in a large 
molecule. 

he are concerned primarily with the effects of vibrational 
disequilibrium on chemical equilibriun and grain nucleation. Thermal 

dissociation is a major process in such phenomena and occurs mainly 
from vibrationally excited molecules in states near the dissociation 
limit (Kondrat'ev, 1^4). Uniroolecular dissociation experiments 
(Rabinovich and Tardy, 1977) show that extremely rapid intramolecular 
energy transfer ('lo'^^s) occurs above the dissociation limit. In 
addition, high vibrational levels are more closely spaced than are 
lower ones. This tends to increase tne importance of near resonant 
transitions in int^rmode energy transfer, A third effect is that at 
high vibrational energies the molecule becomes appreciably distorted 
and a normal mode description is no longer applicable. Recent 
discussions of molecules in such highly excited states are given by 
McDonald (1979) and Rice (1975). 

(Xir treatment of the lowest vibrational levels of polyatomic 
molecules indicates that a considerable degree of disequil ibriun can 
be expected in regions of low pressure and radiation density. Because 
such disequil ibriun increases as the coupling between vibrational 
modes is increased and because such coupling increases as the modes 



becoae inor«asingly axoitad, on* should *xp*ot that tha vibrational 
nodes nost litcaly to b* out of aquilibrlua are those near the 
dissociation Halt. Since these nodes play a aajor part In the 
dissociation of polyatoalc molecules and since our analysis tends to 
underestlaate the laportance of vibrational dlsequlllbrlus In 
polyatoalc systems, polyatoalc molecules are expected to be more 
abundant In regions of low pressure and radiation density than would 
be expected on the basis of equilibrium kinetic or tneraodynaalc 
calculations. 


V. Discussion of Results by Region 

a) Atmospheres of Cool Stars, CircumstelXar Shells and Planetary 
Hebulae . In these regions, H atom collisions will be the dominant 
contributor to translation-vibration relaxation. Figures 2 and 3 are 
applicable in such regions. 

It can easily be seen Iran Figure 2 that CU is significantly out 
of equilibriua at pressures xess than 10 atm. Because we have 
deliberately overestimated the efficiency of collisional relaxation, 
this should be an underestimate of the disequil ibrlun in the system. 
From Figure j we see that biO begins to depart fran equilibrion at 
pressures less than 10~^ atm. If SiO-H atom collisions a. e more 
efficient than predicted by HilliKan and White (I 963 ) then these 
curves will overestimate the disequilibriua in the system, by using 
the proceedure previously discussed, if H atom-SiO collisions were 


found by Mptriatni to aotually b« 10*^ aor* effioimt than pradlotad 
by Hiillkan and Uhlta (1963) then the preaaurea in Figure i would be 
aultiplied by SiO vibrational populations would then depart 

froa equilibritn at preaaurea lesa than about 10 ata. 

Aa can be seen froa Figures 2 and 3 it is likely that CO and SiU 
will only be significantly in disequilibriua in the upperaost regions 
of cool stars. This is consistent with the conclusions of Carbon et 
al. (1976). According to Figure 5 and the previous discussion, this 

behavior will also apply to aost polyatomic species. However, in the 

lower pressure environment more characteristic of circunstellar shells 

-11 -12 

where pressures are expected to be in the range 10 -10 atm. 

(Hagen, 197b) both CU and SiO will significantly depart from thermal 
equilibrium. Similar departures should occur in the low pressure 
grain forming stage of an evolving planetary nebula where pressures 
are of order 10“^^ atm. (Oraine, 1979). these regions the color 

temperature of the stellar radiation field and the kinetic temperature 
of the gas should be comparable during the epoch of grain formation. 

b) Novae and Supernovae . Vibrational relaxation in expanding 
novae envelopes, and to a lesser degree in the outer envelopes of 
supernovae, will be dominated by H atom collisions. Therefore, 

Figures 2 and 3 also apply to these regions. In regions where H atoms 
are effectively absent, such as the innermost zones of supernovae, 
metallic atom collisions dominate. In an intermediate zone, 0 atom or 
CO molecule collisions could oe the dominant interaction. Figures 6, 

7 and b are applicable for Fe, 0 or CO rich environments respectively. 


Figure 6. CU vlDratlonal tcnperatures as a function of total pressure 


in an Fe-atoa gas. 
kinetic temperature 


Each line represents a calculation at the gas 
of the Intercept, ri s 0. 




Figure ?• CO vibrational tanperatura aa a fimction of total prassure 
in an O-atoa gaa. Each lina rapreaanta a calculation at the gas 
kinetic temperataure of the intercept, lil s 0. 
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Figure 8* CO vibrational taaparature as a function of total pressure 


in a CO-Bolecule gas. Each line represents a calculation at the gas 
kinetic temperature of the intercept. W s 0. 




PhotosptMrlo tMpcraturts for novM aro of ordar lO^K at tha 

tiAa of oOaarvad grain oondanaation (Sparks, Starrfiald and ITtiran, 

1978). Iha dust condansas at r 5x10^^ oi (Nay and Hatfiald, 1978) 

12 

Nharaas tha photosphara ooours at r <r 2x10 on (Sparks, Starrfiald 
and ITuran, 1978). Tha uitravioiat flux (i < 2000 it) is coaparabla to 
that at tha surfaoa of a SOOOK star and aiactronic transitions can ha 
ignorad. Tha enargy dansicy in tha infrarad (i > 1 nioron) on tha 
othar hand is only a faw paroant of tha SOOOK black body flux. 
Condansation of grains is ganaraily postulatad at prassuras naar 10 *^^ 
ata. Undar thasa conditions tha systoa should significantly depart 
00(1 aquilibriuB. 

Suparnovaa pnotospharic taaparatures at tha tina of condensation 

(t ^ 10^ sec) are approx iaa tel y 5-7x10^ (Kirshnar et al., 1973). At 

this time the outermost zone of ejecta will be at a radius 1 t >-100 

times larger than tha pliotosphare (Arnett and Falk, 1976). The 

innermost material (at r .r 15 r ^ ^ ) will experience a UV flux which 

phot . 

is comparable to that on the surface of a kOOOK star so that 
electronic transitions can again be ignored. As in the case of novae, 
the infrared flux will be only a few percent I that of a 400QK black 
body. Pressures will be below 10**^^ atm. Condensation would 
therefore occur in a system which is again out of thermal equiiibriun. 

c) Primitive ^Solar Nebula . Figures 2 and j upply to models 01 
the primitive solar nebula, because of the enhanced efficiency of 
non-adiabadic H atom collisions it appears probable that many, if not 
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All, aoltouiar ap«oi«s will b« in ttwraal tquUlbrlui undtr th« 
oondltlont pr*dlott<t in aost ao<t«la of tht nonula (at* for Inataaoa 
CMtron and Flna, 1973; Caaaron, 197B)* Only tht lowaat prtaaurt 
rtdlona (p ^ 10*^ ata) will ahow non-tquilihriua population 
diatributiona. 

VI. Diaouaaion of Htaulta; Koltcular Diatributiona 

a) Holtoular Paraatttra . In syattaa oonbaining niatroua 
aoltoular sptcita, aany ooabinations of aoltoular paraatttra art 
posaiblt. a aoltoult with low Einattin A and rtlativtly hidh 
which would bt tht oharaottriatios of a aptoita with aaali dlpoit 
aoatnt and a low fundaatntai vibration frtqutnoy, would rtaain oloat 
to tqullibrlua tvtn at rtlativtly low prtaaurta. Qonveratly, one with 
high A and rtlativtly low - a aptoita with a high dipole aontnt 
and nort tntrgetic vibration fundMtntai - would dtpart froa an 
ttquilibrluB vibrational diatribution at higher prtaaurta. Thtrtl'ort, 
difftrtnt aoltcular aptoita in a particular region night bt txptottd 
to be out of equilibriuB to different d^rtta. 

For diatonic noltculea it ia poaaibit to eatinate the vibrational 
ttnptrature of a particular aoltoular aptoita if one knowa both tht 
tranaition nontnt (TH) and vibrational fundantntal (w) by the use of 
Figurta 9 and 10. Kigirt 9 ia a plot of the vibrational ttnptrature 
of a diatonic moltoult aa a function of H aton prtaaurt in the abatnct 
of radiation at a gaa Kinetic ttnptrature of 20UUK for linta of 
oonatant tranaition nontnt. Tht noltoular weight and vibrational 
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Figure 9. VlbretioneX teapereture as a function of H atca pressure at 
a gas kinetic tenperature of 2000K in the absence of radiation for 
lines of constant transition ■onent. All other aoleoular parweters 


2000 




\a' 
(Ct V 




ALL PARAMETERS BUT 
TM ARE THOSE OF CO 


Z 


fundancntal art ttoat of CO. Fifura 10 is a ala liar plot for Unas of 
oonatam: u. Again tha aolaouXar waight and transition aoaant ara 
thoaa of CO. For both figuras tha C. . oaloulatad froa tha work of 

1* W 

Hillikan and khita (1963) has baan aultiplied by 10^ and should 
tharafora undarastiaata tha vibrational disaquilibriiai of tha aystaa. 
lha use of the aolecular weight of CO in the calculation introduces no 
error into our estiaate since only the reduced aass of tha colliding 
systw is of iaportance. Since we already assuae that all collisions 
are with atcaic hydrogen this tera is very nearly unity for aost 
aolecules of interest. 

The following procedure will yield an approxiaate vibrational 

teaperature for a diatonic aolecule of Icnown u and Ifl. First find the 

line of constant transition moment corresponding to the aolecule of 

interest on Figure 9 and note the difference (in units of pressure) 

—1 y 

between this line and that of CO (1 x 10 ). Next, add or subtract 

this difference from the pressure scale of Figure 10. Now find the 
line of constant w which corresponds to the vibrational fundamental of 
the aolecule of interest and read the vibrational temperature as a 
function of H atom pressure directly from the previously modified 
pressure scale. 

b) Non-Adiabadic Collisions . Kelaxation of molecules via 

non-ad iabadic electronic or "complex forming" reactive interactions 

(Saith, 1976) is again expected to depend specifically on the dominant 

colliding species. These species can vary in efficiency by factors of 
14 

10 (Glanzer and Troe, 1975) for the relaxation of the same molecule. 
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Flgur* to. Vibrational tMperatwa as a function of H aton prtsaura 
at a gaa Kinatlo toaparatura of 2000K in tha abaanoa of radiation for 
iinaa of constant anargy of tha fundiaantal vibrational moda. All 
othar nolacular paramatars ara thoaa of CO. 



Although it sons liicoly thst stoMio hydrogsn will bt ths doninsnt 
oolliding spsoiss in aost systems of astrophysiosl intsrost, fow 
aessisrncnts hsv* boon aodo of its rolnstion offioionoy for nolooulos 
oxpoctod to occur in those regions. Suoh aoosuroaonts are necessary 
before aore exact calculations of the vibrational populations of these 
aoieoules is possible. 

c) Assumption of Theraal Eguilibriua . Holecular dissociation 
preferably occurs from the upper vibrational levels (Snith, 1976, You 
and Pritchard, 1979) • If these levels are depleted for any reason, 
then the molecule is aore stable with respect to its atomic 
constituents. This is amply demonstrated by the calculations of 
Dalgarno and Hoberge (1979), cited earlier, tOtich found that the 
molecular dissociation rate for CO and, to a lesser degree, for H2, 
decreased by orders of magnitude as the total density decreased. 
Holecular equilibrium calculations, such as those by Tsuji (1964, 
1973), can therefore not be applied to regions of lower pressure. It 
seems clear that new calculations which account for the varying 
degrees of vibrational disequllibriim in regions such as circuastellar 
shells, novae and planetary nebulae are called for. idefore such 
molecular distributioh calculations can be attempted however, one 
needs Einstein A values for groixid electronic state vibration-rotation 
transitions and H atom-molecule relaxation efficiencies for each 
molecule of importance in the region under consideration. In 
addition, a suitable model of the rate of internal energy transfer for 
polyatomic molecules is required. Ideally, state selected reaction 
rates for individual vibrational levels should be used in such regions 
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although, dut to tha oohplaxlty of tha aystaas, aoaa suitahla 
approx iaat Iona alght naad to ha aada. Unfortunataly, alaoat nona of 
thaaa data ara ourrantly availabla In tha lltaratira. 

Tha inplloationa of thaaa raaulta for grain eondanaation ara 
conaidarabla. Tha non-aquUlbriua eharaotar of nuclaation for ooMio 
aystans has bean discussed in several papers (Donn, 1975. 1979a) and 
experiaental evidence applicable to astronoaiioai systaas has been 
obtained (Day and Donn, 1978a, b) . The analysis of the present ' 
chapter shows that the gas can deviate significantly frj a^quil ibri lai 
before condensation begins. Thus, it is wilil^ktf^That classical 
nucleation theory is applicable. A procedure for treating 
condensation on a purely kinetic basis has been proposed by Donn 
(1979b) and is being developed further (Donn et al, 1980). 

It is interesting to note that Arrtienius and De (1973) and De 
(197d) have determined that snail grains embeddr^a in a low density gas 
may be much cooler than their surroundings. The present work has 
indicated that many molecular species are similarly out of equilibrium 
at low pressures. The evidence is now compelling ttiat thermal 
equilibriua in low density clouds is no longer an acceptable 
approximation. The extent to which such an assuaption may be a useful 
approximation must be examined for the specific problem under 
discussion. 
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VII. Concluaiont 


Vibrational disaquilibriui baooaaa inoraaaingly signifioant in CO 
and SiO as tba total density of partioles and radiation deoreaaes. 

Sooie degree of diaequilibriua. uhioti will depend on specific •oleouiar 
properties will occur for aost heteronuolear aoleoules— ^including 
polyatoaic species— in low density clouds such as the expanding shells 
of cool stars, novae and supernovae. Previous work which assisted ttiat 
themodynaaic equilibriua had been attained in such regions must be 
reexasined in light of the iresent calculations. Einstein A values 
and H atoa relaxation efficiences for di and tri-atoaic netal 
hydrides, carbides, nitrides, oxides and sulfides are necessary before 
more accurate calculations of nolecular abundances under these 
conditions can be made. Vibrational disequilibriin has signifioant 
implications for the calculation of many chemical processes and 
characteristics of low pressure, low optical depth r^ions. 
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Cnaptftr 4 


Extinction Iteasureaents of Iran and Magnetite Grains 


Both l! on (Fe) and Magnetite (Fe^O^) have previously been proposed as 
constituents of interstellar grains (Huflman* 1977) • Fe grains have 
been pr'oposed on the basis of coaoic abundance (bchalen,l9j6) , as a 
universal condensate (Lewis and Ney, 1979) • and as a possible cause of 
ttie 4.6 \t~^ mo) featu'^e in the Interstt^iiar Extinction curve 

(Hufllaan, 1977). Fe-^O^^ has been considered because it nas been found 

in oieteorites and because it is a possible cause of the Very broad 
structure in the Interstellar Extinction curve (Huffman, 1977; 

Uiiittet, 19b1). If eitner type of magnetic g^ain exists it would help 
to explain the polarization of starlight which is tnought to be caused 
by g'‘ains aligned by the interstellar magnetic field (Huffman, 1^.'7). 

This chapter reports laboratory measurements of the optical 
prope'*ties of Fe and Fe^O^ grains. Tnese results are discussed in 
terms of the effect of suet particles on the extinction observed in 
the interstellar (nedi>»ra o'* in circuti stellar shells. Special regard is 
given to tne recent reports of anomalous extinction curves which lack 
Che 4.0 featu»'e. This chapter represents a coiiaooravive effoi"!. 
between I>r. dixaes Hecht (NASA/NhC Kesident hesea**ch Associate at 
Uodda'*d Epace Flight tenter) and myself. Dr. Hecht performed 



theorttlcal Hie oaicuiatiaes on amaii and magnetite particles. 
Hy prima'-y responsibility was to obtain tne experimental extinction 
and size distribution measurements. 


II. txperiwental Procedure 

Tne tlie’*faal dissociation oi iron pentacarbonyl , via the reaction 

Pe(CQ),- Fe ■«- bCU, occurs at temperatures in excess ot' idOO^C 
0 

(heilor, 1962). Ttu reaction proauces an extremely supersaturated 
iron vapor wnicn quicKly nucleates (Frurip and bauer, 1977) yielding 
solid iron particles. In ou’* experiments ve«’y small iron grains 
(<r25nm in radius) were produced by tne thermal decomposition of dilute 
mixtures of Fe(tU)^ in either an argon or hydrogen oar»*ier gas. 
Dilution was sucn that 1:10 > Fe(CO),: Car'-ie** Gas > 1:1000. Total 

pressures were from 10-tiO torr. Dissociation occured at temperatures 
in excess of 27b^t and was accomplished as follows. 

Carrier gas, the flow rate of which was regulated by u needle 
valve, was allowed to flow into a flask containing Fe(CO)^. Tlie flask 
was immersed in a slush oath to regulate tne partial pressure oi 
Fe(CO)^, The use of various slush baths yielded a range of Fe(CO)i^ 

pressures between U.0‘j torr (chloroform slustO and 2.6 torr (benzyl 

Alcohol slush) in discrete steps. The gas mixture was then allowed to 

flow into the bottom ol a furnace which was enclosed within u 160 1 
stainless steel bell jar. Constant pressure was maintained within the 


systea by pumping at ttie top of the furnace using an auxiliary 
mechaniudi pump, before each run the bell jar was puaped for at least 
two hours, and in most cases overnight, to pressures less than o X 
lu”^ torr at temperatures in excess of jOU^C. 

Itie furnace is a resistively heated mullite cylinder 6 inches in 
uiametei* and 1b inches high through which four equidistant, 
equato*'ial, one inch diameter, holes had been drilled. Six incn long, 
one inch ID alumina tubes had been inserted into each noie and 
aligned with the four window ports of the bell jar so that the 
interior of the furnace could be viewed. These tubes eliminated 
significant perturbation of toe temperature profile of the system. A 
iOOW Xe arc lamp was placed in front of one of the HgF^ windows and 
qualitative obsevations of the ^oke color and intensity could be 
obtained by looKing through tne window port wiuoh was 90^ to the lamp 
pO'^t. by using a polarizing filter an estimate could also be made of 
th;. polarization of tne light scatte'^ed by tne soioke tnrough an angle 
of to the incident beam. Extinction measurements were made by 
placing d ioonoutiromator by tne port opposite to the lamp port. 

Heasurements were made with a O.dm scanning monochromato'' in the 
range fr>o between 1.4-b. between l^U-ubdrin). A steady rate 

of pa'^ticie production was established within the system and an 
extinction spectrun I , (i) was recorded. Tne :iyst«u was then 
mechanicaiiy puaped to remove ail particles from the gas phase and an 
unobscu'*ed ioiap spectrun, I (X) was taken. The quantity in (I_/i 
IS pro po port! onal to the extinction of the anoke. because of the 
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iMrge rang* In the intensity sf the Xe laap, a*asur«i*nts were aade in 
eight spectral ranges and scaled to sne another in regions cf mutual 
overlap. The results, snown in Fig. 1, will he discussed in a later 
section. A broadband filter (FWHHcMOnm) centered at 4.6w''^(216na) was 
used to obtain the data for < 23bnm). This filter, 

placed before tne entrance slit, eliiainated background caused by 
scattering of visible light within the monochromator . 

Pa'^ticles collected during our runs and analyzed by x-ray 
diffraction within several hours of their production were found to 
consist of ot-Fe, and y-Fe. We believe that the Fe^O^ eicne»* 

formed in transfer to the diffraction apparatus or from iow level 0 ,^ 
contamination which *”eacted with Fe particles on the wails of the 
system, inis could have occurred during the two hour period while the 
system was cooling to room temperature. Granqvist and Buarman (1^76) 
also reported finding Fe^O^ in their Fe particles even under extremely 
clean conditions. Ue do not feel that such contacaination could nave 
affected our measurements of Fe particle extinction for the following 
reason. 


Our wens'jrements involved only tnose grains suspended in the ^as 
phase. Cbse''Vdtions indicated that the residence time of individual 
particles, in tne gas phase , was on the orde»* of minutes even at the 
nighest systen p»*essure used. As tlie particles settled or were pu.iped 
away they we'‘e replcnisned by freshly nucleated grains. The leak rate 
neces3a'*y to oxidize a raeasurabie fraction of these grains on this 
short time scale would De compa'^abie to the rate at which FeCCU)^ was 
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Mj^ure 1 , LxpttriaientMl data points r«pres«nt maasuraiaent of 
extinction vs from 1.5 to b.O in eight ove*‘lapping spectral 
regions (see text). Ihe light scattered at 90^ was lUOt polarized in 
ail spectra, Ove**lapping points are scaled to one an other. Th.i 
absolute scaling of the extinction is arbitrary. Hie solid (dashed) 
line represents calculations using JC(H(;) data, as desc'ibed in the 
text, of L, the extinction per unit length. Fo** tne JC calculation 
the point at is arbitra''ily fit to the data. Fo" hO the point 

at 5.0 is arbitrarily fit. 
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introduced into tne furnace, a leak rate would have made it 
impossible to p>jtnp the system below 10 mic'‘ons mucn less beldw 10 
tsrr. 


Particles allowed to accumulate on tne walls of tne furnace for 

0 

two weeks at temperatures greater than oOO C were fomd to oe pure 
Fe^U^. Analysis by a scanninii electron microscope (5m) indicated 
tnat a typical particle had a radius of 2bna altnomn b-10'ib sf the 
total nad ’^adil as iur^e as ^Una. A con3iae<^aole amount of cluopin(i 
could be seen in the 5111 pnotographs although we feel, because of 
polarization measurements discussed later, that most, if not ail, of 
this occurred on the walls rather than in the gas pnase. 

The extinction spect'-a of the g**ains collected froio ou»* 

syst«ai) was obtained uy fi'*st dispersing them in ethanol and tnen 

despositing tnis on a clean sapphire window. The ethanol evdpo*'ated 

and I . was recorded. Tne window was then cleaned and I was taken, 
ext 0 

Ln is shown in Fig. 2 and discussed later, 

Hnally, measurements of In d./I we*^e made for large ( > bU 
nm radius) Fe particles using a low dispe’^sion (14na/i.ui) monoctroiaator 
equipped witn an optical multichannel analyzer. This allowed the 
acquisition of spectra, in the '^ange of (bjOnu) to 2.b!iu”^ 

(oYbnin), to be coiapleteu in seconds. It tnus became possible to 
obtain I over the entii-e specfal range iii a time small oompa'*ed to 

the rate at whicii the size dist'^ibution of tne cloud was changing. An 
exanpie of tnis is shown in Fig. j. Fig was taken just after a 
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Figure i». i>oiid line is a plot of b Vb x”^ arbitrd»-ily scaled to 100, 
for a "distribution'’ of Fe^O^ particles wne**e N(25na) is 13 tiues 
N(bOniD). Points represent experimental measureiaents of tlie extinction 
of Fe^O^ particles, collected from ou»* system. Tne experimental point 
at t>.0p”^ is arbitrarily scaled to the tneo>*etical curve. 
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figure 3. I^ie S3ila lines are plots of iceasureu extinction uiinus a 


constant (Y^) vs taken wittt an optical uuiiticnannel andlyze**. The 
solid points are a theoretical calculation of L, for a best fit to 
this data, usin^ a distribution of particle sizes (see text) of radius 
6b to 100 mi in 3 mi bins. 

A. Taken after a steady state cloud was f onaed . NC fo'* a < Yb nm 
IS four tiises HC for a > 7b mi. 

b. Taken after cloud n^s settled for 100 seconds. Nc fo” a < 

7b run is ten times NC for a > 7b r«n. 
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steddy state particle distrlbuticn had been established. 'Ihe t'laSK 


csiitainiiiii Fe(CO)^ was then immersed in liquid which eliminated 
particle prsductisn. Fig jb snows In taken 10U seconds 

later aite*' some oi' the la**ge’* partiuies had settled out. bkM 
pictures ot' ttie snoke settled on tne furnace walls showed particles 
t'roin ‘jO to lOUUnu in radius. 


111. Co.nputdtional Procedure 

I'rit FULLUi^lNU IS h SUrtMAHY JF irit CUIPUTiillONAL PHOCLDUHk, UbtU IN THlo 
STUDY hND was WKITILN dY Dh. HSCHT, WHO ALSO PLKFOhMLO THL 
CALCULAnoUS. 

The use oT tile theory to calculate tiie extinction. In (lo/I), nas 
been discussed by many authors (see t'O'* example Itecht, 19^9 or 
Hut'lman, 1^77). uriet'ly an ei't'iciency factor, g, can Ue calculated 
tor a given particle ■’adius, a, anu incident pnoton wavelength, X. 
b(a) the extitiction pc'* unit length, is then equal to w times the 
pa**ticle nusibe** density, N, ana the geometric cross-section, C. Ihe 
total extinction pe** unit length, £, To'* a given wavelength is equal 
to tne sua oT L(a) taken ove'‘ all pa'‘ticle sizes. Itiis can be 
c>upa>’ed with expe’^lment uy no'*maiizing the peaK ol eacn cu^ve ( tneo''y 
anu experLneiit) to a coiikaon value ana onparing tne shape of tne 
resultant cu'*ve as a function o: wavelength. To make these 
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oaicuiatisns, the caaplex diaiectric constant, c > ♦ i c^* 

needed. For Fe these were taken t'ron Johnson and Christy (iy/4, 
hereafter JC) and from Horavec, hife, and Dexter (1976; nereefter 
HU); for these were taken from SclUe^ei, Aiverado end Wachter 

(1979)). 

Caiuuidtions of Ha) were made from i;>M to 6u;fnin(d.1 to 1.7 m~^) 
in b nm steps, for pa*‘ticie sizes from 9 nm to lOOnm, in rudius, in 
5nui intervals. Hie tneory also uiiows tne polarization of tne light 
scattered at to be calculated and this was done ove*” the same 
ranges. 

In S0(ae calculations it was found necessary to calculate L(a) fo” 
a coated gram sucn as 1'e^U^ coated with H^U ice. Calculations were 
made fr>n loU to na(Y.69 to ;j.4o w**^) in '^ma steps. In these 
calculations a mouified Hie theory is usee (Kerker 1969; hecht 1979). 
The c (1) for ice we^e taken f'*Oio Jreenberii (lyOb) and Dressier 
ana bchnepp(196U) . Linear inte*’polations were made uetween points. 

riathij, ham pi, anu horasieck li977; he^'eafter Mhl«) have shown 
that tne Interstellar kxtlnction curve enn ue fit wtien N(a) da, the 
nunber density of particles between a and a^da, is given by 

N(a) da = ka da (1) 

wner,; K IS a constant depending on tne pa'*ticle material, oince in 
tills paper only relative extinction curves are calculated, K. can be 

loo 
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set equal to 1. For Fe and Fe^O^^ ue have calculated E tor a MRN type 
distribution of particle sizes fron 2.3 to 102.3 na. Tnis was done by 
taking 20 bins 3nui wide starting at 3nn and ending at lUOna. Hie 0 
for each bin is constant (i.e., U(3nn) . W(IOnra) ...gdOOiia) and E(a) 
for each bin is integrated using eq. 1. The 20 values of E( a) are 
then added to give a MhN type extinction curve for Fe and 


IV. hesults 

I»*on 

Calculations ; Using values fo*” c front eitner JC or HO it was found 
that t(a) could peak anywhere from 7.Yw”^ (IJOnin) to (bdOiin) 

depending on pa’'ticle siz«;. Host of tnese peaks a'*e extremely broad 
Lhall wiatti at half max imun (riWHH) > lUOrinJ. Morever, tiie position of 
tne calculated peak is extremely dependent on particle size. Using 
tfia constants of JC(MO) a 23 nm particle had a peaK at 6. (3.Cp~^ ) 
and a KwhM of 144ria (106 nm) . For comparison a ^Unu radius particle 
had a peak at 3.0p”^ (3. 1p”^ ) ano a HWhM of 136 nm (126 rm) . Tne HHk 
results, plotted in Fig. 4, snow that using the data from JC, t 
increases witn energy, wniie tne data from Ho predicts a peak nea<* 

6. Uu \ 

Measurements ; Pa*'ticle sizes determined from 6tH pictures indicated 
that individual particles were nea" 23 nm in radius with a few larger 
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Fiiiure 4. Caicuiatisns of k vs for a MRK type distribution ol' >■'€ 
parcicles (:>ee text) . 

A, Taken I'rou Koravec et al.(HU) data for tne dielectric 
constants. 

b. Taken fro.a Johnson and Christy(JC) oata for the dielecfic 
constants . 

The peak tfal’je of L is aroitrarily normalized oo 1UJ. 
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particles present (50niu). Fsmatisn sf larger particles, via 
ciuaping, prsbably did not occur in tne gas pnase. 

Visual oeasurenents of the optical polarization indicated tnat 
tne light scattered at VO was 100( ve’*tically polarized. Froa 
previous wo'‘k this indicates that the ratio of vertically to 
horizontally polarized light was better than 100:1 (Hecht, ivdO) . Hie 
calculations indicate that this is consistant with a gas phase 
particle size of less than 40 nm. 

Fig. 1 snows the results of the laeasureioents from l.b to 4.V 
Also shown are some theo’‘eticai calculations for a two bin 
distribution where N(2bnn) is 1j times N(bOno) . As can be seen the hO 
points lie within the experimental uncertainty while the JC points 
give a slightly poo^'e'* fit. 

Fig. j snows the '■esuits lor tlie UMA eperiment where tne dust 
cloud was allowed to settle. I40tice that the plateau wtiich begins at 
2.4 14 **^ has disappeared after 100 seconds. The model used to fit 
these curves results in relatively fewer- large*- particles being 
present after 100 seconds. This is consistent with larger pa*’ticles 
settling out faster. Ihe plateau is not a small particle resonance 


featu'*e 


Magnetite 


Calcuiatians ; Ine E(a) for can peak over a wide ’*ange from 

.6u~^(17bnm) to 2.0 w~'\500run). The variability in E(a) witn 
particle size is much less severe than with Fe. A 2bm particle has a 
peak at 196 nm and HWHH of 114m, while a 30 m particle peaks at 202 
nm and has a HWHH of 120 m. This suggests that a real absorption 
occurs and it snould be noted that the measured absorption coefficient 
of magnetite does show a peak near 200 nro. 

A liKN type curve was also calculated and is shown in Fig. 5a. 
Noteworthy is the shoulder starting nea*" 5u”' (200m). Also shown is 
a modified HKN curve. Fig 5b, where an additional nunbe** of 40 nu 
particles, equal to .007 the density of 5na pa<'ticles in our standard 
HhN distribution, are added to the dist*'ibution shown in Fig. ja. 

Note that tne snoulder shifts to longer wavelengths. 

Measurements ; bince the Fe^U^ particles were presujabiy Fe pa'-ticles, 
which wee oxidized after having initially settled on the furnace 
walls, the size distribution sr.ould be similar to that of tne Fe 
experiments althoi^h climping could add a signiticant nunbe'* of la^ge** 
particles. Shown in Fig. 2 is the Fe^O^^ extinction data plotted on 
the same g'*apn as a dual uin distrioution of Fe^t)^ particles wne^'e 
N(25 nm) is 13 times N(50 nti) . 
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Figure 5, Calculation of £ VS x“^ f3»* a HHN type distrioution of 


Fe^U^ particles (See text) 


A. Standard HKN distribution 


b. Same as A with some additional 40 m radius particles (See 
text) . 


£ is a'^bitrarily normalized, at its peak value, to 100. 
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V. Discussion 


Tne results ano procedures of the previous sections will be 
discussed in relation to (A) the 4.6 u~\220 nm) peak in tne 
Interstella** Extinction curve, (b) the Very broad Structure in the 
Interstellar Extinction curve, and (C) the recently announced (Sitko, 
Savage, and Heade, I9b1, hereafter SSH; Snow and Seab, l^oO) anomalous 
extinction curves which snow no prominent nm) feature. 

(A)- Huffman (1977) in his review article on interstellar grains 
suggests that Fe might be responsible for the 4.6a~^(220 nm) feature 
in the Interstelid'" Extinction curve. It is known that, for small 
spherical parti'; .es, when is equal to -2 and is near 0, a 
resonance occa'*s in the small particle (less than 50 nm) extinction 
cu''ve. because the MO ana JC data for e.| are in disagreenent Huffman 
stated tnat if their measu-ed values of which were in agreement, 
were too high oy a factor of 2 then tnis resonance might be expected 
to Occur anywhere between 5.0ii~\200 nin) and 1.bbn~\b00 nm) depending 
on tne f'ue value of c^. Ihe '*eason Huftbian gave for the possible 
error in the measurement of was oxide contamination. 

In our exreriiaents we saw no evidence of any resonance peak 
betwieon 4.7bu ^(210nm) anJ l.oou \bUUnia) even when small (less than 
50 nm) particles were present, in fact our data is consistent with 
extinction cu'^ves de'*ived from the hO data and possibly even from the 
JC data. This is because, except wne**e is near u, the extinction 
IS sf'ongly dependent on (tne value ol whicn both g'^oups agree 


upon). Therct'ore, .lie Fe cauld certainly contribute to the 
Intersteilfir Extinction it does not seem to be responsible for the 
4.b|i~\220 niu) feature. 

(b) While ou** experinentdi data did not try to •^esoive the Ve»*y 
h'oad btruccure, the th*>oreticai calculations of E which we made 
Showed some interesting consequences of tiie hypothesis that the Very 
broad Structure is caused by magnetite, (Huffman, 1977; Van breda and 
Wnittet, 19d1). If the approximate 0.02b mag/kpc VBS is caused by 
magnetite than tne magnetite contribution to the Inte'^steilar 
Extinction cu''ve at b.Op~^ (200 nni) would be ahout one fifth tne total 
extinction. This would make magnetite an important component of tne 
mid UV portion of tne Interstellar Extinction especially in •'egions 
where grapnite may not oe p'*esent. However, magnetite is known to 
hove two strong IH aosorption bands nea*” 10 and 2oani. These features 
have been tentatively identified in Cl carbonaceous chondrites which 
are known to contain magnetite (Huffman 1977). If magnetite is a 
major component of interstellar dust grains then these featu''es should 
be seen as either IH aoso’-ption or emission featu''es nea** IR sources. 

(C) Recently SdH and bnow and beau (19b0) have reported 
extinction measurements whicn do not show a 4.0 u“'(220nm) fedtj'"e. 
Hiese specf’a do seem to show a snouide** wtiich begins at somewhere 
oetweeii 4p“U2oO rin) and 5p~U200nm) . Fo'" discussion purposes we will 
divide these spectra into three catagories: (1) those which show a 

Change in tne slope at about o . bg~^ ( ) VUtia) , (2) tiiose which stay 
nearly flat past 170nu) , and the specfun of HD441?y. 
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The extinctlsn sp«ct’‘a of HD163296. HOi164b, HD190073* AB AUR, 
were discussed in BSM and all seem to shsw a tooderate silicate IR 
feature at 9.7i<a. SSH nste that recent ws<^ic l)y h>anpl indicates that a 
rise in extinction at x”' > 5.ba”\x < 170 ria> is typLcal of 

* silicates. We also note that recent calcui-ations by Mathis and 
Wallenhorst (19dl) confinn tnis increase, especially when tne size 

* distributions include aoall (lOtxn radius) silicate particles. 

Tnerefo”e the increase in extinction past b.bp”' ( 170no) is probably 
due to silicate particles. Althoi^h no IR spectra are available for 
HO ^9o47(Snow and bead, 19d0) tne silicate explanation prsbably also 
accounts fo** the increase in its extinction cu'*ve. 

These spectra also show a weak UV o>jap and/ o'* shoulder oe^mning 
near *».i! y”’ (240nrn) in HD d164b, HD I6j296, and HD 190073, and near 
3.0 n”'(200 nin) in Ab AUH and HD 29o47. bucn features are 
qualitatively similar to the magnetite featu'-es shown in Fig. 4. To 
make a quantitative fit tne exact particle size dist''ibutioti is 
needed. It woulu seun to oe more pi^ofitaDle to attempt obse'‘vacion of 
tne IK leatjres of magnetite. If suen features we^'e seen tlien there 
would oe strong evidence fo'" the presence of Ke^O^^. 

r 

The extinction spectra of HD 2o9^31» oD 43o7Y and HD bOljd s«iow 
no increase in UV extinction and no silicate II: feature. Fresunably 
no silicates a»'e present a»' 0 'jnd these stars. Tne shouldc" oOse*’ved in 
HD 2o94j1 and HD bOIdd are qualitatively simila** to magnetite 


Wo 


’ extinction. As aoove, a search for laagnetite Ih feature: would 

provide mo'*e definitive proof ot the p**esence of oaiinetite. 

t 

I Ue do not believe the sui^estion that Fe grains have 

{ condensed around these stars acco'*din^ to the scneme of Levrt.s and Mey 

(l9Yb). This schaae was based on the assumption of equilibriua 
condensation yet it is known f'-om laboratory studies (Fru»*ip and 
bauer, 19YY anu bauer and F»'urip, 19YY) that Fe condensation is a 
sf'on^iy non-equil ibri jm reaction. Also, at low densities around 
stars, the vio*"ational disequilib<”iun problems discussed by Nuth and 
i)onn (lydl) can neone important. Tne»*efore wlule pu'*e Fe 
condensation is not impossible, it is not as likely as Lewis and Ney 
p''oposea (Donn et al . 1b>oD. 

Finally, it snould be noted that hi) 43&YY is p’'Obabiy not a good 
candidate fo” magnetite dust. Itiat spucf'uii snows a sha^p inc*ease 
f-'orn 4 to 4.5n“^ (iJ50 to ktiO nin) with flat extinction before and after 
those energies. This •'esembles what would be expected from the saall 
MgU pa*'ticies proposed by Dulcy, dilia''. and kiliiaias (lyYy). Tnere 
IS even an IK feature near ^0 pn whicii would be expected from suen 
particles . 

The extinction spectrun of hJ441Yy, the iliunination source for 
tne hed hectangle Nebula, is unique. It shows a snoulde** beginning at 
4. ou~^ wriicn mignt be due to magnetite, howeve'', it also snows 

a large inc’ease in extinction wnicn begins near b.ou~^n^d nm) and 
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continues to aDout 6.Bp~^(147no). SBH suggest this could be due to tne 
P'^esence of eitner dust, tnoiecula'*, or atomic species. 

In Fig. 6 we snow a calculation for tne b(a) of a 15 nm 
particle coated with a 7.5 nm H2O ice sneli. Also snown are selected 
^ points of the HD 44179 spectrun taken from SBM and scaled to ttiis 

curve. Inis fiisure snould bo viewed with tne Knowledge tnat tnere are 
• several variables, i.e. the exact values for c of ice between 6 

and 7 p~\l67 and 14]j nm) , tne size distribution of magnetite 
particles, and the thickness of tne ice sneli. hevertneless the 
feature is nearly reproduced. (Note that a sepa<*ate Jist'*ibution of 
uncoated magnetite particles and il2U ice particles could also be found 
so tnat tne feature is reproouced) . 

IH spectra of tnis object a'-e available (Hussell, Soifer, and 
kiline*' lV7o) and no 3.I itn featu*-e is seen, however tne IK 
spectrun appears to ^rise from an optically trim sneli ana in oroe” 
for the coatings to exist they must be in a cool region far away 
frotn tne excitation sou'*ce. From tile strength of tne UV extinction, 
tne j. 1 pin extinction should be about '.UW of the incident light. 
Furutermo''e , since tno UV and IH radiation appea'’ to co.ne fr>n 
^ different ■‘egions ( Conen et ul . , 1975) the ice coated uust may be a 

separate g’-ain population fr oiu tnat '•esponsible for the inf»*a'"ed 
emission. Ihe j. 1 po featu>*e would tneti not be p’*esent as a dip in the 
emission specfun. ligain tne identification ol magnetite would be 
confinned if tne IK magnetite oands could be detected. 
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Kiii""e 6. Solid line is a plot of E(a) vs x”^ for a 1b nio Ke^O^ 
particle which has a 7.5 nn thick H^U ice shell. Tne crosses are 
experimental poihts taken from Sitko et al.(SSH). The e'*ror bars 
represent ou^ estimate of the uhuertainty of selected points. The two 
points at 3.3 and 4.U a^'e arbitrarily scaled to the plot. 
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EXTINCTION (E ) 




VI. ^imiciry 


i‘iiL3 scuJy nrfS aeni3n3Criit.ed cnat soiaii iron p<i'‘clcies arc rasst 
likely n3l responsiole I'^r tne (?20nm) feature 3f Uie 

Interstellar Lxtiiictisn cu'*ve. Ou'* uate is cjnsisteiit wltti Pstn 31 
tne previous laeasureineiits 31' tne aielectric constants 3i iron, 
dltnomn tne values ri M3'‘avec, Rife anij pexter (1^70) yield a 
Silently betrer tic than Jo those of Johnson anu Christy (197^). 
HeaS'j''eiaents of tne '*elativc extinction efficiency oi small magnetite 
grains are quite consistent with the ”.cent detennination of tne 
optical constants of tnin magnetite films oy 6cnleg«i, Alvarado and 
Uacnter (1979) • Calculations indicate tnat if magnetite is indeed 
•"esponsiLiie fo>" tne Very droaJ bt''uctu'*e of tne Interstellar 
Lxtinction Curve, tiien it snoulu also account for about of tne 
total obSs;'’ved extinction at gUJun. If ^''ap-iite we*’e ausetit Uien tne 
cont'*ibution of magnetite to tne line of sigtit extinction could be 
muoti ino''e significant, ae hcvc snown ttiat r- is could actually ae an 
explanation for tne ananolous extir.cr: on curves •■ecenlly reported 
towa'-d seve''<ii stars (diCKO, aaVcge and hiCade, I 9 p 1 ; Ciiow anu deab , 
19o0). 

We nave deni onsf'ated that tne ex pe''u.ientai appa''atus aesc'iCed In 
ttJt p'-evious ;.-ection is suitable l or tne production of ve'^y small iron 
pj'‘ticles. irie s*ze dl3t'‘iuution and concentration 01 tness pa»"ticle3 
Call be controlled quite easily for ••elatively long b«?'"iOds of time so 
that iacasu'’e(iients of tne optical extinction of tne Sinoke Can ue made. 
lnco'’po'*jtion ol an oHh into tne systeo allows data aquisition on a 
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aucn shorter tiaescale. If this were coabxned with a real time 
particle collectisn system it shsuld be psssibie ts sbtain nigh 
ressiutisn extinction spectra of very well characterized particle 
distributions. Sucn a system may become available at GSFC within a 
year. 

It should also be noted that uhe flow system used to produce ! e 
iron particles can easily be adapted to produce other t>pes of 
particles such as amorphous silica, silicon carbide, iron ca'*oide and 
iron silicate oy the use of suitable gaseous precirsers. These include 
acetylene, i>'on caroonyl , silane and water vapor. 5uch a system 
snould be able to establish a significant body of data fo** tne optical 
p'^operties of small amorphous grains of astrophysical interest. 
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ChaptT b Tht Tf iiTatur» Utptndtni Nucitation of SiO 


Th« nucieatlon of refractory fflaterlals: netais, metal oxides, 
carbon ana its compounus is a problem of considerable Importance in 
astrophysics tor which no suitable framework c experimental data is 
yet available. Host nucieation experiments involve either simple 

1 j -3 a 

metals * or relatively low temperature organic compounds.''* Few 
experimental studies have been performed involving more complex high 
temperature species. * Yet recondensation of such species from 
meteoritiu materials could affect ttie earth's ozone balance^ or be a 
source of particulate pollution from the utilization of low grade coal 
resources, ae have tnerefore initiated a systematic study of the 
condensation of highly refractory species under carefully controlled 
laboratory conditions. Ttiese studies were initiated primarily to 
understand the formation and physical properties of the i.iterstellar 
dust 


Section II contains a rompiete description of tne apparatus and 
procedure used in this work. Section III presents typical data 
obtained in the system as well as the methods by which such data can 
be analyzed in order to determine values of Yi surface free energy 
of the condensate or J, ttie particle flux observed in the system. 
Section IV presents a comparison of the results obtained in our new 
system with those from the literature. In particular, we compare the 
value of Y Obtained by our analysis with that which might be expected 
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of SIO. (X K5). also dlsouas tha applioabillty of an analysis 
basad on tha usual thaorlas of hoao^anaoua nuolaatloii for a coapound 
wtiioh doas not oondansa to a solid which ia thansodynaaioally stabla 
undar tha conditions at which it foras. SiU condansas to SiO^ rathar 
than to tha pradictad SiO^ in agraaaant with tha work of Day and 
Donn.^ Wa also compare our rasults with tha racant work of ^taphans 
and bauar^ on tha condansation of SiO at tha much higher tamparaturas 
obsarvad in thair shock tuba-lasar scattaring system. Wa find that 
our rasults, obtained at YOO K < T < 1000 K, extrapolate fairly wall 
to those obtained in thair system at K < T < Z2U0 K. Section V 

presents the conclusions wtiioh wa believe can be drawn from this 
study. 

11. EXPtKlHbNTAL PRUCEDUHE AND EQUIPHENT 

A schematic representation of the apparatus used in this study is 
shown in Figure 1. SiO is evaporated from a resxstively heated 
ceramic crucible 1/2" inside diameter, 1-1/2" deep. Tlie temperature 
of the crucible is measured by a Pt/Pt lit Kd thermocouple encased 
within a steel sheath and further shielded by a ceramic tube which is 
placed within the top of the crucible, because it was found that the 
SiU tended to form a cap over the crucible which then served as the 
source of SiO vapor seen in the system, tha thermocouple is situated 
so as to measure the temperature of the cap rathar than the 
temperature of the interior of the crucible. It was found that this 
temperature could be kept within ib** K for extended periods even when 
the temperature of th«: ambient gas varied significantly. 
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ScriMatlc Diagram of Experlaantal Apparatus (saa taxt 
for datallad dascrlptlon) . 


loo 


I50i STAINLESS STEEL 
BELL JAR 




The SIO evaporator la enoloaed within a 12" high, inaide 
diaaeter oylindrioai, resistlveiy heated muiiite furnace. The furnace 
ia used tc control the aabient teaperature at which condensation 
occurs. At right angles around the equator are four 1" inside 
diaaeter d" long alumina tubes. The temperature profile of the 
furnace has been determined as a function of evaporator temperature 
and total pressure. The primary uncertainty in our experiment is the 
effect of the evaporator on the temperature and position at which 
nucleation occurs. This uncertainty increases at higher evaporator 
temperatures and lower total pressures since our observed temperature 
gradient increases under these circumstances. Our overall estimate of 
tne uncertainty in the temperature at which particle formation occurs 
is i 20 " K. The range in ambient temperature within which 
measurements were possiOle lay between about 70U and lUUU'' K. 

The furnace Itself is enclosed within a 1501 stainless steel bell 
jar equipped with U ports, each of which is aligned with one of the 
tubes leading into the furnace. The system is pumped via mechanical 
and liquid nitrogen trapped diffusion pumps. We attain pressures in 
the 10 ^ torr range even when the furnace is baked at temperatures 
near 750 '' K. The system is typically baked out for about lb hours 
after it has been opened to the atmosphere before condensation 
experiments are begun. The total leak rate of the vacuum system is 
less than 25 um/hour. uur experiments are typically carried out at 
total pressures in the range of 2 to 50 torr. 


A hl4h pr ssurt X«non arc laap, aountad on tha ball Jar, la 

baaaad through a HgF^ window and through ona of tha aluaina tubaa. 

Ubaarvationa at 90* to thia baaa giva qualltatlva data on tha 

condenaata auoh aa color and dagraa or polar! aatlon of tha acattarad 

light, wa can alao gat aoaa Idaa of tha spatial aitant of tha 

particla formation in this way. A 0.3 matar monochromatar is mountad 

in lina with tha baau on tha opposita sida of tha bell jar. 

Extinction maasuramants at 20U nm are used tu determine the onset ot 

condensation. In addition, the extinction spectrum of the condensates 

u 

can be determined in situ from 19;> nm to gUO nu. 

The vapor pressure of HiO as a function of temperature was 
determined thermodynamically. The experiments were done under 
reducing conditions (in H.. gas) and P (SiO) is therefore given by 
equation (1) over the evaporator . A simple model of the 

P (iiiO) * (760 torr) exp[li».«l (1 . 
eq T 

diffusion of the SiO vapor from tne crucible to the region in which 
condensation is observed to occur, will reduce the partial pressure of 
SiO in this region by about a factor of 10. If convection is 
important this pressure could be further reduced. 

Our measurements were obtained as follows. First, the system was 
cleaned and ttie evaporator filled witn union Carbide Select (^rade 
silicon monoxide. The system was then evacuated and pumped overnight 
with a diftusion pump while the furnace remained at a tem(.>erature of 
approximately j00''C. before a run, liquid was added to the trap 


and tha systeia was punped for at least two more hours. The pusp was 
then Isolated fron the system and gas introauced. After the 
temperature of the ambient gas reached equilibrium, the temperature of 
the bio evaporator was Increased until condensation was observed. 

This was indicated by a decrease in the intensity of light transmitted 
through the furnace and by the visual observation of scattered light 
at yo'* to the beam of the Xe arc lamp. 

We now increased the temperature of the ambient gas until the 
scattering (or extinction) disappeared. At this point we reinitiated 
condensation either by lowering the temperature of the ambient gas or 
by raising the temperature of the biO evaporator. This cycle was 
repeated nunerous times — usually until the SiO evaporator was 
emptied. 

We usually found that an increase (decrease) in the temperature 
of the ambient gas tended to increase (decrease) the temperature of 
the SiO crucible and vice versa. This change could be compensated for 
by a slight decrease (increase) in the voltage across the crucible. 
Although care was taken to hold the temperature of the crucible 
constant during a change in the furnace temperature, tne temperature 
of the ambient gas was allowed to float freely as the temperature of 
the crucible was varied, because of this, the temperature of the 
furnace generally increased during a run and yielded data at several 
condensation temperatures (see Figure 2). 


Figure 2 . Typicei Date vs, Tiae . Shown ere a) the current fro« the 
photoaultiplier tube, b) the anbient temperature of the gas, (T > 
450*'C) , and c) th<^ teaperatuare of the SiO evaporator (T > 925**C) . 
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50 torr 




Th«raocoupl« tMp«ratur«a war# aonltorad via atrip chart raoordai 
and the data read off at one minute intervale after ooapletlon of a 
run. The PM current waa aimilarly treated uaing a aeparate, 
aynchronized recorder. The oonditiona under which condenaation began 
(ended) were uaually determined viaualiy by the onset (disappearance) 
of scattering during a run. the final data analysis was made more 
self-conii latent by the use of the PH current which is a direct, 
quantitative measurement of the particle extinction in the system. 

Wc arbitrarily defined the c;;set of particle formation to be the 
period at which extinction was observed. It should be noted that 
the transition from Oi extinction to extinction usually occurred 
on a time scale of only about 10-15 seconds, because the tempe'*atures 
of both the biO evaporator and the ambient gas cnanged more slowly 
than did the measured extinction whenever avalanche nucleation 
occurred, tnese temperatures were quite well determined and changed 
only a total of 10-^0 degrees between the time just prior to 
nucleation and the point at which significant extinction (^10:1) was 
observed . 

111. KEbUL'^b AND AKALYblb 
a. first order Analysis 

Figure 2 shows typical data collected during one experimental 
run at 50 torr. bhown are the furnace temperature, the evaporator 


teaperature and the current from the photoauitiplier (PH) tuoe as a 
function of time. One can Inaeoiateiy notice that the PM current 
gradually decreases with time. We believe that this is due to a thin 
film of condensate which accuoulates on the HgF^ windows during the 
course of the experiment. One also notes that the PM current varies 
approximately inversely with evaporator temperature and directly with 
ambient furnace temperature. Since the current from the PH tube is 
proportional to tne transmitted light at 200 nm. as the particle 
concentration increases, the PM current decreases. As expected, an 

increase in the evaporator temperature increases the partial pressure 
of SiO and therefore the rate of particle formation. This, in turn, 
increases the extinction. Conversely, a decrease in the furnace 
temperature causes an increase in the extinction, since such a 
decrease increases the supcrsaturatiun of the vapor and therefore also 
increases the rate of particle formation. 

/n estimate of tne particle density can be obtained from 
extinction measurements using the relation 

I = Iq exp(- N« a^i (2) 

where I is the transmitted intensity, is the incident intensity, N 
is tne average particle density, a is the radius oi the particles, i 
is the path length and the extinction efficiency, is given by 

equation (3) for small particles. 
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Th« lignt soattvrtd lYoa thu partioivs at right anglaa to tha 
InoiOant bta« waa aatraatly biua ana ^1001 poiariiaa, Ma aatiaata 
that tha avaraga partiola rauiua (a) iiaa batwaan 20-40 na (aaa 
Karkar, 19bV). Ilia pathiangtii ovar which axtinotion oocurraO waa 
aaiooa lass than 2 ca and probably not Mora than b oa in langth. Our 
axtinotion aaasurananta wara aada at i a 200 ma. If wa uaa tiia 
optical oonatanta for at 200 nm a 2.6; > O.U^^ and 

ohoosa it axtinotion as tha point balow which wa dalina particla 
formation to caasa, wa find that log N a <j.7 t 1.7. 


Tha tima scala ovar which nuclaation occurs is not wall definao. 

A lowar limit to tha nuclaation rata can ba obtainad by obsarving tha 
time naadad for a stabla cloud of pcrticias to disparsa altar tha :>iu 
evaporator is turned off. In a typical axparimant, a cloud causing 
12t axtinotion disappears i< It extinction) about b seconds after 
power to the i>i0 evaporator is shut olf. Ttiis yields a minimum 

u • 4 _ 1 

formation rata ol about 10 particles cm s 


A maximua rata can be estimated from kinetic argmants in tha 
following manner, first, assume that biU laolvculas stick on every 
collision and grow to ^20 nm radius particles. Lacli of these 
particles contains ^10^ molecules If we assume p ^ \ for a fluffy 
particle, for an partial pressure on the order cl 10 utm and an 
aiablant gas twmparaturc ol approximately yoo K, the blO-SlO collision 




trcquanoy la on tha ordar of 10^ a*\ Tharafora a typical partlola 
foraa in laaa thwn a. OOaarvationa indicate particle foraation 

(aore than IS extinction) can occur laaa than 1 aacond after a 
dacraaae in furnace taiaparature. Thia supporta uur aatiaata and 
iapliaa a foraation rate of at laaat 10^^ partioiaa oa*^ a'\ Thia 
calculation ignorea the fact that the coiliaion oroaa section for a 
growing cluster increases as the 2/3 power of the number of aoleoules 
already incorporated within, as well as the fact tliat particles aay 
only nucleate into clusters containing as tew as 10 -10'' molecules 
before coagulation becomes an important growth process. If we include 
such factors our nucleation rate could easily be as high as 1U^^-10^^ 

.4.1 2 i 

particles cm *' s since each 20 nm particle could represent 10 -10 

i 2 

clusters each containing 10 -10 luonomers. 

The supersaturation was calculated in the following manner, 
hirst, we assume that the temperature of the SiU cap determines P(SiO) 
in the region just above the crucible and that this pressure can be 
determined using equation (1). Next, we assune that the SiO 
concentration gradient above the crucible can be modeled as a simple 
diffusion process, hinally, we assuae that the equilibrium 
partial pressure of SiO at ttte point at which nucleation is observed 
to occur is determined by the temperature of the ambient gas through 
equation (1). Nucleation always occurred approximately 1/2 to 1 inch 
above the crucible. Tlierefore, as an example, if the crucible 
temperature was 1j21>K and nucleation occurred at an aiabient 

I) 

temperature of b7bK, the supersaturation would be 10 . 
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Figure 3 is • plot of the calculated supei saturation (S) of SiO 
vs ambient temperature (T) at times when extinction is just observable 
(between 0-1ik). Data obtained wnen extinction is "turned off" due to 
an increase in ambient temperature is indistinguishable from that 
obtained when it is "turned on" due to a decrease in the ambient 
temperature. The experimental data from which the points in Figure 3 
were calculated was obtained at total pressures of 50« 35 and 20 torr. 
Although the data obtained at 35 and 50 torr is quite similar, the 
points obtained at 20 torr seem uniformly high. This could be the 
result of a nimber of factors. 

First, it is possible that diffusion at 20 torr could reduce the 
partial pressure of SiO in tne region of nucleation by as much as a 
factor of X).75 as compared to data obtained at 35 or 5u torr. This 
small correction would suffice to lower all 20 torr data taken at 
temperatures less than 925 K into agreement with higher pressure 
results. Second, it should be noted that even the 50 torr data at 950 
K is slightly higher than expected from a simple extrapolation of data 
taken at lower temperatures. This sut^ests the possiblity that as the 
crucible temperature increases, the region of nucleation moves away 
from the crucible. This would mean that the measured ambient 
temperatures in the region of nucleation could be too high by as much 
as 10-15 K while ihe estimated pressures for SiO could be high by as 
much as a factor of 6 due to the additional diffusion of the monomer 
away from tne crucible before nucleation occurs. From previous 
temperature calibration measurements, we know ^:hat the temperature 
gradient away from the crucible increases as the temperature of tne 


Figure j. Supcreaturation va, Awbient Tenperature . Shown is data 
taken at 3U torr (filled circles), 35 torr (crosses) and 20 torr (open 
circles)) fcr SiO evaporator temperatures of 1325** K (region A), 1473** 
K (region b) , 1443<> K (region C) and less tnan 1370<* K (all other 
data) . These represent points at which condensation could Just be 
noted (j 10 ' ' - s”^). 
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crucible itself increases and as the total pressure decreases. 
Therefore, the uncertainty in the exact region of nuoleation could 
account for the observed discrepancy of the low pressure, higher 
temperature data. This explanation seems even more reasonable when 
one realizes that the data in region "A" of Figure i was taken at a 
crucible temperature of approximately 13<?5'*K, that in region **B" at 
147^**K, that in region "C** at and all other uata at a crucible 

temperature below It should be noted triat the slope of the 

data in region "A" is approximately the same as that in regions "U" 
ana "c" and also the same as the slope of the low temperature portion 
of the curve. These regions are merely displaced to higher T and/or 
B as the temperature of the crucible Increases. 

b. Analysis via Nucleation Theory 

Nucleation theory predicts that the particle formation rate 
(J) should be related to S ana T by the following relationship^^ 

2y - 16 I 

J = — V exp (M) 

ira 3k^ T^(in:>)^ 

where y is the surface free energy, m is the mass of the condensing 
molecule, V is the voline of the condensing molecule, is the nunoer 
density of the monomer ana k is boltxman's constant. If the variation 
in is small, it is easy to see that for constant J, an increase in 
T necessitates a corresponding decrease in 5. 
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Equation (4) can be rearranged to the following expreaslon 


i 

I 
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I 

I ' 
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♦ m (- 


T'^(inS) 


2 




1/2 


) s in (N^^) 


(5) 


where * t(2V^)/(*«) J and s 16 w A plot of 

(T^(inS)^)~^ vs in(N^^) s^iould be a straight line if y remains 

approximately constant over the temperature range of the experiment. 

An average value for y can be obtained from the slope of this line if 

V can be estimated. This in turn allows one to calculate J 

2 

independently f'om the value of the in(N^ ) Intercept. Figure 4 is a 

plot of (T^( in!i)^)~^ vs in(N^^) in which the calculated 

super saturations have been corrected empirically for the temperature 

of th'i crucible. As can be seen, the slopes of the three sets of data 

obtained by a linear regression analysis lie between 4.7 x 10^^ and 

5.2 X 10^^ with an average slope of 4.9 x 10^ \ If we use a value of 

2.4 X 10“^^ cm^ for V, this implies a value of 500 ergs/cm^ for the 

surface free energy of the condensate. If we then use this value for 

Y« assuoe a value of 7.3 x 10 g for m and use the value of the 
2 

overall in(N^ ) intercept, given in Figure 4, we calculate J to be 
about 5 X 10^^ particles/ cm^. This is in reasonable agreement with 
our previous estimate for J from the optical extinction measurements 
discussed above. 
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I 


Flaure 4 . (T^ (^ns)^ )"^ vs, in(ll ^)^. Shown are the siopa, Intaroept 

and correlation coefficient for runs at 30 torr (filled circles), 33 
torr (open circles) , 20 torr (crosses) and the combined data as 
obtained from a linear regression analysis. Also shown is the 
resulting best fit through all of the data. 
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1tl 

A ■odll'loation to baokor-Doring thaory praOlota t^.at tha 
nuolaation rata darivad t'roa claaaical thaory should ba aodit'iad by an 
additional faotor of X/b such that 

3 2 

X 2y V X N , 16 I 

jt , _ . j , — L. exp C — -) (b) 

i> la S 3 k^ T-^CinS)^ 

whara X is a factor which accounts for statistical aachanical 
contributions to tha free anergy of foraation of the clustar due to 
previously ignored vibrational and rotational uegrees of freedoa. If 
we assume that X is approximately constant over the taaperature range 
of the expariaent (725-950K). then a plot of CT^(xnb)^)”^ vs xn(k^^/S) 
will again yield a straight line, the slope of which is equal to 
The xn(k^ /b) intercept, however, now equals LtnCJ/^^ ) - |nXj. 

Hecalcuiation ol the data plotted in Figure <i according to equation 

2 

(6) yields a value of 6b0 er^s/cro for y> If we assume that our 

q 1b 

(tarticle production rate must lie between about 1U end 1U particles 
-3 -1 

ca a as discussed previously, then we find that ll.b < tnX < 26> 
IV. U16CUSS1UK 

We nave derived values for the surface free energy of our 

3 2-1 2 

condensatenominally bi^U^-from plots of (l'’(inS) ) vs either tn(N^ ) 

2 2 

or /S)) of bUU and b*>U ergs/cm respectively, tariier 

calculations ot y based on extrapolation of measured values for 
silicious slags to concentrations of 1UU> SiO^ yield values in the 
region 1bU-27i ergs/cm^. a value of 269 ergs/ca‘ was derived 
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17 

froa ■••sur«atntt of th« h«at of aoiutlon of Morphous filioio aoid 
whilt eillot, GloiMr and Raaakrlahna^^ glva a valut of 350 arga/oa^ 
for liquid ailioataa in ganaral. Blandar and Kats^^ aatiaate tha 

■j 20 

aurfact fraa anargy of aolid ailioataa to ba 650 arga/oa . Bruce 
givaa tha following foraula for tha avaraga aurfaca fraa anargy of 
g-criatobalita aa a function of teaparatura 

Y « (925 - 0.193 T) arga/ca^. (7) 

It anould ba noted that y thia aliotropa of SiO^ variea only 

20 

alowly with teaparatura. truce expacta y for tridyaita, which haa a 
lower denaity, to be lower and y for quartz to be higher. Since it ia 
probable that our aaorphoua condenaate haa a lower denaity than any 
cryatalline form of SIO^. ita aurface free energy ahould be atill 
lower. It ahould be noted that Hhillip^^ feela that the optical 
propertiea of compounoa of the type SiO^ (0 < z < 2) vary continuously 
throughout the series. If this applies to other properties as well. 

then the value of y for Si,.0.. should be between that for silica (> 300 

2 2 
erg a/ cm ) ana that of silicon metal 726 ergs/cm ). Therefore, the 

2 

values which we derive (,r50U < y < v*b50 er^s/cm ) seem quite 
reasonable. 

In spite of the fact that we derive reasonable values for both y 
and J using techniques based on nucleation theory, one must still 
question whether such results could ue fortuitous. The fundamental 
question in this system is, "hhat does 6 mean?" We evaporate 6iU ana 
nucleate Ttus has previously been noted by Day and Bonn and 


it oonflnied in th«s« sanpltts by tht presence of e ohsraoterlstio 

21 

f«atur« in the infrar«sd aoaorption spec trial at 11.5 pn • Although 

22 

the structure of has not been definitely determined , 

aeasurement of the Si concentration in our sample via x-ray 

fluorescence confirms this stoichiometry. One possible structure that 

has been suggested for this compound is that of a typical SiU^ 

22 

tetrahedral network in which every fourth 0 a cm is missing . 

If the sample is annealed under vacuin at 75u K for periods on 

the order of 10-20 hrs the 11.;> sa peak disappears and is replaced by 

a peak at 12.5 mbi which is characteristic of amorphous silica. It ia 

not unusual for samples collected in our system to show both 11.5 pa 

and 12.5 absorption bands. It is known that SiO is the primary 

silicon bearing species in the vapor when either SiO or iiio^ is 

vaporized under reducing conditions. Nevertheless, supersaturation 

is defined in terms of the partial pressure of a gas phase species in 

equilibrium with a solid phase. If the solid phase wtiich actually 

nucleates is itself unstable at the temperatures at which nucleation 

occurs, then a thermodynamic formulation of the problem could be in 

serious error. Under such conditions, nucleation might easily be 

controlled by such kinetic considerations as which poljAorph first 

forms a metastable phase even if this phase is not ultimately the most 

23 24 

thermodynamically stable under the conditions of the experiment. * 

25 

In fact, Ostwald concluded that the most stable phase will generally 
not be the first to nucleate. 
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A furtbtr problM arises if one uses our values of Vt T and S in 

order to calculate n^, the nuaber of aonoaers in a "critical" sized 

c 

cluster. This is given by standard Becker-Doring theory as 

n^j ■ 32 1 V^/3(KT tnS)^. (b) 

Typical values of n derived froa our data indicate that the critical 

c 

sized cluster at a temperature of 750K is the aonomer. At 900K tne 
critical nucleus becomes the diaer while at a teaperature in excess of 
930K n^ is the triaer. At such aaall values of n^ the applicibility 
of the liquid drop approxiaation upon which classical nucleation 
theory is based becomes extremely doubtful. This in turn leads one to 
doubt the validity of a calculation of y based upon nucleation theory, 
even if the value derived seems reasonable, since nucleation theory 
itself may not be applicable to this sytem. 

In light of the above arguaerts, it might seem more appropriate 

to report the results ol nucleation experiments on very reiractory 

species in the manner adopted by Bauer atiu co-worwers. * ** ITiey 

present plots for log vs T^ of the conditions which yield avalanche 

nucleation. is defineo as the critical partial pressure of the 

condensing species (in torr) ana T^ is tne condensation temperature* . 

.Stephens and bauer^ have recently completeo a study of the nucleation 

of Bi(>^ using a shock tube at ambient temperatures between ^ }2'30 K 

and <r 4UU0 K. Figure b is a plot of log P vs T under conditions of 

c c 

avalanche nucleation which combines our data (700 K < T < 1000 K)with 
that of Stephens and jauer^ in the range 1250 K < T < 2200 K. Their 


Fi gure 5» L ok P^q vs. Aabient Teaperature . bnown is the Critical 
Pressure vs. Condensation Temperature at wnich avalanche 

nucleation was detected in this work (filled circles) and in the shock 
tube experiments of Stephens and bauer^ for oxygen to silicon ratios 
of 3.1 (squares), 1.6 (open circles) and 1.4 (triangles). Also shown 
is an exponential regression best fit to our data (correlation 
coefficient of .*)5) extrapolated to that of Stephens and Bauer 
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data was obtained at oaloulated 0/Si ratios of 3. 1, 1.6 and 1.4 
respectively and indicates a tendency for nucleatioti to occur at lower 
at a given teaperature as the ratio 0/Si increases. Also shown in 
Figure 4 is a linear regression fit to equation (9) 

m « 0 »nT - m a (9) 

wnich used only data obtained in this work (7UU K < T < 100U K) and 

was extrapolated to higher temperatures. We find that tn a s 70.52, b 

2 

* 9.66 with a correlation coefficient (r ) of 0.95. The agreement of 
our work with that of Stephens and Bauer seems quite good. This is in 
spite of the uncertainty in the exact polymorph of silica which 
actually condenses in either system. Stephens and Bauer report 
finding amorphous SiO^ in high 0/Si experiments and a combination of 
Si metal and amorphous SiO^ in low 0/Si runs.^ 

V. CONCLUSIONS 

We have demonstrated that controlled, icproducible nuclcation 
experiments involving refractory oxides can be performed in the 
apparatus described in this report and that the nucleation of SiO is 
extremely sensitive to the temperature of the ambient gas. We have 
found that the vapor phase nucleation of SiO in yi<ilds the solid 
phase SiO^ (where X I.5) rather than the thermodynamically favored 
product SiU^* This is in agreement with the earlier work of Day and 
Oonn.*' 
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Analysis of our data within the fTaaework of the Becker-Dorlng or 

filander-Katz theories of hooogeneous nucleation yield values for the 

2 

surface free energy, y, of our condensate of 500 ergs/ cm and 650 
2 

ergs/ cm respectively. This is within the range which could be 
expected of a silicate compound, the composition of which lies between 
that of amorphous silica and silicon metal. However, because of the 
difficulty in the formulation of a realistic definition of the 
supersaturation ratio, S, for a material which preferentially 
nucleates into a form which itself is thermodynamically unstable, 
standard theories of homogeneous nucleation which rely on 
thermodynamic argunents may not be applicable. In such a case, a 
kinetic theory of nucleation, such as that proposed by Bauer and 
Frurip,^^ may be the preferential way to describe the condensation 
process . 

Finally, we find excellent agreement between our results, 
obtaineu at relatively low temperatures and high supersaturations, and 
those of Stephens and bauer^ obtained at much higher temperatures. 

Both experiments den<onstrate that the vapor phase nucleation of SiO 
requires a very high degree of super saturation. 
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Chapter 6 


The Wucleatlon and Infrared Spectra of Aaorphous Mg-SiO Saokea 

binary nucieatlon processes are expected to be of considerable 
importance in the chemically complex astrophysical environments in 
which interstellar grains are likely to form (Donn, 1^76; Donn et al, 
I9b1). Unfortunately, theoretical treatment of this problem is much 
more difficult than for the case of homogeneous nucleation. 
Consequently, few theories of binary nucleation exist (Wilemski, 1975 
a,b; Keiss, 1950; Hirschfelder , 1974) and these have usually been 
developed to treat nucleation in relatively low temperature systems. 

Since most theories of binary nucleation build upon the classical 
theory of honogeneous nucleation, all of the problems inherent in this 
framework are intensified. This is especially true for the vapor 
phase nucleation of highly refractory species (Stephens and Bauer, 
1981; Nuth and Donn, 1981a). It is interesting to note that many of 
the treatments of the formation of grains which have appeared in the 
astronomical literature (Draine, 1981; Yamamoto and Hasegawa, 1977; 
DeGuchi, I960) have been based upon homogeneous nucleation theory. 

Few quantitative studies of the binary nucleation of refractory 
materials have been published (Stephens and bauer, 1961) althoi^h the 
work of Day and Donn (1976) indicates that considerable 
supersaturation is necessary before condensation occurs in such 
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systems. In a pravious study (Nuth and Donn, 1981a; hareaftar Papar 
1) wa raportad maasuramants of tha vapor phasa nucieatlon df SiO in 
gas as a function of tamperatura for tha ranga 750K < T < 1000K. kia 
report hare measurements of the binary nuoieation from the vapor, of a 
Hg - SLO - system in the range 750K < T < 1010K. 

Section 11 presents a brief description of the apparatus used in 
this study as well as an outline of our experimental procedure. 

Section 111 presents the experimental data for the nucleation of the 
Hg -S 1 O-H 2 system as compared to that for the "homogeneous" nucleation 
of 810-^2 (Paper 1). Section IV presents the infrared spectra of 
representative samples of condensate both as collected and after 
annealing at 10U0K and l^iSOK. Section V discusses the relevance of 
this study to the condensation of refractory grains in various 
astrophysical environments. 

11 Experimental Apparatus and Procedure 

llie experimental apparatus used in this study was similar to that 
reported in Paper 1 with several relatively minor modifications. Hg 
was introduced into the furnace system by placing a crucible 
containing the metal (and a tnennocouple) within the furnace. Ihe 
temperature of the ambient gas the*’efore controlled the partial 
pressure of within the furnace. Several methods to independently 
vary the Hg pressure were tried, but all failed to produce a unifonnly 
mixed vapor and reproducible results (Mork to overcome this limitation 
is in progress). 
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The Hi partial pressure within the system osuld be msnltsred 
Independently by measuring the strength sf the Hi 265*2 nm resonance 
line In absorption, because the line was saturated under most of the 
conditions of the experiment only relatively qualitative Information 
could be obtained. This did show however that the partial pressure of 
Hg within the fu^^nace was reproduclbly controlled by the furnace 
temperature. The actual hK pressure was calculated empirically using 
the data In the Handbook of Chemistry and Physics (53nl Edition). 

This is plotted In Figure 1. As discussed In Paper 1, the partial 
pressu**e of SiO In the region of condensation was determined t>om the 
thermodynamic relationship given by Schick (I960) combined with a 
suitable model for the expansion of the SiO vapor from the crucible. 

The experimental procedure used in this study differed in two 
major respects from that reported In Paper 1. First, in ou*' previous 
study, the temperature of the SIO evaporator was held approximately 
constant while the temperature of the ambient gas was varied to 
initiate or halt condensation. In this study, because of the 
dependence of the Hi partial pressure on the temperatu'^e of the 
wnblent gas, we attempted to hold the furnace temperature constant and 
va*'ied the temperature ol the SiO crucible. This procedure tends to 
increase the error in the measurement of the partial pressure of SiO 
necessary for nucleatlon since it increases the uncertainty in the 
temperature of the SiO crucible at the time avalanche nucleatlon 


begins 


Figure 1. The partial pressure sf m^nesiun as a funotisn of 


temperature. 

r 
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S«oond, b«oauM w« us«d tht X« Arc lap - Monaotraatar systa to 
Bsnitor tha oanoantration within tha furnaoa, tha tlaa at which 
nuolaatian bagan ar andad ws datarainnd by aya fra tha praanoa ar 
absanoa af light aattarad at 90 ^ ta tha baa. Aa in tha praviaua 
study, light aattarad fra thase grains was axtraaly blue and 100} 
polarizad. fcxpariaent shawad that far tha Hg-SiO systa, this aathad 
was, an tha average, as sensitive ta the onset af candensatian as ms 
nanitaring tha extinctian af the newly farnad particles at 200 na. 
Unfortunately, this technique certainly was neithe'' as precise nar as 
quantitative as that used in tha p**aviaus study (Paper 1 ). All runs 
were dane at a total pressu’*e af 35 tarr af H2. 

Ill Results 

Data collected in this study was analyzed using the same 
correction factors determined fran our previous wark (see the 
discussion in section 111 af Paper 1 ). Figure 2 is a plat af the 
partial pressure of SiO at which avai anche nucleatian occurs as a 
function af ambient temperature far the candensatian af "pu’*e'' HiO 
( 2 a) and far the Hg-SiO system ( 2 b). 

Similarly, Figure 3 is a plat af the supersaturatxan af SiO (with 
respect ta SiO solid) necessary ta initiate avalanche nucleatian far 
the SiO-H^ system (Figire 3a) and far the Hg-SiO-H^ system (Figure 
3b). In bath af these plats, it is obvious that although the presehce 
at Hg lowers the required SiO pa*'tiai pressure ta induce avalanche 
nucleatian at temperatures less than ^ 900-950K, it appeal's that is 
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Figure 2a. The partial pressure of SIO at which avalanche nucleati^n 
is observed in the SiCMl^ systea as a function of the teaperature of 


the Kibient gas 


Ql- 



JURE (K) 


Figure 2b. The partial pressure sf SiO at which avalanche nucleation 
is sbservad in the MS-SiO-H^ systea as a function of the t^peratur'.* 


of the aabient gas 
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Figure 3a. 
nuoieatisn 


The supersaturation sf SLO vapsr at which avalanche 
ia abserved in the SiO-M^ system as a functisn sf the 


temperature sf the ambient gas 


ln(S) = 34.4-0.0287T 
( r2 = 0.98 ) 
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Figure 3b. The super saturation of SiO vapor at which avalanche 
nuoleation is observed in the Hg-SiO-H^ system as a function of the 
temperature of the ambient gas. It should be noted that the 
supersaturation of this system with respect to a mineral such as 
olivine is likely to be much higher. 
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of little, if any, impsrtance for the onset of condensation at higher 
teaperatures. This is in spite of the fact that much nore is 
present at temperatures in excess of 950K in our system, than at the 
lower tefflpe<*atures whe'*e it is effective (see Figure 1). The 
implications of this finding will be discussed in more detail in 
section V. 


IV Infrared Spectra of Typical Condensates 

The infrared spectrin of an average of the material collected 
from three separate nucleation experiments is presented in Figure 4. 
Also shown are the spectra of samples of this material annealed under 
vacuua for 1, 2, 4, 8, 16.5 and 30 noun's at 1000K. All spectra were 
obtained by dispersing the sample in finely ground Kbr powder and then 
pressing the mixture at 10,000 psi into a clear disk. 

Figure 5 shows the spectrin of the material collects, when pure 
SiO vapor was condensed in gas (Paper I). Also shown in Figure 5 
IS the specfiu of a sample of annealed at 1250K for 30 

minutes. It should be noted that the condensation of SiO vapor in 3b 
torr of gas at temperatures in excess of 750K produced a mixtu'^e of 
amorphous SiO^ and ^> 12^3 d**ains (Pape^* 1). This is obvious from the 
spect^'in of the material since it contains peaks at both 11.4 and 12.5 
micons as well as the more familiar peaks near iu and 20 microns. 

The peak at 11.4 microns (attributed to Si^O^) disappea>'s when the 


Figure 4. The Infrared spectrin of an averi^e Mg-SiO moke, as well as 
the spectra of anoke samples annealed in vacuo at 1000K for 1, 2. 4. 

6. 16.5 and 30 hours. 

A. Unannealed amorphous magnesiiu silicate anoke 

B. ^oke annealed at 1000K for 1 hour 

C. ^teoke annealed at 1000K for 2 hours 

D. ^oke annealed at 1000K for 4 hours 

£. Smoke annealed at 1000K for B hours 

F. !:teoke annealed at 1000K for 16.5 hours 
U. Smoke annealed at 1U0UK for 30 hours 
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Figure 5. The infr .red speotrua of anorphous silica species. 

A. The spect**ui of a sample of amorphous quartz prepared by annealing 
Si^O^ SBoke in vacuo at a temperature of 1250K for 30 minutes. 

B. The spectrvm of a smoke produced in the SiO-H^ system at a 
condensation temperature of 750K which shows both the 12.5 micron 
feature characteristic of amorphous quartz and the 11.2 micron 
feature associated with Si^O^. 

C. The spect«*ua of 'pure' ^^2^3 smoke preparea at room temperature. 


MICRONS 



1 


1400 1200 1000 800 600 

WAVENUMBERS 




saaplt is annealed In vaouiai for ^ alnutes at 1250K. This is 
indioative of the t^ansforaatlon of the into aaorphous quanta. 

V Discussion 

Our results will be discussed in terms of both their implications 
for nucleation theory in general and their relevance to condensation 
processes in the circunstellar shells of cool stars. 

a. Nucleation Theory 

Paper 1 emphasized that, although the data for pure SiO-M^ 
nucleation could be analy :ed in terms of Qassical Nucleation Theory, 
this analysis led to some logical inconsistencies. For instance, 
using the value for the surface free energy of the (^0)^ condensate, 
one calculates that the monomer, dimer and trime** constitute the 
critical sized cluster at 750K, 900K, and 950K respectively. An 
additional problem stems from the fact that the condensate in this 
system is found to be Si^O^ rather than the thermodynamically stable 
species SiO^. This causes a problem in that a logical and physically 
meaningful definition of supe** saturation is difficult to formulate 
Since the usual definition involves the free ene^-gy of formation of 
the stable species, while the application refers to the cluste''s which 
actually nucleate. 

Because theories of binary nucleation are generally only 
extentions of Qassical Homogeneous Nucleation Theory, we might expect 
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•lailir probiws ts ««rg« if suoh th« 9 rl«t w«r« applied to our 
rtsults. Unfsrtunataly, additional problaia arisa sinoa thaoriaa of 
binary nuolaation ara raally fomulatad to account for tha 
prec pitatior of nearly ideal liquid aiKturaa froa vapors of two 
coapletaly aiscribie esnponants. In suoh a systasi it is typical that 
both species are at least slightly supersaturated . In our systaa, 
although SiU is always supa** saturated , the pressure of aagnasius never 
exceeds the equilibriua vapor pressure. In fact, the aagnesiiu 
concentration most likely remains below the equillbriun value due to 
delays inhe^’ent in the diffusion of the vapor from the source into the 
region of condensation. 

Furthemore, we actually condense an doorphous magnestun silicate 
(Day and Donn, rather than a mixture of magnesuitm metal and 

silica. In this sense, it would be more appropriate to use a theory 
of nucleation developed for systems in which the monomer fora of the 
condensate must first appear via chemical reaction in the vapor. Such 
a theory has been developed by katz and Donohue (1962). In order to 
apply this theoretical treatment to ou** data we need to know the 
surface tension of our clusters as well as various reaction rates for 
the formation of ou^ "monomers'* from a magnesiua rich, SiU vapor 
diluted in H^. Unfortunately we don't know the exact composition of 
ou** "critical" clusters. Furthermore, we really don't even know the 
exact composition of the important vapor phase monomers - i.e. Hg-SiU, 
Mg-(SiO)^, Hg-SiU-Hg, SiO-Mg, etc. -> which condense. Although we do 
plan to study this problem in the near futu«’e using a quadropole mass 
spect*‘oneter to monitor the p*‘e-nucleation vapor composition as a 
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funotlan af taapcraturt, prasaurt and initial ohMieal ooMponition, wa 
can now aiaply stata tha qualitativa iaplioationa af aur raauUs. 

Inspaction of Figuras 2 and 3 indicatas that tha prasenoa af 
magnasiiai vapor reduces the concentration of SiO needed to initiate 
avalanche nuclaation for temperaturas lasa than about 925K. This is 
to be aipected if Hg-SiO clusters are uaportant nucleating species. 

At taoperaturas in excess of 925K« such clusters apparently bacoae 
less stable than those involving only SiU, even though the 
concentration of m^nesim in tha vapor is 2 or 3 orders of 
■agnititude greater than that of SiO (see Figs. 1 and 2). Mucleation 
unde** such conditions is therefore initiated by the precipitation of 
(SiU)^ iron the vap.«r. 

Inspection oi the infrared spectura of ou** "nagnesiuD silicate" 
(Figure 4) shows suostantial differences from that of snoke 

(Figure 5) nucleated froui a SiO-h^ systen (Paper 1). This inoicates 
that although fomation of (^0)^ clusters aay initiate the nucleation 
process in our h^-SiO-H 2 systeo at high teaperatures, magnesius is 
apparently still incorporated into the grains during the g^'owth 
process. 

D. Astrophysical Inplications 

These results have two aajor implications with •'egard to tne 
formation of refractory grains in the regions surrounding cool stars. 
First, condensation in such sources may oe initiated by the formation 


3 f cluste<’s, althot^h the conpasitlon of the resultant grains 

would then be subject to fflodifications by subsequent reaction with 
other refractory constituents also p*’esent in the gas. Second, it may 
be possible to observe distinct chemical differences in the infrared 
spectra of grains newly formed in stellar sources, especially compared 
to those grains which may nave undergone significant reprocessing such 
as those in molecular clouds. Ihese implications will te discussed in 
more detail below. 

Grain Formation ; It has been shown that the vibrational population of 
SiO molecules in circuosteliar regions is characterized by a 
temperature significantly below the kinetic temperature of the ambient 
gas (Nuth and Donn, 19ti1b). It is well known that the rates of gas 
pnase reactions are very sensitive to the vibrational states of the 
reactants (Polanyi, 1973: bartozek et al, 1981). It is quite 
possible, therefore that nucleation in stellar sou*'ces could still be 
initiated by Hg-biO clusters if such clusters could be further 
stabilized by radiative cooling. Of course, there is the additional 
possibility that Fe-SiO or any nunber of combinations of the more 
abundant •'efractory elements might also be important sources of 
condensation nuclei. This possibility can however be checked by 
laboratory studies similar to those reported in this work. 

Infrared Observations ; Figure 6 is the spect”un of OH26.5 +0.6 from 
the work of Forr-est et al (1978) taken in Ap'^il/Hay 1973 and again 
hay/June 197b. beveral points should be noted. Fi’*st, in the 
spectrun obtained in April/Hay 1973, it is possible to identify minima 
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Flgur* 6. The lnfrare<( spectrm 9f OH 26.5 ♦ 0.6 obtained over a two 
year period after an Infrared brightening (Forrest et al.,19Y8) end 
which shows features at 11.2, 12.5, 18-19, 20 and 22-2'j aiorons 
suggestive of the features in Si 20 ^, amorphous quartz and partially 
annealed Hg-SiO aiokes. 
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at i>10, 11 aid 12.5 BicronSi Ihls could easily be attributed to a 
aaall population of Si^O^anorpnous Si02 grains (see Figure 5) 
conbined with a more ninerous population of partially annealed Hg-2iiO 
grains. There is even the suggestion of the peak at >r9.2 microns due 
to amorphous quanta, although this is much more speculative. 

A second point to note is that the spectrun of this object 
changes with time in a manner suggestive of the annealing of our 
Hg -SiO smoke (Figure 4). In othe<^ words, definite structure begins to 
appea** within the 10 micron band suggestive of the formation of a more 
C'*ystalline solid. The irregular scatter in this spectrun is not due 
to random noise but is real (Merrill, 1981) and due to the pointwise 
procedure used to obtain the measurements. 

The third interesting point is the similarity of the 16-25 micron 
spectrun of this object to that of our amorphous magnesian silicate. 
OH26.5 ♦ 0.6 may show weak features at <r18-l9, ^0 and <r22-23 microns. 
At least two of these components can be reproduced in the laboratory 
(Figures 4 and 5). Future experiments using additional condensible 
species (i.e. Fe or Ni) may allow other residuals within these rather 
broad features to be better identified. 

llie previous rema'*ks indicate that much more infonnation about 
the chemical nature of the grains in circunstellar sources can be 
Obtained from infrared observations than has hitherto been expected. 

In this regard, high resolution continuous spectra in the 8-1 3 and 
16-25 micron regions of grain forming circunstellar shells could be 
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•speoialiy interesting since such material represents new inte'*stellar 
grains which have undergone little reprocessing or annealing. 
Comparison with the chemical state of "average" interstellar material 
might reveal the process by which such fresh grains are annealed, 
homoginlzed and mixed to yield the normal population of "silicates" 
seen in dense clouds. 

Lven in such homogenized materials it should be noted that some 
st’*ucture within the 10 micron featu'-e is evident in the literature. 
For instance, the spectrun of the Tr«peziun shows a peak near 11 
microns (Forrest, Gillett and Stein, 1975) while that of the bN-KL 
complex shows a definite feature at 12.5 microns as well as possible 
st'*ucture within the 10 micron peak itself (Uillett and Forrest, 

1973). 

VI Conclusions 

We have demonstrated that the nucieution ba>’rier in the system 
^lg-SiO-b^ at temperatures less than <r925K is lower than that for 
SiU-H^ alone. At temperatures in excess of <r925K however, nucieation 
in the Hg-SiO-4i2 system is unaffected by the presence of magnesiun. 
This could have interesting applications to the problem of 
condensation in the atmospheres of cool stag's. 

We have also demonstrated an interesting correspondence between 
small residuals seen in the infrared spectra of various astronomical 
sources (b microns < i < 25 microns) and the spect'*a of specific 
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laboratory oondansad silicates. Ue suggest that more careful infrared 
observation of such sources could reveal valuable infonaation about 
the chenical state of the grains. Thjs could in turn lead to a better 
understanding of the processes which tranfona newly formed dust into 
the relatively hanogeneous material ooserved in the general 
interstellar mediun. 
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